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Chapter 1 


Vector Analysis 


1.1 Vector Operations 
The Dot product of two vectors is defined as 
A x B zABsin0Of 


where 0 is the angle they form when placed tail-to-tail. Note that A.B is itself a scalar (hence the 
alternative name scalar product). The dot product is commutative, 


A-B=B-A 


and distributive, 
A-(B+C)=A-B+A-C 


if two vectors are parallel, then A - B —AB. In particular, for any vector A, 
A-A =A? 


if A and B are perpendicular, then A-B =0. 


The cross product is defined as 
A x B — ABsin Oi, 


where fi is a unit vector pointing perpendicular to the plane of A and B. The cross product is distributive, 
A x (B—- C) 2 (A x B) - (A x C) 
but not commutative. In fact, 


(Bx A) = —(A xB) 


geometrically, |A x B| is the area of the parallelogram generated by A and B. If two vectors are parallel, 


their cross product is zero. In particular, 
AxA=0 


1.2 Vector Algebra: Triple Products 


(i) Scalar triple product: A-(B x C). Geometrically, |A-(B x C)| is the volume of the parallelepiped 
generated by A, B, and C, since |B x C| is the area of the base, and A cos Ois the altitude. Evidently, 


A-(Bx C)  B.(CxA) 2 C. (A x B) 
for they all correspond to the same figure. Not that “alphabetical” order is preserved. The “nonalpha- 


betical" triple products, 
A-(CxB)  B.(Ax C) - C. (BxA) 


have the opposite sign. In component form, 


N 


= 


w > 
z 
w > 
< 


A-(BxC) = 


N 
fs 
Ob > 


N 


(ii) Vector Triple Product: A x (B x C). The vector triple product can be simplified by the so called 
BAC-CAB rule, 


A x (B x C) =B(A - C) — C(A - B) 
1.3 Position, Displacement, and seperation vectors 


The location of a point in three dimensions can be described by listing it Cartesian coordinates (x, y,z) . 
The vector to that point from the origin is called the position vector: 


r=x+yŷ+zĉ 
its magnitude is given as 


r=Vxe+yt+2 


m m xR +y 4- zZ 
ro a/y +z 


is a unit vector pointing radially outward. The infinitesimal displacement vector, is given as 


is the distance to the origin, and 


dl — dx& + dyŷ + dzĉ 


the seperation vector from a source point to the field is given as 


its magnitude is 





è * r $ / . 
and a unit vector in the direction from r to r is 


/ 
Z r-r 
n = 


3 | 


= / 
r-r] 


2 


in Cartesian coordinates, 


Å = (xx) + (yy) + (z-z)? 
q — Vox -»? - (- 2? 
> (x—x")&- (y -y)$ - (z - z)£ 
n = 

eps 


v/(x x) - (y - y)? ( 


1.4 Problems with Solutions 


Problem # 1 
The Equation giving a family of ellipsoids is 


x2 y? z 


I d aea 
Find the unit vector normal to each point of the surface of these ellipsoids. 


Given a three-dimensional surface we can find what the normal vector is by using 


z Vu 
Uu = 
|| Vu || 


and knowing that 

v Ju, Ou. Ou, 2x. Yoy Z; 
u = x SS ek 
ox dy” 3z“ qm Tas 


and for the magnitude of this vector we find 


B py g 
Vu ||=2 
| ve ie T p^ * c^ 


thus the unit vector û is given as 


Problem # 2 

For each of the following vector fields A (1) determine whether A is constant, (2) find the divergence 
and the curl A, and (3) find the components of A in cartesian (expressed in x,y, and z), cylindrical (ex- 
pressed using p, and z) and spherical coordinates (expressed in 7,0, and 0). In all cases a,b, and c are 
constants. 

(a). A — a£4- b$ -- c£ 

A is constant because any point in the field will give you the same unit vector at that point. 

The Divergence of A is given by 


và 9A 9A | OA E 200 dE 
o Nax Oy dz) Ox Oy daz 


The curl of A is given by 


u 


VxA— 


a Ho w 
SYouxe 
6 Awa 


in cartesian coordinates the components of A are given by 


A,y=a Ay=b A,=c 


in cylindrical coordinates we must use a transformation matrix to go from one basis to another, i.e 
y 


6 cosó sind 0 & 
o |} ={ —sind cosd 0 y 
Z 0 0 1 Z 


where the 3x3 matrix 1s the transformation matrix that lets you represent cylindrical coordinates in 
cartesian coordinates. Thus we find 


6 = cosox+sinoy 
6 = —sinod&+cosoy 
= 4 
thus the components are given by 
Ap = -A-acosQ-Fbsino 
Ap = $6-A=-—asingd+bcoso 
Az, = Z:-A=c 


In spherical coordinates we must also use a transformation matrix. 


»^ 


sinOcosÓ sinOsind — cosO 
— | cos0cosó cosOsinó —sinO 
—singd cos Q 0 


© D =» 
N> 


where the 3x3 matrix is the transformation matrix that lets you represent the spherical coordinates in 
cartesian coordinates, this yields 


sin 8 cos ox + sin O sin þŷ + cos 02 
cos 8cos OX + cos Osin by — sin 0Z 
— sin OX + cos by 


=> 
|| 


© D 


thus we find the components to be given by 


A, — -A =asinĝcoso + bsinð sinp +ccosð 
Ag — 0-A =acosĝðcosp+ bcosəðsino -— csinð 
Ag = 9-A=—asing+ bcoso 


D D 


(b). A — ap -- bó -- c£ 

A is not constant because as your position changes, so do f and à. The 6 component is the distance 
from the z axis and the direction is always in the fdirection. 

The Divergence of A is given by 


] 2 lob ac 
V.A ——— t Ip 
pop P4 595. a: 


The curl of A is given by 
b 
VxA=-2 
p 
in cartesian coordinates we must use 
X cosó —sinó 0 6 
y |={ sind cosd 0 o 
2 0 0 1 e 


where the 3x3 matrix is the transformation matrix that lets you represent cylindrical coordinates in 
cartesian coordinates. Thus we find 


& — cosóp —sinQó 
¥ = sindp+cosod 
2 = 2 
thus the components are given by 
Ax = X-A=acoso—bsingo 
Ay = y-A=asingo+bcoso 
A, = BAS 


in cylindrical coordinates the components of A are given by 
In spherical coordinates we must also use a transformation matrix. 


sinO 0 cosO 
= —cos0 0 -—sinO 
0 1 0 


D) & =» 
^p > 


where the 3x3 matrix is the transformation matrix that lets you represent the spherical coordinates in 
cylindrical coordinates, this yields 


"b 
| 


sin 0 + cos 02 
— cos 0p — sin 0Z 


= 9 


thus we find the components to be given by 


D D 


A, f- A = asinð + ccos O 
Ag = Î-A = —acos0 —csin0 
Ao Q-A—b 


(c). A=a?+b6+co 
A is not constant with the same argument as above except that 7,0, and all change according to where 
the point is. 


The Divergence of A is given by 








T0, ] 587. 1 dc 
Meses zat 4t cungag up) T BON 
r r 
The curl of A is given by 
1 fa.. ob]. 1f læa ð a 1f dal » 
VxA = ane ine) 4 tu andas 00 OP HZ - Q 


A 
— |tcotOt — 70. 20 
r r r 


in cartesian coordinates the components of A are given by 


» 


sinOcosó cos0cosó —sind 
sinOsinó cosOsinó cosd 
cos 0 —sinO 0 


N> ct 
Il 
D D p 


where the 3x3 matrix is the transformation matrix that lets you represent the spherical coordinates in 
cartesian coordinates, this yields 


sin@ cos df + cos 9 cos QÓ + — sind 


sinOsinóf + cos @sino6 — cos od 
cos 0f -- — sin 6 


» 
| 


Ny t€ 
Il 


thus we find the components to be given by 


> 
= 
» 


: À — asin0cos$ 4- cos0cos $6 4- —csin 0 
A — asinOsino -- bcosOsinó — ccos $0 
:A— acos0— bsinO 


Ay = 
A, = 


N> te 


in cylindrical coordinates the components of A are given by 


0 sinO  cos0 0 f 
ô |= 0 0 1 6 
2 cos0 —sin@ 0 Ô 


where the 3x3 matrix is the transformation matrix that lets you represent the spherical coordinates in 
cylindrical coordinates, this yields 


6 — sin0f--cos00 
6 = 6 
2 = cosOf—sin0d 


thus we find the components to be given by 


Ap = -A-asin0--bcos0 
Ag = 6-A=c 
A, = 2:A-acos0— bsinO 


and in spherical coordinates they are given by 


A =a Ag =b Ap =c 


Problem # 3 
Show that 


[vrav=$ fan 


where A is the area of the closed surface bounding the volume V. Hint: multiply each side by a constant 
vector and use calculus theorems. 


If we multiply both sides by a constant vector field A we get 
A. [ vf aV — A: d f-da 
V A 
if we work on the left-hand side, along with using equation 5 from the front of Griffith's we find 
] A v.nav s I v. (FAjav — [ f(v-A)av 
V V V 


we know that the right hand term goes to zero because the gradient of a constant vector field is 0. 
Invoking the Divergence theorem on the left hand term we get 


] V: UA)dv — d rA da 


we finally find 


A. [ vyav f /^ da 
A. [ vfav 


A- | rav zs A- $ fda 


[vrav= d ta^ 
V A 
Problem # 4 


Show that the velocity (dr/dt = r) is expressed in cylindrical coordinates as pp + p66 + 28 and in 
spherical coordinates as rf + r66 + rsin600. 


ll 
>So 
> 
c 
a 
m 


thus we have 


We know that the velocity can also be written as 


and we know that in cylindrical coordinates 


x = pcosd <(pcosd) = pcos — ópsinó 
d : 
y — psino 4; (P sino) — psinó d ópcosó 
d 
f= 2 ge? 


and we also know that in cylindrical cooridinates 


£ = cosop—singd 
$ =  sindp+cos od 
Z 4 


thus we can write the velocity as 


r = (Ppcosd—dpsind)(cosop — sindd) + (psind+ op cos) (sindp 4- cos 00) +22 
= pcos” 6 + dsin’ 6 + P sin” op + bcos” oo + 22 
= 


in Spherical coordinates we can apply the same technique.and we know that in spherical coordinates 


d : I 
x rsin0cosQ 4; V sin9cosq) = r'sin0cosQ + Orcos O cos ¢ — rsin Osin p 


y = rsinĝsino < (rsin@ sind) = Fsin@sing + OrcosOsing + drsin®cose 


z rcos®@ < (reos 9) — ?cos0 — ÓsinO 


and we also know that in cylindrical cooridinates 


f£ — cosQf sin Qó 
$ = sindp+cosod 
£ Z 


thus we can write the velocity as being (minus some algebra) 
d d d T aa 
2 — i+ 5+ 7 =|7# +706 +rsin0d6 
Problem # 5 
Here we will practice the use of field differential operators and visualize some fields. 


(a). Find the gradient of the scalar potential (x, y, z) = oxy. Provide a clear sketch of the countour 


lines of @ in the € — ? plane and a representation of its gradient field. Such a field is known as a radial 


8 


quadrupole field, and is used in the focusing of charged particle beams and of dipolar (electric or magnetic) 
particles. 


The gradient of $ is given by 


Vo = a EZ + $0» + 5») 


| 
= 
= 
+ 
2 


The sketch is given as 


(A) 


(b). Provide a clear sketch (in the € — ? plane) of the vector field expressed in cylindrical coordinates as 
E — app where a is a constant. Calculate its divergence. If E is an electric field, what charge distribution 
generates it? 


The gradient of E is given by 


The divergence of E is given by 


which is the direction of most change. The charge distribution that produces this charge can be found 
by using Gauss’s law 


V-E=  — 20 => [Gene = 2080 
0 


The sketch is given by 


AAAA RS 
VVVVY 


^e 


(B) 


(c). Consider a large sheet placed on a turntable that is set rotating at angular velocity Q. What vector 
field v describes the local velocity on this sheet? Provide a clear sketch of this field. What is the curl of 
this vector field (now expressed as a vector in a three-dimensional space, i.e including a coordinate for the 
z axis)? If a magnetic vector potential A is defined with the functional form of v (just a different constant 
out front), what magnetic field B =V x A does it represent. 


The vector describing this system is given by 
v =rQÔ 
where r is the distance from the origin and 9 is the angle between any lines originating from the origin. 
The curl is given by 
Vxv= 19 (aye =|202 


The magnetic field is given by 
B=V x v= 202 


The sketch is given by 


10 


(C) 


(d). Now generate new vector field by taking your vector field v obtained in part (c) above, and adding 
(1) a constant field, and (2) a radial quadrupole field from part (a). Provide clear sketches of these two 
fields and compare with that of v. Show explicitly that the curl of these three fields is the same. Note that 
these three fields all represent forms of a magnetic vector potential that describes the same magnetic field. 


We can define a constant vector field as 
k — a& 4- D$ 4- c2 
which can also be represented as 
k = (acos -- bsin09)6 -- (—asinó -- bcosQ)ó -- c£ 
so now we can add these two vectors 
v+k=| (acoso +bsin@o)6 + (—asind + bcos -- rQ)Ó -- c£ 
From part (a) we found 


Vo = a(y& +29) 
= a(ycos +xsind)p + a(xcosd — ysind)d 


adding up these vectors yield 
v+ Vo = a(ycos + xsino)6 + [rQ + a(xcoso—ysing)d 


To find the curl we use 
Vx (vt+k) =Vxv+Vxk 


the right hand term goes to 0 and the left hand term is 


Vxv—202 


11 


also 


V x (v+Vo) =Vxv+Vx Vo 
the right hand term is also 0, so we get 
V xv =202 
we have just showed explicitly that all the curls are the same,. For the firs sketch we can define a 
constant vector field as 
k— u$ 


and then add these two vector fields, as show in the diagram part (a), the second diagram (b) is the sum 
of our original vector field added to the quadrupole field given by (a). 














——| A < 
<] A AAA = 
=< A N 
, —— MAA L 
P A oA IMA P 
L—— Á> 
(a) 
<— ee + | 
<4 = 4 
Aut ONE. A 
UY E 
D'E PENNA "E 
> =_~< P 
(b) 


Problem # 6 
(a) Show that 
[re x A)-da = [ia x (Vf)] da f fA. 
s 5 P 
Using equation 7 from the front cover of Griffiths’s we find 


[IO xA) da= fV x (SA) da-- [IA x (Vf)] da 


invoking Stoke’s theorem on the left hand term yield 


[vx (fA)] da f. [Aa 


thus we find 





[fv xa)-aa= [ax (Vf)]-da-- d 


f^. dl 
s p 
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(b) Show that 
] 8: Ads [ Ax (Vx Base dA x B) da 
V V S 


Using equation 6 from thr front cover of Griffiths we find 
[ BA [ v Amas f A (xmas 
Using Gauss's theorem on the left hand term we find 
[v V-(A x B)d t= (A xB) -da 


Thus we have 
] B: ac | Ax(YxB)dt+ $(AxB)-da 
V V S 


Problem #7 

(a) Let Fy =x?2 and Fz = x + yŷ +z2. Calculate the divergence and curl of Fı and F2. Which one can 
be written as a gradient of a scalar? Find a scalr potential that does the job. Which one can be written as a 
curl of a vector? Find a suitable vector potential. 


The divergence of F; is given by 


and the curl is given by 


VxF = 


o YHo x 
o How 


The divergence of F» is given by 


9 9 ð 
ME Hcet eee 


and the curl is given by 


VxE-— =0 


x Ho x 
< Yve 
N Hone 


We know that F can be written as the curl of a vector potential because it is divergen-less or “solenoidal’, 
to find the vector potential we do 
F; =VxV 


We need to find the vector potential V, so we can find it using 


R  2ĉ 
= |a M LE TENES 
F(-VxV- ox Oy X XE CIN E 
Ay Ay A; 


thus A, is given by 


and the vector potential is given as 


x) 
V= ak+ {y+ be 


where a and b are constants. We know that F2 can be written as the gradient of a scalar potential 
because it is *curl-less". we find the scalar potential by using 


F = -VV 
x2 y? z2 
IDEE qe qut. 
I$ 


thus we know that the scalar potential is 


2 2 2 
Ea e 
EE 


(b) Show that F3 — yz& -- zx$ 4- xyZ can be written both as the gradient of a scalar and as the curl of a 
vector. Find scalar and vector potentials for this function. 


First we must show that F5 is is both divergen-less and curl-less, to show divergen-less and find the 
vector potential we use 


0 0 0 
V-F3= 5, (92) T3, 09 +5 Gy) =0 


thus we must use 





eae ae 
9A, OA ðA; OA ðA, OA 
erli jI GD GBG 
pe ws dy az ox az Ox dy 


Solving for Ax „Ay, and A; we find 





" 0A, OA,N. zy? 
R = [—-——-J£ = 
2 ( y oz s 
ðA; OA yz 
9 — — —— V A — 
oa ( ax aj * 2 
ðA, 9A : 
ni = (2-2: A= 
x y 2 
Thus the vector potential is 
2 2 2 
Yin Ye, Na 
Wee eq. 79 eee 
Qt aU 


For the scalar potential we do 


VxF3= 


Age» 
Y Hore 
3 Ro» 
I 
© 


by inspection. The scalar potential is given by 
F3 = yzk + zx + xy = —VV 
and we find 
V = —xyz 
thus the vector can be written as 
F3 =| —V (xyz) 
Problem # 8 
Check the divergence theorem for the function 
v=? sin@f + 4r’ cos 06 +r’ tan 00 


using the volume of the “ice-cream cone” shown on page 56 Fig. 1.52 (the top surface is spherical, 
with radius R and centered at the origin). 


The divergence theorem states 


[vevav = $ v-aa 
V A 


the left hand side can be solved by taking a volume integral of v. 





l ð, 3z, D 053 
— 4 Les 
7sin8.38 T secos EE ODE ane) 


[ (cos 0)? — (sin @)7] 


1 
V.v — —(rsin0) 4 
r r 
4 
= 4rsin0 + —— 
sin 0 


1 2 
— Arsin0--4r ss = 2(sin0)] 2,4, 16089) 
sin 


sinO 





and the integral is given by 


2x R 1/6 
i V.vdV — f dà [ 4dr [ (cos@)2d0 
V 0 0 0 
which yields 





Now for the area integral must look at the the cone and the surface separately 


a x 20 
f vda = i ] R^dà f *(sine)?a0-+ + | , ^? f " sin(1/6) cos(n/6)4r3dr 
n V3 


12 8 


2nR* 








+7R4 ($) 
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Problem # 9 
A little more practice on vector analysis: consider the vector function F(r) = 6 where @ is defined as 
conventional for spherical coordinates. 
f F.dl 
d 


where C is a circle of radius p in the x — y plane, centered at the origin, and where the integral is 
evaluated with dl oriented counter clockwise. 


(a). Calculate the line integral 


from the diagram we can see that 
dl = pdoo 
given that F(r) = à we find 


fr-a= ["6-(pag6) = [pao = [270] 


(b). Calculate the surface integral 


[0 x8)-da 


where the surface H is a hemisphere that is above and bounded by the curve C used for part (a). The 
surface integral is calculated with da oriented outward. 


from the diagram we can see that 
da = p* sin@d0dot 


and to find what the curl of F we do 


1 0 x 
VSN B [ag sino) £- EJL s t9 


in our case r — p thus we have 





VxF= s [coros — 


thus the integral is given by 


1 A n/2 2x 
VxF)-da- f - tor — -p*sinédedot =p [ odo | ao =[2x 
d x F)-da nro f — 9] - p^sin or =p i cos : $ 


(c). Calculate the surface integral 


] Vx 9) da 


where the surface D is now the disk in the x — y plane that is bounded by the curve C used for part (a). 
The surface integral is calculated with da oriented upward. 


from the diagram we can see that 


da — pdrdbz = pdrdo(cos6#— sind) = 
and to find what the curl is we do 

1 d lod ^  cotO 14 
V F = — i t—|—— = r—— 
m rsin@ EI r - s ? r m 5° 
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in our case r — p thus we have 
1 h 
VxF= pee 


thus the integral is given by 


UE = f, zloto- 6] -pdrdó(cosOf — sin0Q) 
D 


-ff peewee J dus 
sin ð 


= T af dr =| 2mp] 


(d). Verify that Stoke’s theorem holds for both surfaces H and D. 


Stokes theorem states 
jra- ] (Vx E): da — 2np 
S 


thus Stoke's theorem holds. 


Problem # 10 
Let H(r) = x*y& + y?z$ -- z2x2. Find an irrotational function F(r) and a sonisoidal function G(r) such 


thatH = F+G. 
We know that an irrotational function can be described as 


VxF=0 
and a sonisoidal function is described as 


V.G-0 


thus we know that 
VxH=VxF+VxG=VxG 


and we know that 








& $ d 
VxH- 2 » 2 = -y'—2°9 — 172 
xy yz zx 
thus we can find G by 
& $ 2 
0G, x dG, OG, 
VxG= 2 2 2 - (F -P)e- ($- A) Gt ) 
G, G, G, y z x z x y 


looking at the following terms we find 


0G, 2 9G, Er 2 9G, 2 


dy ” z | ox 
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thus 


3 3 3 
y Z x 
G,——— G,——— G,——— 
7 g o9 3 
thus we find G to be 
ü 2. x. y». 
ee eode an 


since we know what H and G are we can easily find F 
3 3 3 
Pee HA 


and it is seen that 


VxF = 0 
V.G = 0 
and by construction we know that 
H=F+G 
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Chapter 2 


Electrostatics 


2.1 The Electric Field 


2.1.1 Coulomb’s Law 


The fundamental problem electromagnetic theory hopes to solve is this: We have some electric charges 
41,92;93;--- (call them source charges), what force do they exert on another charge, Q (call it the test 
charge). The force on the test charge Q due to a single point charge q which is at rest a distance r’ away 
is given by Coulomb’s Law 


| qQ 3i 
— ———r 2.1 
Arey r? Ind 
The cosntant € is called the permitivity of free space. In SI units, where force is in Newtons (N), distance 


in meters (m), and charge in Coulombs (C) 





C2 
=8.85x 10° 
£0 x NS 


In words, the force is proportional to the product of the charges and inversely proportional to the square 
of the seperation distance. As always 


where r" is the seperation vector from the location of q to Q, r is the magnitude and r’ is the direction. 
The force points along the line from q to Q; it is repulsive if q and Q have the same sign, and attractive if 
their signs are different. Coulomb’s law and the principle of superposition constitute the physical input for 
electrostatics- the rest, except for some special properties of matter, is mathematically elaborationof these 
fundemental rules. 


2.1.2 The Electric Field 


If we have several point charges q1,q2,...-,gn at distances r1,/2,...,% from Q, the total force on Q is 
evidently 


1 
F =F, +F)+...=— (29s + 2a) 
4n£0 V rj r5 
or 
F— QE 
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where 


(2.2) 





E is called the electric field of the source charges. 


2.1.3 Continuous Charge Distributions 


Our definition of the electric field Equation 2.2 assumes that the source of the field is a collection of 
discrete point charges q;. If, instead, the charge is distributed continuously over some region, the sum 


becomes an integral 
1 1 
E(r) = —— | sfd 
(r) ATEQ ; ms 4 


If the charge is spread out along a line, with charge-per-unit length A, then dq = Adl' (where dl' is an 
element of length along the line); if the charge is smeared out over a surface, with charge-per-unit are 
©, the dq = oda’, (where da’ is an element of area on the surface) and if the charge fills a volume, with 
charge-per-unit volume p then dq — pdt’ (where dt’ is an element of volume) 


dq — Adl' — oda! ^ pd 


thus the electric field of a line charge is 





Rae / Mr) cap 


B ATEQ r2 


for a surface charge 





E(r) = — / S) egg 


9 4T£0 r 


and for a volume charge 
V TQ 
E(r) = —— | ——fdt 2.3 
(r) ah a (23) 


Equation 2.3 itself is often referred to as Coulomb’s law, because it is such a short step from the original, 
and because a volume charge is in a sense the most general and realistic case. 


2.2 Divergence and Curl of Electrostatic Fields 


2.2.1 Field lines, Flux, and Gauss’s Law 
The flux of E through a surface S is given as 
P= [ E-da 
S 
it is a measure of the "number of field lines" passing through S. This suggests that the flux through any 


closed surface is a measure of the total charge inside. This is in essence Gauss's Law. In the case of a 
point charge q at the origin, the flux E through a sphere of radius r is 


1 q 2.5 a 4 
f£ [ (48) -(Psinededge = z 
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For any closed surface, then 
1 
fE :da — —Qenc 
£0 


where Qenc is the total charge enclosed within the surface. As it stands, Gauss's law is an integral 
equation, but we can readily turn it into a differential one, by applying the divergence theorem 


f 5-4 [iv es 


rewriting Qenc in terms of the charge density p, we have 


Qenc = | par 
V 


[e ma- [(8) dt 


and since this holds for any volume the integrand must be equal 


so Gauss’s law becomes 


Vana” (2.4) 
£0 


Equation 2.4 is Gauss's law in differential form. 


2.2.2 The Curlof E 


Pll calculate the curl of E, as I did the divergence by studying the simplest possible configuration, a point 
charge at the origin. In this case 


1l q 
p- | 1e 
ATE, pt 


what if we calculate the line integral of this field from some point a to some other point b 


b 
fea 


in spherical coordinates, dl = drf + rd06+rsin Odo, so 


1 1 
Arey r2 


eal 


where r, is the distance from the origin to point a and r; is the distance to point b. The integral around a 
closed path is evidently zero (for then rg = rp) 


fE-dl=0 


m 


Equation 2.5 holds for any static charge distribution. 


E-dl= 


therefore 


and hence applying Stoke’s theorem 


2] 


2.3 Electric Potential 


The electric field E is not just any old vector; it is a special kind of vector function, one whose curl is 
always zero. There is a theorem that asserts that any vector whose curl is zero is equal to the gradient of 
some scalar. Because the line integral is indepenent of path, we can define a function 


V(r)z -f E-dl (2.6) 
O 
this is called the electric potential. This allows us to write the electric field as 
E= -VV (2.7) 


Equation 2.7 is the differential version of Equation 2.6, it says that the electric field is the gradient of a 
scalar potential. If you know V, you can easily get E by just taking the gradient of V. 


2.3.1 Poisson’s Equation and Laplace’s Equation 


We found that the electric field can be written as the gradient of a scalar potential. The question arises: 
What do the fundamental equations for E, 


V-E=2 VxE=0 
£0 
look like, in terms of V? Well. V. E — V-(VV) = —V°V, so apart from the persisting minus sign, the 
divergence of E is the Laplacian of V. Gauss’s law then says 
p 


VV == 2.8 
a (2.8) 


this is known as Poisson’s equation. In regions where there is no charge, so that p = 0, Poisson’s equation 
reduces to Laplace’s equation 
VV =0 


2.3.2 The Potential of a Localized Charge Distribution 


I defined V in terms of E. Ordinarily, though it is E that we’re looking for. The idea is that it might be 
easier to get V first, and the calculate E by taking the gradient. In general, the potential of a point charge 
qis 
l q 
V(r) 2 —- 
(r) ATE r 


1 1 
=—— | Z4 
y) s 4 


Oe a l B 


this is the equation we are looking for, telling us how to compute V when we know p; it is, if you like, the 
"solution" to Poisson's equation, for a localized charge distribution. 


or for continuous distribution 


in particular, for a volume charge 
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2.3.3 Summary; Electrostatic Boundary Conditions 


The electric field always undergoes a discontinuity when you cross a surface charge ©. In fact, it is a 
simple matter to find the amount by which E changes at such a boundary. Gauss’s law states that 


1 1 
E:da=—Qene = —oA 
$ a c; dene a 


where A is the area of a pillbox used on the surface. Now the sides of the pillbox contribute nothing to the 
flux, in the limit as the thickness goes to zero, so we are left with 


o 


Jd: el? E 
E — Ebelow ^ B 


above 
Conclusion: The normal component of E is discontinuous by an amount ©/£ọ at any boundary. In par- 
ticular, where there is no surface charge, E + is contnuous, as for instance at the surface of a uniformly 
charged solid sphere. The tangential component of E, by contrast, is always contnuous. 


I _ pil 
E; bove a Ed elow 


the boundary conditions on E can be combined into a single formula 


oO. 
Eabove ~ Ebelow = ear 

0 
where fi is a unit vector perpendicular to the surface, pointing from “below” to “above”. The potential, 


meanwhile, is continuous across any boundary, since 


b 
Vabove ~ Mbelow ^ — f E-dl 


as the path length shrinks to zero, so too does the integral 


Vabove = Vbelow 


However, the gradient of V inherits the discontinuity in E; since E = —VV, implies 


On, 
VVabove — VVbelow — S 
or more conveniently 
WVabove E 9Vpelow ee 
on on Eo 


where 


oV 
^ et VE oi 
on " 


denotes the normal derivative of V (that is, the rate of change in the direction perpendicular to the 
surface). 
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2.4 Work and Energy in Electrostatics 


2.4.1 The Work Done to Move a Charge 


Suppose you have a stationary configuration of source charges, and you want to move a test charge Q from 
point a to point b. How much work will you have to do? The work is defined as 


w= f F-a -o [ &-a- ovt) -v(uj 


notice that the answer is independent of the path you take from a to b; in mechanics, then, we would call 
the electric force “conservative”. Dividing through by Q, we have 


in words, the potential difference between two points is equal to the work-per-unit charge required to carry 
the particle between the points. In particular, if you want to bring the charge Q in from far away and stick 
it at a point r, the work you must do is 


W = QiV(r) —V(~)] 
so if you have set the reference point out at infinity 


W =OV(r) 


2.4. The Energy of a Point Charge Distribution 


How much work would it take to assemble an entire collection of point charges? Imagine bringing in the 
charges, one by one, from far away. The work needed to assemble a configuration of point charges is given 
by 


(2.9) 





this is also the amount of work you’d get back out if you dismantled the system. In the mean time it 
represents the energy stored in the configuration. 


2.4.3 The Energy of a Continuous Charge Distribution 


For a volume charge density p, Equation 2.9 becomes 


1 


(The corresponding integrals for line and surface charges would be f AVdl and f oV da, respectively. This 


can also be written as 
€ 
=e E?^d« 
2 Jspace 


The energy can be thought of as being stored in the field or stored in the charge. We can also say 


€ ; 
3E 2- energy per unit volume 
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Because electrostatic energy is quadratic in the field, it does not obey a superposition principle. The energy 
of a compound system is not the sum of the energies of its parts considered seperately, there are also cross 
terms 


Wiot = Wi + W2 +e f Ei *Endt 


for example, if you double the charge everywhere, you quadruple the total energy. 


2.5 Conductors 


2.5.1 Basic Properties 


In an insulator, such as glass or rubber, each electron is attached to a particular atom. In a metallic 
conductor, by contrats, one or more electrons per atom are free to roam about at will through the material. 
A perfect conductor would be a material containing an unlimited supply of completely free charges. There 
are some general rules when dealing with conductors 


e E=0 inside a conductor 

e p —O inside a conductor 

e any net charge resides on the surface 
e A conductor is an equipotential 


e E is perpendicular to the surface, just outside a conductor 


2.5.2 Surface Charge and the Force on a Conductor 


Because the field inside a conductor is zero, boundary conditions requires thatb the field immediately 
outside is 
o 
E = —û 
£0 


consistent with our earlier conclusion that the field is normal to the surface. In terms of potential 


un ay 


these equations enable you to calculate the surface charge on a conductor, if you can determine E or V. In 
the presence of an electric field, a surface charge will, naturally experience a force; the force per unit area, 
f, is cE. But there's a problem here, for the electric field is discontinuous at a surface charge, so which 
value are we supposed to use? The answer is that we should use the average 


1 
f = oEaverage = 59 (Eabove + Epelow) 


the force per unit area is also given as 
E 2, 
f— — o^ 
2€0 
this amounts to an outward electrostatic pressure on the surface, tending to draw the conductor into the 
field, regardless of the sign of o. Expressing the pressure in terms of the field just outside the surface 
£0 2 
P=—E 
2 
whic is identical to the energy-per-unit volume stored in the fields. 
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2.5.3 Capacitors 


Suppose we have two conductors, and we put charge +Q on one and —Q the other. Since V is constant 
over a conductor, we can speak unambiguously of the potential difference between them 


V=Vi- =- [Ea 


We don’t know how the charge distributes itself over the two conductors, and calculating the field would 
be a mess, if their shapes are complicated, but this much we know: E is proprtional to Q. Since E is pro- 
portional to Q, so also is V. The constant of proportionality is called the capacitance of the arrangement: 


Ce 
V 


Capacitance is a purely geometrical quantity, determined by the sizes, shapes, and seperation of the two 
conductors. 


2.6 Problems with Solutions 


Problem # 1 
Electrostatics with cylindrical symmetry 


(a). Find the electric potential V(z) a distance z above the center of a circular loop of radius R which 
carries linear charge density A. Obtain an expression correct for all values of z (not just z > 0). 


e Total charge on ring: Q 





e Charge per unit length: A — Q/2xa 


e Charge on arc: dq 


We know that 


1 d 
v=- F 
4TE0 r 
since we know that dq is given by 


dq—AXds— ARdó r— V R22 


thus the above expression for the potential is 


2x ARdO ra R 


-maul VETSCUGRJETS 


(b). Find the electric potential a distance z above the center of a circular disk of radius R, which carries 
a uniform surface charge o. 
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e Area of ring: 27ada 
e Charge on ring: dq = o (27ada) 
e Charge on disk: Q = a(t R?) 





We know that 
dq — oda — o2nrdr 


putting this into the expression for the potential we find 


1 i dq V 1 R mor d 
=> - — = — —— — ar 
4TE0 r 4n£oJO0 yz +r 


using a u substitution we find the solution to the integral to which is given by 


u=z +r du=?2rdr (r20:u—5 Z,r2R:u— ZR) 


2 2 
Oo z^--R 6 
v=-2 | udu =|-— [V2 +R? -;| 
4£9 J2 2€0 


(b). Using the potential and the symmetry of the problem, determine the electric field a distance z 
above the center of the circular disk for part b. Check that your answer makes sense: What do you get 
for |z| >> R and for z >> R: Show also that he magnitude of the field at z is proprtional to the solid angle 
subtended by the disk as seen from the point z. 


which gives us 


We know that the electric field can be defined as 


For R < z we can use a Taylor expansion 


(1 2- x)" zz 1 4- nx 





so we find 
o z o quce o R? 
Es——|———1 ze— +s —1|/ 2% —--———2 
2€9 | Vz2 + R2 2€ z 4£0 Z 


and we know that 
q= oOnR? z=R 


thus at R < z we should see a point charge i.e 


E=-— 3? 
Are R2 
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For R > z we can use a Taylor expansion 
(1+x)" zz 1 4- nx 


so we find 


xccl 3 
Z Z 6 Gd Z 
Ecc epe Eats cT 
a ex) E 29 m RU | 





6 z 
E--——|———-1|2- 
2€0 lcm | 
this is only for R > z. 


To show also that the magnitude of the field at z is proprtional to the solid angle subtended by the disk 
as seen from the point z. We will say 


20 0 
o [ ao | sin0d0 = —2ncos® 
0 0 
but the limits of 0 are given by 


0 and 9 — cos! a) 
(Jame 


which gives us 


Z 
o= | ae 

and 

E = I — constant 

Q ATEQ 
thus 

E = kQ 

Problem # 2 


Two spheres, each of radius R carrying uniform charge densitys +p and —p, respectively, are placed 
so that they partialy overlap. Call the vector from the positive center to the negative center d. Show that 
the field in the region of overlap is constant, and find its value. 


We can solve this problem by finding the electric fields for 1 sphere and adding the spheres using 
superposition. It is well known that the electric field of a uniform sphere is given by 








OQ 4 
fE-da x45 a ne 
£0 3 
since the electric field is constant we find 
4 1 p(r)r 
EA — -p(r)n?— E-— f 
3 Pr) f £0 389 5 


so we find that the two spheres are given by 








Les Pre. g | Pr. 
3£9 3£9 
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and so the total electric field is just the sum of these two (superposition), we also know that the vectors 
are pointing towards each other and that the vertical components of the vectors cancel by symmetry. We 
can find the total electric field as 


' x | 3£0 3€0 3€0 i r 


but since we know that the total distance between this two vectors will always be the same we can 


write the r_ vector as —r_which would give you just give you d (which is constant) and thus your total 
electric field would be constant. i.e 


d =r} +r- =r} r- 


Problem # 3 

Consider the so-called “screened Coulomb potential” of a point charge q that arises, for example, in 
plasma physics 
—r/À 





q e 
V(r) = —— 
e 4n£90 r 


where À is a constant (called the screening length) 
(a). Determine the charge distribution p(r) that produces this potential. Sketch this function in a 


manner that clearly describes all of its characteristics (i.e what's the best way of representing this three- 
dimensional charge distribution? Use it and explain what your potting. 


We can solve this problem two ways, you can take the Laplacian of a Laplacian or just use the following 
identities 


VENA EVEEN) 
and we also know the following relationships 
E-—-VV V.E-V.(-VV)- e 
0 


thus we know that 
P(r) =€0(V-E) = €0(V-(—VV)) 


so we need to find what the divergence of the electric field is, by finding the gradient of the potential 
first and taking the divergence of this 


Q q e^ q Q e^ e ^ r 
V = — —— = — —|f 
E or E r 47£0 Or r r? [i $ A 








thus we find A 
et rl. 
letting 
E^ r f 
foe IE and 8—3 
we find 


f f r 
V.E— e^ (1 jV LV pu -) 
7 T go $ ltz 
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but we know from Griffith page 50 that 


so we find 





The plot is given by 


where the plot is a countour map that describes the distribution with respect to the radius. 


(b). Show by explicit calculation over p(r) that the net charge represented by this distribution is zero. 
Explain how this observation is easily derived using the integral form of Gauss’ law. 


the total charge Q can be found by 


—rfA 
= n er deu l Te 
Ones I p(r)dt an B 1 i (1+ x) And? (r)dt Re n? is 


which gives us 
re ^ 
[ear 4n — af 7 aa ot 


where the term on the right hand side is just the Gamma function and the delta function integral is just 
4n by definition. So we find that the total charge is zero 


oo —rf[A 
re 
i re me 
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Q(tot) — J p(r)dx — 2L [An — 4n] - 0 
space 


this tedious calculation could have been solved by using the divergence theorem and Gauss's law. i.e 


n» E)dV = fE-da= Ge 
£0 


[O Dav = f (av -o- fan - e 


but we know t 





Since we now know that we have a +g charge due to the delta function and a -q charge from the integral 
of the potential contained within a surface S then the surface has the same number of field lines coming 
in that are going out and thus the surface integral goes to zero because the divergence goes to zero thus 
implying that dene is Zero. 


Problem # 4 
(a). Show that the electrostatic energy of a uniformly charged solid sphere, with total charge Q and 
radius R is 3Q7/(20meoR). 


It is assumed that we know what the electric field is in the inside of the sphere and on the outside of 


the sphere is. i.e 
Hx d 4 


Ej, — — — rf = ——f 
" — 4n£oR? ou Amegr2 
and using the following expression for the electrostatic energy 
€ 
wWw=2/ Ed 
2 space 


thus we must consider two integrals, the contribution from the inside and the contribution from the 
outside. Knowing the limits of integration we find 


2 

q i d 
w= 2 a 

Se nn 0 ire [a arj- S60 = lata] 


3q? 
20nxego9R 





W = 


(b). Use the result above to compute the electrostatic energy of an atomic nucleus with Z protons and a 
total of A nucleons, using an approximation for the nuclear radius as R = (1.2 x 107 !5m)A" 3. Give your 
results in units of MeV times Z2/A!/? 


this is just a plug and chug problem, figure out your constants and plug away, i.e 


3(Zq)* V 
mue c. D cec [D20MENE 
20n€o(1.2 x 10-5mJAU2 AFUE 


(c). Calculate the change of electrostatic energy when a uranium nucleus undergoes fission. 
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We know that U7*> decays as 
U?5 532 ky. L1 Ba +3n 
thus we know that 
AW = Wreac — Wu = |(Wkr + Wg) — Wu 
thus 


2 2 2 
AW =0.720MeV ( 2 2 ) e 


923 * 1411/3] — 2351/3 


which is a large amount of liberated energy. 


| — —347.02MeV 


Problem #5 
Evaluate the following integrals 


(a). 
f (2 49-a+a2)83(r—a)dt 
space 
where a is a fixed vector and a is its magnitude. 


We know that the delta function is zero everywhere except at r = a thus we must evaluate this integral 
at a i.e 


[ (9-9) — fa) 
thus this yields 
/ (E r-a c a2)8 (r — aac — 3a? 
space 


(b). 
i Ir — b|283(Sr)dt 
V 


where V is a cube of side 2, centered on the origin, and b — 4$ 4- 32. 


We know that we must evaluate this function at r — 0, which yields 
J ir — b|283(5r)dt = | — b|? = 
V 
(c). 
n (44 P(r-©) +c4)83(r—e)dt 
V 


where V is a sphere of radius 6 about the origin, and c= 5X+ 39 + 22. 


We know that the delta function is zero everywhere excpet for where it is defined, so this means if the 
limits of the integral are not within the specified value at which the delta function is equal to 1 then the 
integral is zero. 


c= y 52 +32 +22 = V38=6.16 — 6.16 » 6.0 


thus this implies that the magnitude of c is outside of the range and so 


[Atreo 4-c*)8 (v — e)dx — [0] 
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(d). 
I r-(d—r)83(e—r)dt 


where d = (1,2,3),and e = (3,2,1), and V is a sphere of radius 1.5 centered at (2,2,2). 


From our experience with delta function we now know that this will be given by 
i, r.(d — r)8 (e — r)dc — e (d—e) =[—4] 
V 


Problem # 6 
The farad is actually an enormous unit of capacitance. To illustrate this, treat the Earth as a conducting 


sphere and find its capacitance. 


The above figure describes the earth capacitance problem, Where A; is the radius of the earth and also 
R» > R; and also a positive charge Q on the outside surface and a negative charge Q for the earth. 


We know that the capacitance is defined as 


Q 
C= = 
V 
where Q is the charge and V is the potential. The potential is given as 
Ro 1 Ro 1 1 1 
VSS E-dr- -—— | £ (cosn)dr — —— —— — 
Ri 4TE0 R r 4Treo | R2 Ry 


because the electric field points in the —£ direction we must consider the cosz term. Thus we find 


1 1 
es ek 
Amey |Ro Ry 
and so the capacitance is given by 
R4R54n£9 


C-————— R R, R(,—6,378k 
RoR, 2> Kı 1 m 
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so we find 


C — Rj4n£o — | 700 uF 
Problem £ 7 


Consider a system of n conductors that are prepared with total charges of Qj , with ic(1,2, 3,. . . 
n) labeling the conductors. The surface charge densities on each of the conductors are described by the 
functions o; , for which, clearly 
Qi — 7i Oda; 
Si 


with the integral being over the surface S; of conductor i. 

(a). Find an expression for the potential on conductor i. This you may obtain as usual by summing 
over the contributions to the potential from all the surface charges on all the conductors (including the 
i — th conductor). 


Since we know that there i;;j and j; conductors, to find the potential of the i,, conducter from all the 
jin conductors, we must sum over all the j; conducters contribution to the potential, i.e 


1 Qj j ojdaj 
jal? Si 


i 
4T£0 S; fij lij 


where r;; is the distance from the i, conducter to the j,, conductor. 
(b). Let us define a coefficient pi; as 








1 fidaj 
Pij = = 
ij 


4T1£9S j4S; T 
Here, r;; denotes the distance from one point on conductor i to the point on conductor j at which the 
integrand is evaluated (as we integrate over the surface of j ). In this integral, we make us of a relation 


5; - (cf; - iy; 
J 


in which S$; is now the area of the j,; conductor. The dimensionless function fj describes how the actual 
charge density differs from the average. With this definition, show that the following set of equations hold: 


Vi = piuQit+pi2Q2+...+ pinQn 
Vo = p21Q1 + p22Q2 +... + P2Qn 


Va = PniQ1 + Pn2Q2 bus + Pano, 


The coefficients p;; are known as coefficients of potential. They express how the potential on conductor 
i will vary (linearly) upon changing the charge Q; on conductor j . We note here that the coefficients no 
longer contain any reference to the specific potentials or charges placed on the conductors - they simply 
reflect geometric properties of the conductors and their placements. 


(c). By considering the quantity p;i;Q;O; , show that pj; — pj;. This nice relation expresses the fol- 
lowing fact: if a charge Q on conductor 1 brings that conductor to a potential V , then the same charge Q 
placed on conductor i would bring conductor j to the same potential V. 
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This can be shown by 


iQj jdaj (Qi id dj 
pijQiQ; — oig; Hu e Ea = pjiQjQi 


Ane S; S; Fij 4T£0$; Si Vii 


and using the fact that 


This will shows the symmetry of the problem, i.e whatever conductor i sees in the presence of con- 
ductor j of charge Q will be exactly the same if conductor i of charge Q as seen by conductor j. 


(d). Consider a situation with only two conductors that are used as a capacitor. What is the capacitance 
for this system expressed in terms of the coefficients of potential p11 , p12 and po» ? 


if we know that 


Vi = piuQit+pi2Q2 
V p21Q1 t p22Q» 


we also know that 
p12 = p21 


thus we can find the potential by using 


AV =V - Vi = (p2 - p11)Qı + (pz - pi2)Q2 


and we we assume that Q1 = — Q2 we find 


c= 1 


AV — pucpi-cpo 


Problem # 8 
Suppose that a parallel plate ca-pacitor has rectangular plates but the plates are not exactly parallel. 
The separation at one edge is d — a and d + a at the other edge, where a < d. Show that the capacitance is 


given approximately by 
E0A a? 
Ce 1 
d ( 2 


where A is the area of the plate. A figure is given by 
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Thus we know from the diagram 


A — 4lw 
and knowing from Gauss's law 
f E-da = — 
which give us 
odA o 
EdA = — =>E=—(-2z 
2€0 2€0 (o2) 
to find the potential we must use 
z5 d E-dl 
and knowing that 
dl — —2ds 


where ds is in the direction of the field. Thus we find the potential to be given by 


Oo , B 
v-- | 3) Cà) - = -s 


which S is the the distance between the plane that is a function of x which is given by the equation of 
the line in the form a 
S=d+—x 
w 


thus we now know that the capacitance is given by and denc = OdA 


ee d w dx 
] af uu e 
oe —w —w (d F 2) 


using a u substitution, in which we 2 usi = the solution. We know how to do u substitutions at 


this level Cans n a (1+5)] —]n fa (1-2)]| 


which if we do a bit of algebra with a knowledge of how the natural log functions work along with the 
Taylor expansion for the natural log function which is 


2 3 


x x 
1 Diya Se 
n(x4-1) zx p Us 


we find 
lw|a a a? a qd a 
C —2€ —|-|---~=5-=— 
a $- 2d) 3 | | d 24 all 
the even terms cancel and thus to third order we find, with a bit of algebra 


1[2a @ 1 A£0 í 
n E xm "e og [t+ m- am | "ER 


Problem # 9 


(a) Two grounded conducting planes meet at an angle of 0 — 60? at the origin. A point charge q is 
a distance ro from the origin along their angular bisector. What image charges are needed to satisfy the 
boundary conditions? (Hint: their arrangement is going to be highly symmetric) 
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We can see from the diagram will require 6 total charges and 5 image charges. This is necessary to 
meet the boundary condition implied by the use of images. We can see that there is strong dependence in 
the angle between the plates. You must balance the charges. where the total charge is zero. 


What is the r depedence of the potential along the angular bisector for r >> ro? 
We can find the r dependence is given by 
q |1 1 1 1 1 1 
B Jeep 
4T1£0 rı Y? r3 r4  Y5 re 


if we pick an arbitrary point (0,0,z) we can find a relationship for the distance between these charges 
and that point. The diagram is given by 


092 | , 


r6 


<---------> 
1 





knowing that the distance between each pair of charges is the same because it is a equilatiral triangle 
and by symmetry we find 
r4=rs andalso r= r3 


and also from the diagram we can see 
r=r—Ro re=r+Ro 
Thus we find that the potential is given by 


q 1 2700 1 


= 4mtég9 |r —Ro ro ra r+Ro 
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combining the first term and the last term gives us 


<1 q 2Ro 2:172. 
^. 4T£0 r? +R? m2 14 





From the diagram we can also infer that 
2 2 0 
SE Pe ae uu 
r2 | 0 (r 2 ) 


3 Ro V? 
ri = Fo (r+ 2) 


thus we find that the r dependence on the potential is given by 


q 2Ro 2 2 
is | a reer ree a ee 
0 
0 3R + (r- 3) iR + (r+ 3) 


(b) Repeating the above problem but for different angles 0 between the grounded plates. Show that if 
0 — 2r/n, where n is any even integer, then n — 1 image charges are required to ground both sides of the 
wedge. Show that if, however, n is an odd integer (e.g. n — 3) then an image solution cannot be obtained. 


6-60. 
o o 
a ‘ 
6-180 
9-90 o 
Z^ 


o 
4 -4 


and by symmetry we know that 





We can see that there is a dependence on n by 


360 360 360 
2 sjes S loea a = —1)= 
a T44- T-N 6(6-1)=5 


60 
from these three examples we can see that 


= =n total number of charges 


and 
n — | image charges 

we can also see that O must be even to allow an even number of charges that to be used in order to 
satisfy the boundary conditions. 

Problem # 10 

A conductor is formed with a grounded conducting plane spanning the x — y axis (z = 0). Upon it is 
attached a hemispherical conductor of radius Ro that extends into the half-plane z > 0 and whose center is 
at the origin. A charge q is placed at the point (x,y,z) = (0,0,z) with z > Ro. 
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(a). What is the electrostatic potential V(r,9) with r and 0 being conventionally defined spherical 
coordinates (by symmetry, the potential must be independent of the azimuthal angle ©)? Hint: there will 
be three image charges. Answer: 


z — R? 


z4/ z6 +R 


T] 


A diagram is given by 


V(r,0) 





-q 





Thus we know that 


Ce ee ee 
4TE0 | r1 r2 r3 r4 


Using geometry and symmetry, along with the the example in Griffiths’s we find that the relationships 
for the all the radius are given as 


rı — rsim 0-4 (rcos0 — z)? 
RI 

rj? = ?sin20 + (reso - 2) 
z 


and by symmetry we find that the other two radius are given by 


z 
rà? — rsin 04 (rcos0 4- z)? 


RA? 
rn? = sino + (reos. 2) 


so we know that the potential is given by after cleaning it up a bit 


q 1 Ro Ro 1 


ee en e e T ML Led 
(n ) 4TE0 J/r2 +72 —2rzcos8 2 2R5rcosO Rj 2 2R5rcos0 Rj vr? - z2 4-2rzcos0 
ZAT T ZAJ T T m 


Z 


and we know that the total charge distribution induced on the plane is given by 


and so we find that the derivative of the potential with respect to 0 is given as 


dV rR? sin O rRQ sin O rzsin® 
do ^4 Bf 2 ax eae mer 3/2 
* (7? - a i =) x (P 4 aes js 3) (r? - z2 — 2rzcos0) 





rzsinO 
q 
(r? +z? — 2rzcos oy? 


thus we know that 


pies eh Rosin E Risin zsin 
xd e G __ 2RjrcosO gyn 1C i zu a aye (24.22 — 2rzcos 0)! 
Z 4 Z 
q zsinO 
Gs) (r? +z? — 2rzcos — 





and since we are integrating over the plane this means @ = 5 and da = rsin@dodr and 


dp= oda 
space 
which will give us 
S T ef — — — — 
Ip = 
An Ro é HONOR 2 (248) (r? + 22)3/2 (r2 + 22)3/2 
which becomes 
mT o Rdr c 
ARV On 2 E (2x22 


3 zt 
0 2 (e+? +4) lg jaa 
solving this integral using u substitution yields 


_ Ro(-2) 


-1 
gu “ED TP X 
T ^ RA 


simpyfying this term yields 


z R z z? +R? 
Qp — —d | —7———— ae Ot 
JR c2 5 J2cRe z2-R, 
and knowing that 
qh q 
dptdp- -q ———-1-— 
—q —q 
thus 
di. 
—q 





Chapter 3 


Special Techniques 


3.1 Laplace’s Equation 


The primary task of electrostatics is to find the electric field of a given stationary charge distribution. In 
principle, this purpose is accomplished by Coulomb’s law, in the form of Equation 2.3 
1 f 
E(r) - — | —p(r)ac 
Unfortunately, integrals of this type can be difficult to calculate for any but the simplest charge configura- 
tions. We can attempt to solve for the potential using 


1 D oue 
V(r) = us | oe (3.1) 


Still, even this integral is often too tough to handle analytically. In such cases it is fruitful to recast the 
problem in differential, using Poisson’s equation 


1 
V?V = —p 
£0 
which, together with appropriate boundary conditions, is equivalent to Equation 3.1. Very often in fact, we 
are interested in finding the potential in a region where p — 0. There may be plenty of charge elsewhere, 
but we're confining our attention to places where there is no charge. In this case Poisson's equation reduces 
to Laplace's equation 
V?y «0 
or, written in cartesian coordinates 
3V à 3V 2 PV — 
o? ay? oo 
this formula is so fundamental to the subject that one might almost say electrostatics is the study of 
Laplace's equation. At the same time, it is a ubiquitous equation, appearing in such diverse branches of 
physics as gravitation and magnetism, the theory of heat, and the study of soap bubbles. 


3.1.1 Boundary Conditions and Uniqueness Theorems 


Laplace’s equation does not by itself determine V; in addition, a suitable set of boundary conditions must 
be supplied. The proof that a proposedset of boundary conditions will suffice is usually presented in the 


4] 


form of a uniqueness theorem. There are many such theorems for electrostatics, all sharing the same 
basic format-TI'll show you the two most useful ones 

First uniqueness theorem: The solution to Laplace's equation in some volume V is uniquely deter- 
mined if V is specified on the boundary surface S. 


Proof: Imagine a region and its boundary, there could be "islands" inside, so long as V is given on all 
their surfaces; also the outer boundary could be at infinity, where V is ordinarily taken to be zero. Suppose 
there were two solutions to Laplace's equation 


vVyw-20 Vwz20 


both of which assume the specified value on the surface. I want to prove that they must be equal. The trick 
is to look at the difference 
V3 =Vi —V2 


This obey’s Laplace’s equation 
V?y, « V?y, - V?v, 


and it takes the value zero on all the boundaries (since V; and V2 are equal there). But Laplace’s equation 
allows no maxima or minima-all extrema occur on the boundaries. So the maximum and minimum of V3 
are both zero. Therefore V3 must be zero, and hence 


Vi =V2 


The uniqueness theorem is a license to your imagination. It doesnt matter how you come to your solution; 
if a) it satisfies Laplace’s equation and (b) it has the correct values on the boundaries, then it’s right. 


3.1.2 Conductors and the Second Uniqueness Theorem 


The simplest way to set boundaries for an electrostatic problem is to specify the value of V on all the 
surfaces surrounding the region of interest. And this situation often occurs in practice. In the laboratory, 
we have conductors connected to batteries, which maintain a given potebtial, or groundm which is the 
experimentalist’s word for V = 0. 

Second uniqueness theorem: In a volume V surrounded by conductors and containing a specified 
charge density p, the electric field is uniquely determined if the total charge on each conductor is given. 
(The region as a whole can be bounded by another conductor, or else unbounded.) 


3.2 Multipole Expansion 


3.2.1 Approximate Potentials at Large Distances 


If you are very far away from a localized charge distribution, it “looks” like a point charge, and the potential 
is, to good approximation, (1/47€9)Q/r, where Q is the total charge. We can expand the potential as 
1 C 1 INn / / / 
V(r) = qus 2 cnm fo P,(cos0')p(r’)dt (3.2) 


n=0 r( 


or more explicitely 
V(r) = poe T f oae 7 [ r'eostfotr at S [y P aisigz. t p(r)d« 4- 
4n£o |r r? r? 2 2 p 
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this is the desired result- the multipole expansion of V in powers of 1/r. The first term n = 0 is the 
monopole contribution. The second (n = 1) is the dipole term; the third is the quadrupole term; the 
fourth is the octopole term and so on. As it stands Equation 3.2 is exact, but it is useful primarily as an 
approximation scheme: The lowest nonzero term in the expansion provides the approximate potential at 
large r, and succesive terms tell us how to improve the approximation if greater precision is required. 


3.2.2 The Monopole and Dipole Terms 


3.3 Problems and Solutions 


Problem # 1 

Consider a cube made of 6 conducting plates of size a x a , that encloses the space —a/2 < x,y,z <a/2. 
The plates at z= + a/2 are held at potentials +Vo , respectively, while the other four sides are all grounded 
(V 2 0). 


(a). Solve the boundary value problem for V (x, y, z) inside the cube. 
First we must solve the Laplace equation 


eV gv gv 


ox? = oy? a dz? 
As always, we look for solutions that are products 
V (x,y,z) =X (x) ¥(y)Z(z) 
Plugging this into the Laplace equation we find 
Lax... tay. Vez 
Sot Sa Sa 0 
X dx Y dyt Z dz 
thus we know that 
1 X | dy Id 
xq cU yon ga 
which implies that 
C1 - c2 c3 —0 


we can now solve this second order differential equation. One of these coefficients must be negative, 
this can only be c3 because we know that this must be an oscillatory solution. In this case we find 


c= -(cyte2) =P eask c3=-(P +h) 
and so : I 
X 


du 
d dy 2 


— =-(P+R)Z 
dz Np 


43 


the solutions to these differential equations are 


X 
ce 
| 


— Asin(Ix) 4- Bcos(Ix) 
Y(y) — Csin(ky) + Dsin(ky) 
Z(z) = Bev p pe V Hs 


Boundary conditions 


our boundary conditions are given as 


(i) V — 0 when y 2 £a/2 
(ii) V = O when x 2 £a/2 
(iii) Vo = +Vo when z= +a/2 
from (i) we can see that 
Y (a/2) = Csin e -- Dcos (i5) —0 Y(—a/2) 2 —Csin (i5) a dyes (5) 20 


and from (ii) we can see that 


X(a/2) — Asin ($) + B cos ($) =0 X(—a/2) = —Asin (e) +Bcos (KS) =0 


from the first and second boundary condition we find 








C=0 A=0 
and also 
pk ae pane game 
a a a? a? 


where n and m are positive integers. Since the given boundary conditions imply an odd potential we 
can see that 


Z(-1) » -z()  |E VP *P* e ge Ps] — — [ge PP pu VP 


and so we find that 
F=-E 


and this now gives 
Z(z) =E le Heu d = E2 sinh( V12 +k2z) 
combining the remaining constants we are left with 
V (x,y,z) = Ccos (=x) cos (5) sinh( v12 + k?z) 
where C has absorbed the other constants, we can also write this as 
V (x,y,z) = Ccos (=x) cos (5 sinh (Sv +z) 
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we can now write this as a double sum over the integers n and m 
V (x, y,z) zs 25 n,m Sinh (= vnm m cos(nnx/a) cos(mny/a) 


and finally the last boundary conditions (iii) gives us 


V(x,y, ta/2) = y y C; msinh (7 vn +m?) cos(ntx/a)cos(mny/a) = +Vo (3.3) 


n=|1m=1 


to determine the constant Cp „n we multiply both sides by cos(n'mx/a) cos(m/my/a) where n' and m' are 
positive integers and integrate, this is known as Fouriers trick. For the left hand side we find 


a/2 / a/2 / 
2, sinh (sv? +m?) Cun BR cos (=x) cos (=) esf cos (55) cos (25) dy 


we also know from Fouriers trick and orthogonality that n — n' m — m' and so we get 


a/2 a/2 
sinh (2v n? +n?) Gf cos? (=x) dx | cos? (=) dy = Cnm sinh (2v n? m) - 
a 


—a/2 a —a/2 4 


2 


and now for the right hand side of equation 1 we find 


a/2 / a/2 / A 2 
zv f cos (=) ax | cos (25) dy = Vo : 5 
—a/2 a —a/2 a nmi 


and finally we can put these two solutions together to get 





; nī . mu 
sin (=) sin (=) 
a a 








2 2 
4 
Cn.msinh é vn? +m?) e us Vo- = 5 sin (=) sin (=) 
' 2 4 mx? a a 
this yields 
16V, 1 
Cim = E 2 —————— — —- for n and m odd 
nM” sinh (8v n? 4- m?) 


if n and m are even then C, ,, is zero. We can now write the potential as 





(b). The potential is obviously zero at the center of the cube. What is the electric field there? 


To find the electric field we use 


z ð ð ƏV 
Ee me ouo 
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and since we are evaluating the electric field at the origin we know that 


oV _ oV zz 
ox Oy | 

this leaves us with 

cosh (z n? Tm) 


sinh (Ev n? 4- m?) 





Š oV 16aV 1 TX T 
Ex yz) = 77t — Y — eos (75 cos (22y) 


v/ n? 4- m?2 


and so the electric field at the origin is given by 


16Vo 1 i 


E(0,0,0) =|+ 


GV m=1,3,5 ™” sinh (8v? Fm?) 





Problem # 2 


(a). Suppose that on the surface of a sphere of radius R there is a surface charge density o(0) — 
60 cos? 0, where 09 is a constant. What is V(r,0) inside and outside the sphere? What is the electric field 
E(r,0) inside and outside the sphere? Confirm that the field is discontinuous at the surface of the sphere 
in accord with the given charge distribution. 


From equation 3.65 in Griffiths we know that 


oo 


V(n8)2Y (ar + x) P;(cos®) (3.4) 


1=0 
this is the general form of a potential in which we assumed azimuthal symmetry in spherical coordi- 
nates, this was derived in Griffiths on page 137-139. Looking at the inside of the sphere when r < R we 
find that B; = 0 for all / otherwise the potential will diverge. Looking at the potential outside the sphere 


we can see that A; = 0 otherwise the potential would not satisfy the boundery condition that V — 0 r — oo. 
From this we can see that 


V(r « R,9)— Y Air P(cos0) 
i=0 


p 
V(r > R.8) — Y 7 Pi(cos®) 
i=0 


we know that the potentials are the same at r — R and continous so we can set these expression equal 
to each other 


co 


A|R'P;(cos®) = ¥ ——P,(cos@ 
L 1R' P;(cos®) L R 1 (cos 8) 


and we can see that 
Beak giu 


since there is a discontinuity in the radial derivative we can use equation 2.36 


OVout Win mE 1 
( a ) "oue 
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from this we find 


co oo fe 1 
p (+1) atts x5 i (cos0) a 'P;(cos®) = E 
substituting the above expression for B; and some algebra yields 


1 
(21 + 1)A)R'|P,(cos®@) = c 0D cos? 0 
0 0 


Ms: 


V 


we can use our knowledge of Legendre polynomials to solve this problem 


12 P 0) \ 2 
cos?0 — (ra(cos0) + 5) ao (Pies 0) 4 LT 


3 


looking at only the / — 0,2 terms in the sum and using this on the above expression yields 
lo 20 
AoR" ! Po(cos0) 4- SASRP»(cos0) — 3 = Po(cos 6) + 3 = Pa (cos 0) 
0 0 


this allows us to solve for Ao and A; as 


OoR 200 |. | 
A LSU. Cangius 
De Beh. Ot ASE 


so the potentials for inside and outside the spheres are given as 





OoR 269r? 
R,0) 2 —— 4 ——P 0 
V(r < R,0) 36; + TSe0R (cos 8) 
ooR? 200R* 
Vin = RO0)— ——P. 0 
OE Rma rese d OP) 
We can find the electric field by 
y oV lov. 
E = -VV = — | —î+ -— 0 
ac Las | 


for inside the sphere we find 


> 4 roo 
Ein = —Win Ne T ( 


15 R£o 


3cos?8—1Y , roo3sin(20) 
2 5REo 


knowing that the derivative of the 0 term is 


oS. do il 3 
= Sl 20] = ~sin2 
703 °° 0 2462. 4- cos 26] 5 sin 0 
which can be simplified as 
2 roo |3 . 2 
Ein = —VVin = = sin(2 -1 
Vi 15 Re EET 0)0 — [3cos^ 0 — 1] | 


for the outside we find 


4 1g | g* A 
s. (3eos'0 — D) f+ 5 vig, (sin20)6 


Eout = —VVout = 





to show that these fields are discontinous at the surface, one would just have to show that at r = R E; x 
Eou and this can be seen from the above expressions. 


(b). Now repeat the above exercise for the surface charge density o(0) — oo sin? 0. Note that in this 
case 











2Ro 2Ro 2 
V(r0) — 3 — 156) (g). Ph(cos0) r<R 
f 2Ro, 2Roo (R 

7 E 158r (4) P;(cos8) r>R 


For this problem we can just start with 


o 
lona JAiRI!P(cos9) — $999 = (1 — cos?0) 
1—0 £0 £0 


and only looking at the / = 0,2 case in the sum we find 
2 1 
AoR" ! P,(cos0) 4- SASRP»(cos0) — E Po(cos0) — 3 Pi (cos 0) — 3 Po (cos 0) 
0 
setting like terms equal to each other gives 


200R 200 
pS A= 
3£0 15R£0 


and so finally we find that the potentials are given by 


2Ro9 E 2Ro9 (> i P3(cos8) 


R —— 
WSIS eic sen 





and also 








2R? 2Roy (R 
V(r2R,0)— Ue a(z ) P»(cos9) 


3€or 1589 


and the electric field is given as 


Ein = —Win = 





and for the outside we have 


2R°09  600R* 1 ooR* x 
UE cL = (cos Bos D)e f+ 5 T sin 200 
0 


Eou F —V Vou =| ITI 
' ] ( Beor? 152074 
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Problem # 3 

Solutions to Laplace’s equation in two dimensions using complex functions: One sometimes encoun- 
ters a boundary-value problem that is essentially two-dimensional, in that, by the nature of the problem, 
we can assume the electrostatic potential is constant along one Cartesian direction. In this case, it may be 
possible to find the solution to our boundary-value problem among the many solutions obtained by anal- 
ysis of analytic functions of a complex variable. For this, let's consider the complex number ¢ = x + iy, 
where x and y are the remaining two Cartesian coordinates for which we wish to determine V (x, y). We let 
F(C) be an analytic function of ¢- meaning that it is a well-behaved function with well defined derivatives. 


(a). Show that F'(C) satisfies Laplace’s equation, i.e. that 


2F PF 
d PE tu 


2p 
VF= ax oy? 


To show this we begin with 
F(E) =F(x+iy) =x+iy 


we also know that 


OF OF oz _ OF 
ox az ox Oz 
PF _ 9 [ar] _ a [ar] _ vr 
ox) ox[|O0z|] oOz|ox] 2 
and now for y we find 


OF c her GP 
oy |— azoy z 


dF . ð [F] Of0F) » F FF 
— LI i . — — SS ——— = * os —À RÀ LL 
dy2 dy | dz] az | dy oz? oz? 
because we know that the differential operator is commutative and we can say that 
00 dda 
oxdy  dyox 


so we have shown that all nice well behaved complex functions satisfies Laplaces equation. 


oF OF _ F "a 


2 = — — 
Ve do ea o2 o2 


=0 

Given this fact, we have indeed found two solutions to Laplace’s equation in two dimensions. Express- 
ing F = g(x,y) +ih(x,y), then the real-valued functions g and h are both such solutions. 

Deriving from a single complex analytic function, the functions g(x,y) and h(x,y), known as conjugate 
functions, possess an important relation to one another. By considering first derivatives of F (G) one can 
demonstrate the validity of the Cauchy-Riemann relations: 


dg Oh og oh 


ox Oy Oy 
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(b). Show that the electric field obtained by setting V (x,y) = g(x,y) is everywhere orthogonal to that 
obtained for V (x,y) = h(x,y). Based on this observation, one can state that the curves h(x, y) = ho describe 
the field lines for the potential V (x,y) = g(x,y), and visa versa. 


To show that these two solutions are orthogonal we must take the dot products of the two solutions 


AC e z zai OBR ES 
E(g) — -VV—-Vgisy)- see xs 
Eh) = -VV- -Vh(sy) - — | orga Og 
ox oy 
the dot product must be zero for these two solutions to be orthogonal. 
4 =,  [0g, 0g.| [0h, Oh.| dgdh dgoh _ 
Ble) - [sse sl- [ut t| "as ms 


but since we know that the Cauchy-Reimann equations are 
dg oh Oh og 
ox oy Ox oy 
using this we find 
dgdh dgoh _ 
oxox ðxəðx — 
(c). One application of this method: Find the potential V (x, y) in the wedge-shaped region between 
two plates (infinite half-planes) that intersect at the z axis at an anglea < 1/2, and are both at the potential 
Vo . You might consider the real or imaginary parts of the complex function F (G) = AC? + B and attempt to 


match the given boundary conditions. Other nice solutions obtained by this method are the field produced 
at the edge of a semi- infinite charged conducting plate, and the fringing field of a parallel plate capacitor. 


Since we know that we can write a complex function as 


F(&) — F(x-- iy) 2 Ar"e"?-.-B — & =re'® =r{cos(o) +isin(6)] 
where p is now the angle defined in the complex plane and r is the distance from the origin, thus we 


find 
F(&) =Ar?|cos(po) +isin(po)] +B = Ar? cos(po) +B +Ar?isin(po) 


since we know that both the real solution and imaginary solution satisfy the Laplace equation we can 
solve for the real part and we are guranteed that this is a unique solution. 


RV(r,0) =Ar? cos(po) +B (3.5) 


Imposing boundary conditions we find 


V(nró — 0)=Ar?+B=Vo 
V(ró-—4«) - Arcos(po) 4 B — Vo 
Subtracting the first equation from the second we find 
2nt 


Ar? — Ar? cos(pa) =0 p= Es 
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now if we add the two solutions we find 
Ar? +Ar? cos(pa) + 2B = 2Vo 
using the expression for p we get 
2Ar?+2B=2Vo B=Vo—Ar? 
inserting this into equation 3 we find 
Ar? cos(po) + Vo —Ar? = Ar?(cos(po) — 1) + Vo 


we can also see that 


x=reos(pp) cos(po) =" r= Vx +? 


V(x,y) 2 |AxG? -y?) 0 DP? — AQ? - y??? E vo 


and so we find 


Problem # 4 

Solve Laplace’s equation by separation of variables in cylindrical coordinates, assuming there is no 
dependence on z (cylindrical symmetry). [ Make sure you find all solutions to the radial equation; in 
particular, your result must accomidate the case of an infinite line charge, for which of course we already 
know the answer. ] 


We can begin by writing the potential as a product 
V(s,0) = S(s) ®() 
and solving the Laplace equation in cylindrical coordinates 
1a / ƏV 1 oV 
VV =-—[s— == =0 
s Os (sx) + Dog 
substituting our potential into the Laplace equation and multiplying every term by s? yields 


s d ( dS(s) 1 dé( s [| d?S(s)  dS(s) 1 d®(o) 
sz ms (: ds JE 








DH) de S(s) 








(>) dg? 





and letting 


S(s) 








s [( d’S(s) | dS(s) 1 dó(Q) 
— ——— ——— — — 0 
(s ae ris a pi Blo) d$ P2 Pit p2 
if we now let pı = A* and p2 = —A” „we chose these constants because we were motivated by knowing 
the angular solution should be oscillatory. We can now write 


2b sls) dS(s) 


a2 dP(>) . 
x OM age 





we can see that the angular solutions are given by 


(>) = A’el + Ble“ = Asin(Ad) +Boos(Ao) A> 1 
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where A # A’ and B Z B'. Now to solve the radial differential equation, which is a non-linear homo- 
geneus second order differential equation. We can write it in a more suggestive form 


s*8"(s) +s8'(s) —A7S(s) =0 
and now we can assume that it has a polynomial form 
S(s) 25" S'(s) 2 ms"! S"(s)— m(m— 1)5"? 


substituting this into our equation gives 


s?m(m —1)5" ?-E sms" - Mg" = 0 
m(m—1)+m—’2 = 0 
m = +A 


so we find that the radial solution is given by 
S(s) 2C -Ds^ X21 


we have only found the solutions where A > 1, what about the solutions where À — 0? We can solve 
this by solving the same differential equations for this special case. 


LO) _ 4 (459) =0 








do? ds 
for these solutions to hold we must set the first derivative equal to a constant 


d®(o) 


ie O(o) =co+k A=0 


but since we know 
&(6) — &(0--2x) 


tells us that c — 0, which gives 





P(g) =k 
where k is just a constant. Now for the radial part 
dS 
s e) =d S(s)=dinstj A-0 
s 


we have now found al! solutions to the potential which can now be written as 


V(s,0) =|k+dIns+ y [Cs + Das ™ [A sin(A$) - Bj, cos(4.)] 
A=1 


To show that our solution accomidates an infinite line charge. We know that the electric field of an 
infinite line charge is given by 


A 
S 
2TSE0 





E= 
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and the potential is 





v=- f É-ds= m 
0 2T1£0 


and since we know that by symmetry, and also the fact that a line does not any well defined angles we 
get 


[Cs + D^] [A5 sin(49) 4- Bj cos(A9)] = 0 


mM: 


A=1 


we find 
V(s,0) =k+dIns 


Problem # 5 
Find the potential outside an infinitely long metal pipe, of radius R, placed at right angles to an other- 
wise uniform electric field Eo . Find the surface charge induced on the pipe. 








Given the solution found from the last problem 


V(s,$) =k+dIns+ 3» [Cs + D4s- [Aq sin(Ad) + By cos(Ad)| 
A=1 


the surface of the cylinder is equipotential and thus we can just set the potential equal to zero. We also 
know that the potential does not go to zero at large x, far from the cylinder the potential is 


V = —Eoọox+C x-—scos0 


We need to apply the boundary conditions 


OV = 0 when s=R 
(ii) V — -—Eoscosó when s R 


From boundary condition (1) we find 
V(R,0) =k+dinR+ Y [C R^ ES D4R^^| [Ay sin(Ad) + By cos(A)| = 0 
A=1 


and from condition (ii) we find 


V(s >> R,o) =k+dIns+ Y. [Cs + Ds [A sin(AQ) + By, cos(A)] = —Eoscos 6 
A=1 
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where the term D,,s~* becomes negligable for s `> R, from conditions (i) and (ii) we can see that 
k=0 d=0 A,=0 
and also that A = 1. This gives us 
V(s,0) =Cs+Ds~')Bcos(0) 
we can now absorb B into C and D we get 
V(s,0) = (C's +D's~!) cos(6) 
applying the boundary conditions again gives us 


D' 


(i) C'R4- - 0 D'——cR 


(ii) C'scos(ó) = -—Eoscosó C'-— —Eo 


this works because we can ignore the D' /s term because it is negligable. 


n 


V(s,$) — (-àv- = cos(Q) 


To find the induced charge we can use 


Eus oV 
VOU E 
we know tha 
ƏV R? 
ae cos Q |-- E 
R2 
Oo = -—&£9cosQ -E-E s=R 
s 
o = )|2£9ÉocosQ 


Problem # 6 
A spherical surface of radius R has charge uniformly distributed over its surface with a density Q/4xR? 
, except for a spherical cap at the north pole, defined by the bounding cone 0 — a. 
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o-Q/AxR^2 


(a). Before we get started here, show the following recursion relation for Legendre polynomials to be 
true: 


dPra(x) E (x) | (21+1)P)(x) 2 0 (3.6) 


dx 
(for / > 0). This can be obatined by considering the Rodriguez formula. (Griffiths 3.62) 


The Rodriguez formula (Equation 3.62 in Griffiths) is given as 


PO) = gr (4) e-» 


taking the first derivative with respect to x on P;+1(x) yields 


d 1 d\' a Lee 
Sp an e rod esie ced eres s 24) 4 91,2 (4? — 1)/-! 
dx i0) 7 STET (=) a0 E ag) er ee 
in the above calculation we used the following relationships 
1 l 


CESR comi 





taking the first derivative with respect to x on P;_;(x) yields 


d 2 Eaa COP dA d ess 
Se (| Ee eee N A 
ax 10 7 3073 (=) ea” Sonn a Y 
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and finally expanding the last term gives 


1 


l l 
QI 1) — QLED)377 (5) (2—1 = (QD) (=) Gat) 1) 


putting all this together yields 


(2) [e-wt«aee - 7 -e- n eroe iy] = 0 


l 
(5) [ne o? - 1) - g? - n7! 2b? - 7 e 2 = 14] — 0 

x 
0 = 0 


(b). Show that the potential inside the spherical surface can be expressed as 


Q 1 r! 


V= d Tat [Pi+1(cos a) — P)_1(cos a) pri” (cos 0) 
where, for / — 0, we should substitute Pj. ; (cosa) — —1. 


Integrating every term in Equation 1 gives 


P. gi eS 
/ pude ar 


Knowing that the surface charge density and the general expression of the potential inside a sphere are 
given by 


[Pri (x) — Pr_1(x)] 





= Q 
V(r,9) — Ý Ar’ P (cos8) o= 
d 4TR? 


motivated by example 3.9 in Griffiths we know that the coefficients A; are given by equation 3.84 


1 


A 2gRl-1 


T 
i c(0)P;(cos0) sin@d@ 
0 


we also know from Fouriers trick 


b b a 
| Preosesineas = — [ nydx- | Piddi 
a a b 


We find 


1 


A 
| 2egl-1 


Qa 1 T 
[ o(8)P;(cos 6) sin@d@ + ——— f o(0)P;(cos0) sin 0d 
0 2eR! 1 Q 


the first integral is just zero because the surface charge density is zero there and the second integral is 
given by 


AI 





1 T : o(0) cos 
= zl, o(0)P;(cos0) sin0d0 = TUE I. P;(x)dx 

where the integration limits are given by the change of variable of the integral. From the relationship 
from part (a) we find 


due Q 


———— |P},,;cosa—P);_; cosa] = P;.,;cosa—P)_; cosa 
IFT. j34 js ] 141 =A | 


8re(21 - 1) RH | 
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where Pj,1(—1) — Pj-.1(—1) — 0, putting this into the general expression for the potential inside a 
sphere yields 
r' P;(cos®) 


D+ REA [Pi +1 cos & — Pj. cos ar] 


Q oo 
V(r <R,8) =| — 
(r F L RA 
(c). What is the potential outside the sphere? 
We know that the general expression for the potential outside of a sphere is given by 


a -B 
V(r>R,0) =) "1; Pí(cos0) 
1-0 


but we know that the coefficients B; (Equation 3.81) can be expressed as 


B, - AR? 

thus 
V R,8)— ———BR 0) 2|——»? ———— B ( — Pj. 1cosQ 
(r > R,®) d JH ! (cos) Bre 2S Gri rol 14.1 COS yt 





(d). Find the magnitude and the direction of the electric field at the origin. 
We know that the electric field is given by 
E=-W =-— 
or 
Since 0 is not well defined at the origin and we know that the electric field points in the 2 direction and 


we do not compute the derivative with respect to O. The electric field is given by 


Q c Ir! Pjcos 0 z 
E; = ——— » — — P cosa -— P cosalf 
nO 8ne e (i+ RP Hi at Cosel 


and so we know f = Z, and 6 = 0. We can also see that the above expression tells us the / = 1 in order 
for the electric field to not blow up at the origin or go to zero. 


E = — l [P2 cos a — Po cos q] 
= "zu (Seta 5) -11 
= a zos a- 1]2 
E = Ta pill 0s? ol 
and the magnitude is given by 
|E| pl! cos? a] 


(e). Discuss the limiting forms of the potential (b,c) and the electric field (d) as the spherical cap 
becomes either (1) very small or (2) so large that the entire remaining charged area becomes a very small 
cap at the south pole. 


We can see from the figure that if & gets very small the sphere would be almost represented as a 
spherical conductor and so 
a<<1l avery small 


the potentials and the electric field would behave as 
Vir « R,0)z Vo 
1 
V(r » R,0)«- 
F 
E(r < R,0)z0 
and if the cap becomes very large the sphere would be represented by a point charge at the south pole 
of the sphere. The potentials and electric field would behave as 
1 
V(r « R,0)e- 
F 
1 
V(r > R,O)«- 
F 


1 
E(r < R,0) = = 
r 


Problem # 7 

Many nuclei have electric quadrupole moments. The quadrupole moment of a given state of a given 
nucleus is typically quantified as Q = Q33/e where Qj; is defined as in Griffiths’ Problem 3.45, and the 
index “3” refers to the Z axis. Now, consider a nucleus centered at the origin finding itself in a cylindrically 
symmetric electric field with a gradient JE, /dz along the Z axis at the position of the nucleus. 


(a). Show that the inhomogeneous field modifies the energy of the quadrupole by 


e (OE, 
w= -§0(5) 


For this, you might want to find a suitable form of the electrostatic potential near the position of the 
nucleus. Be sure that this potential satisfies the Laplace equation. You should let V = 0 at the origin. 


We know that energy of a continuous charge distribution can be written as 
W = I p(r)e(r)av 
where ®(7’) is the externel potential energy, where the general form can be expressed as 
Ix S 


ó(r)— beu / (P, (cos0')p()dc' 


we can expend this using Taylor expansion to find 


(0) = (0) +r. V&(0) + 5 L798 (0) 
hj 
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where we know that the last term is the quadrupole term, which is the one we are concerned about. 
Plugging this into the work function yields 


1 
w-z f ptr) Pend jo (0)dv 
and since we are only concerned with Q33 we know that i = j = 3 , thus we find 
2 


wat foar LI 


but we know that 
o (0) = 9E. 
oz © Q 


putting this into the above expression yields 


= 1 0E; 2.2 / / 
w= -3 [pleat 


using the relationship from problem 3.45 


5 S Q33 5 
03 = BEP- EP > S= |p) 

using this expression we find that the work is given by 

Q33 0E; 

W = -= —— 

4 oz 

since we know that 

Q33 = eQ 


we find 
Ms B 
4 Oz 
(b) Nuclear charge distributions can be approximated by a constant charge density throughout a 
spheroidal (egg-shaped) volume of semimajor axis a and semimajor axis b. Calculate the quadrupole 
moment of such a nucleus assuming that the total charge is Ze. Given that the !*?Eu(Z = 63) nucleus has 
a quadrupole moment of Q = 2.5 x 10~7°cm? and a mean radius R = ae — 7 x 107 Pcm, determine the 


fractional difference in radius (a — b) /R. 
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Using the definition of the quadrupole moment 
Qij = [ox — (r')*8;j]p(7 dv’ 
and since we are only looking for Qz; we find 
Q..=(r) [ D —? - "spa: 


This gives us the quadrupole moment in cartesian coordinates, it would be easier to solve this integral 
in spherical coordinates, we can do a linear transformation to do this, the equation that models our spheroid 
is given by 


De udo 22 
X^ cas. Lm 
putt 
letting 
(= yee. wes 
b b a 


where the derivatives are given by 
dx —bdu dy=bdv dz= adw 


and so we know 
dxdydz = ab*dudvdw 


we also know that the cartesian components in spherical coordiantes are given by 


u = rsin@cosd 
v — rsinOsinQ 
w — rcos0 
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using the above relationships we get 


a..=p(r) | | fwa- (bv)? ~ (bu) "Jab ae 
which gives 
2n T pl 
p(r) [ f f [2a2r2 cos20 — b?r? sin? 0 cos? 0 — b?r? sin? @ sin? d]ab?r? sin d@dodr 
0 0 JO 
simplifying this gives 
2n pt prl 
Qo str) f ^ [ [2a?r? cos? @ — b?r sin? @]ab2r? sin @d@dodr 
0 JO JO 
and using Mathematica yields 


Q« — p(r) nab? (a? s p^) 


but we know that the charge density is given as 


| 4 | 34 arb au 
pu am. VAS 


subsituting this expression gives 
2 
Qz = zla —b)(a+b)) 
solving for a — b/R gives 


a—b I 50 = 5Qzz 
R 2Ze(a--b)R | AZeR? 








since we know that 
Qz: = Qe 


we find that the fractional difference is given by 


a—b 50 


i IZE ~ 0.01 








Problem # 8 

Consider a dipole with moment p at a distance zo above a grounded conducting plane that is taken 
to define the x — y plane. p is taken to lie in the x — z plane, oriented at an angle 0 to the normal (z) 
of the plane. Find the torque N on the dipole. What are the equilibrium values of 0 ? [Answer: N 
-— p? sinO cos 0$/32n£0z;, )] 
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[: : [^ - 
dipole 4. 9 dipole a 
p P 
i E 


We know that the torque from an electric field on a dipole is given as 
N=pxE 
and we also know that the energy of a dipole is given as 


o 1l 
UP Aner 
we can solve this problem by using the image of the dipole, we know that the components of the dipole 
and the image are given by 


[3(p-£)£ — p| 


Pdip = (p sin 0,0, pcos 9) 
Pim 
Pim = —Pdip 


(p’ sin, 0, —p’ cos) 


we also know that f points in the z direction, thus 
p.$—-p.Z—p'cos0 r=2z9 


putting this into the electric dipole equation yields 





1 1 
[3p' cos 02 + [p’ sin 8& — p’ cos 64]] = ——-— |p’ sin 0X + 2p’ cos 63] 


Eim 73 
321€ £3 


B Ane (2z9)? 


and thus we know that the torque is 


p $£ -f Z 
= 5 sinO QO cos@ 
3272) | sin@ 0 2cosO 


p^sinOcos0 , 
Negew— om 
32nez, 
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N=pxE 


which yields 


The equilibrium values of 0 are given when the torque goes to zero, i.e 
nc ] 
0— c ud — 0 integer 


Problem # 9 
Using Eq. 3.103, calculate the average electric field of a dipole, over a spherical volume of radius R, 


centered at the origin. Do the angular intervals first. [Vote : You must express f and 6 in terms of &, $ and 
z before integrating. Compare your answer with the general theorem Eq. 3.105. The discrepency here is 
related to the fact that the field of the dipole blows up at r = 0. The angular integral is zero, but th radial 
integral is infinite, so we really don’t know what to make of the answer. To resolve this dilemma, let’s say 
Eq. 3.103 applies outside a tiny sphere of radius €- its contribution to Egye is then unambiguously zero, 
and the whole answer has to come from the field inside the £—sphere. 


We know that the average electric field and the electric field from the dipole can be expressed as 


; Eji dt A 
E; = Ia dip — E (2oos Of + sin 08) 


we need to change to cartesian coordiantes so that we can evaluate this integral, i.e 


"mb 


= SsinOcosó& 4- sinOsinó$ -- cos 02 


> 


— cosO0cosó& -- cosOsinó$ — sin 0 
thus Ej, in cartesian coordinates gives 


ee 


T [3 sin@ cos Ocos OX + 3 sin O@cos @ sin oy + (2cos” @ — sin” 0)3] 


Egip = 


and so we find that the integral can be expressed as 


2t (Rt Ap? r (2x (Run? : 
p [f 3 sin^ 0 cos O cos ọ ~ P o ji 3 sin^ O cos Osin A 
E4;dv — — —————— — —— drd0d — —— — — — ——drd0d 
space PPS M h, r P Oe ane 0 JO 0 r í 9y 
2n fX 2 2 
p [f (2cos^0 — sin^ 0 x 
— L—— — — — — drd0d 
n] I | r rado 


using Mathematica yields (for the angular integrals) 


T Ej; ,dt = 0 
ang 


and we know that the integral for the radial component (at the origin) yields 


E i dT = œ 
2 ME 


E; —0:oo 


thus we find that 


which seems very weird? There must be a delta function in the center of the sphere. 


(b). What must the field inside the £—sphere be, in order for the general theorem (3.105) to hold? 


[Hint : since € is arbitrarily small, we're talking about something that is infinite at r — 0 and whose 
integral over an infinitesimal volume is finite.][Answer : —(p/3€0)8°(r)| 
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[Evidently, the true field of a dipole is 


o 1 l PEN l 3 
"me pow £)£ —p| — ——pó9 (r) 


E;; 


You may well wonder how we missed the delta-function term when we calculated the field back in 
Sect. 3.4.4. The answer is that the differentiation leading to Eq. 3.103 is perfectly valid except at r — 0, 
but we should have known (from our experience in Sect 1.5.1) that the point r = 0 is problematic. 


Since we know that the average electric field inside the sphere is not zero we can write the electric 
dipole as 


Egip = Ape ole E)É—p] Ker) 


where K is a non-zero coefficient that we need to solve for. we can solve for K by solving the average 
electric field in terms of a delta function 


5 _ [K&S (r)dt_ 3K 
PT V — Anm 


setting this expression to equation 3.105 we find 


3K p p 
fee ml Apc dum 
An RP 4TE0R? 3£9 


thus the correct expression for the electric dipole is given by 
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Chapter 4 


Electric Fields in Matter 


4.1 Problems and Solutions 


Problem # 1 Forces on dipoles 

(a). A dipole p with fixed magnitude is placed in an external electric field E(r) that varies with position. 
Show that there is a force on the dipole, and find how this force depends on the orientation of the dipole 
and on the variation of the field. [Answer: F = (p- V)E] 


So we know that 
F, = qE, F_ = —qE. 


and we also know that the total force is given by 
Fio = F} +F- = q(E4 — E) = qAE 
where the components of E are given by 
AE; = (VE;)-d AE,=(VE,):d  AE,—(VE;.d 
where d is the displacement vector. Thus we know 
AE = (d - V)E 
thus the force is given by 
F=(p-V)E 


(b). Find the force and torque on a dipole in the field of a point charge. Let the charge q be at the 
origin, and the dipole p — p(sinSo& -- cos$o2) be at the point (0, 0, zo ). Also, find the force on q, and 
verify Newton's third law. [Answer: the force on q is 

Pod ix 


= — sin of + 2 cos &oz 
4n£0zj ( Šo S0 ) 





We know that the potential of a point charge can be written as (in cartesian coordinates) 


q 1 


AE ASF 


Ve mE 
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and we know that the electric field can be written as 


Rely aes |g Bg 
Atte lo ds (pma (x2 - y? 222 


now that we have the electric field we can find the force 
I 0 0 
F = | psinëo—E + pcos 60 —E 
Ox Oz 


which yields 


F(0,0,z0) = Falsingos — 2cos602] 


and the torque is given by 


i in 
N(0,0,29) =pxe =| Psingo 0 posto |-|-i5hing 
0 0 dnt 0 


To verify Newton’s third law we need to show the force on q from the dipole is equal and opposite of 
the force on the dipole from q, thus we can start with the defenition for the electric field of a dipole 
1 


Egip(r) = Tep OP Êf- p] 


and we know that the force is given by 
F= qKaip 
we also know that 


^ 


f=-Z p= p/(sinEoX+cos€oz) (p-f)F = pcos Ez 


thus we can see that the force is 





Pq A ; A A Pq 
F = —|3cos€oz — sinEGox — cos Ez] = 
4nezp | Šo Šo ud 4nez 


[2 cos Eoz — sin ox] 
thus we can see that 
F = qKEaip = —qE. 


We know that the equilibrium values for 0 are given when 


9-7 n0 


Problem £ 2 Capacitors with dielectric slabs 

The space between the plates of a parallel capacitor (Figure 1) is filled with two slabs of linear dialectric 
material. Each slab has thickness a, so the total distance between the plates is 2a. Slab 1 has a dialectric 
constant of 2, and Slab 2 has a dialectric constant of 1.5. The free charge density on the top plate is o and 
onthe bottom plate —o. 
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Capacitor with two dialectrics 
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Figure 1 


a). Find the electric displacement D in each slab. 


Since we know what the free charge is on both the top plate of the capacitor and the bottom plate of 
the capacitor we can apply the integral form of Gauss's law to find the displacement, using 


D-da= Of 
top 


we can find what the displacement if we place a Gaussian pillbox encompassing the top plate of the 
capacitor and Slab 1 


Qy — ojA D;:.da— 0;A da — A(Z) 


top 


where Z points in the direction of the electric field, which is downward. So for Slab 1 we find the 
displacement to be 


D; = 6 /(2) = o(2) 


and for the displacement of Slab 2 we can place a Gaussian pillbox encompassing the bottom plate of 
the capacitor and Slab 2 


Qr — ojA D2 -da = ofA da — A(2) 


bot 
thus 


D2 = of (Zz) = —o(2) 


and the magnitudes are given by 


b). Find the electric field E in each slab. 


given that we know what D is in each slab we can easily find E 


D, Oo 
E;— — &=2e E; — —2 
1 a 1 0 1 289^ 
and D 5 
2 Oo, 
E,—— &=<€& E, = -— 
2 S 2 ze 2 38 


and the magnitudes are given by 


c). Find the polarization P in each slab. 


We know that the polarization is defined as 


P= e0X-E 
where X, is defined as 
€ 1 
—=€=E9(1+X%e) -— Xe—-£&-l1l Xil  X27z 
£0 2 
for Slab 1 and 2 we find 
o. o. 
P; = —Z P2 = ——Z 
2 3 
and the magnitudes are given by 
o o 
P=- P= = 
nn EE 


d). Find the potential difference between the plates. 
we know that the potential difference is given by the following integral 


oa 20a 


a 2a 
Av - M -vi - 1f FE, -dl— E> -dl| = — 
0 a 2€9 389 


7oa 
AV =|—— 
zl 


e). Find the location and amount of all bound charge. 








given that dl = a2 we find 


given that we know the polarization we can find what the bound charge is by using 
0, — P- n 


where fi points outward with respect to the dialectric, so for Slab 1 we find 
o 


Op = =F Slab 1 top surface 

Op = S Slab 1 bottom surface 
and for Slab 2 we find 

Op = -2 Slab 2 top surface 

Op = ; Slab 2 bottom surface 


this is illustrated in Figure 2. 


Capacitor with two dialectrics 





+0 





total-top(1) 
LIS -0/2 


Slab 1 
+0/2 
-0/3 











Slab 2 i | 
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Figure 2 


f). Now that you know all the charge (free and bound), recalculate the field in each slab, and confirm 
your answer to b. 


Using Figure 2 we can find the total charge o; for Slab 1 and Slab 2. For the top of Slab 1 we add up 
all the charge in the dashed box labeled total-top(1) which yields 


o o 
=O--=- Top of Slab 1 
0; —6 275 op of Slab 


for the bottom of Slab 1 we add up all the charge in the box labeled total-bottom(1) which yields 


6 -2-245-9--5 bottom of Slab1 


and now for the bottom of Slab 2 we add up all the charge in the dashed-dot box labeled total-bottom(2) 


which yeilds 


2 
ous 3-9 — 50 bottom of Slab 2 


and finally for the top of Slab 2 we add up all the charge in the dashed-dot box labeled total-top(2) 


which yields 


oO 0 OG 2 
E —-— lab 2 
3*5 Te 3° Top of Slab 


and the magnitude of the electric field is given by 


oOo; = —- 


_ oO o __ OF 20 


= =) 2 = = —_ 
£0 2£0 £0 3£9 


E, 


which is the same thing we found in part b. 
part 2: problem 2 


Suppose you have enough linear dialectric material, of dialectric constant €, , to hal f-fill a parallel- 
plate capacitor (Figure 3). By what fraction is the capacitance increased when you distribute the material 
as in Fig.3 (a)? How about Fig.3 (b)? 
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Dialectric 





























Figure 3 


For Figure 3 (a) we can see that the displacement is given by 
fD-da=0;=04 D=0 


this is the displacement inside the dialectric, with this we can now find what the electric field is inside 
the dialectric 


D o 
Bl Edia — = 
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we also know that the electric field between the dialectric and the top plate of the capacitor is given by 


Oo 
Fair = £o 


given that the defenition of capacitance is given by 


Cip = Vip 


we are now in a position to find what the potential is between the capacitor 


d/2 


Vus | Ead ris Egia dl = S 
0 


od od od/il 1 |. od 
= —(e-+1) 


2&9 2€ 2 \e €) 2 
thus the capacitance is given by 


| Q cA 24e 
Viz Viz  d(& +1) 


and since we know that a capacitor with no dialectric has a capacitance given by 


A£0 
C = — 
d 
thus the ratio of the capacitance with dialectric to capacitance without dialectric for Figure 3 (a) is 
given by 
Ciz 28 — | 2g 
C €g(€,+1) | (€-+1) 





From Figure 3 (b) we can see that the area has been decreased by 1/2 and the dialectric now fills this 
entire 1/2-space, thus we can treat this as two separate capacitors one in which the 1/2 area is filled with 
the dialectric and another 1/2 area that has no dialectric. We find that the total capacitance is given by 


Ae Atp> A 
Cua 7 3g 33 7 af 


and we know that the capacitance with no dialectric is given by 


£ J- £0) 


A€£0 
d 


Cy; £-£9 |£-4l 
C . 229 2 
Problem # 3 


A spherical conductor, of radius a, carries a charge Q. It is sorrounded by linear dialectric material of 
susceptibility Xe, out to radius b. Find the energy of this configuration. 


C= 


thus the ratio gives 





We know that the total work done, as we build up the free charge from zero to its final confiuguartion 


is given by 
1 
=- [| D. Ed 
X [PRU 


We can find D by using 


2. 
D-da= D= 
j da — Qr m r>a 
we know that P=D=E=0 inside the sphere. Now that we know D we can easily find E, 
D Q 
E=——=-—<~ b 
e pme |SS 
and Q 
— ——f b 
Ameor a 


now that we have D and E we can now find what the energy of this configuration is 


Q : 
w-£ [f anys am f cepe ind 


knowing that dQ — 4n when integrating over all solid angles allow us to symplify this as 


Q? E =2 j! 1 1 f1 
1, m Ad iMt zw Rs 
-ZE p ugar um Bl el 
we can also write this as 
o [/1 1 1 Q 1 
So AA rọ r—1 -— 
BTE G >) m ;|- Snc ite 5 


£y —l-cXe  £-—Eo(l-Xe) 


but since we know 


we find 
Q? 1 1 
-Ji 1 


Problem 4 
The dielectric constant of air is 1.00059. From this, determine the mean polarizability of atmospheric 
molecules. Compare this result to the atomic polarizabilities listed, for example, in Griffiths Table 4.1. 


We can write the polarization as 
P — x,£9E 


The polarization in a dialectric is also defined as 
P=np 


where p is the dipole moment and n is the particle density. We also know that the dipole moment can 
be defined as 


p= OK other 


where œ is the polarizability and Esther is the total electric field, the applied field and the enduced field 


Eother = E- Eseif 
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we need to find out what the self induced field is. The external field produced by the sphere itself is a 
pure dipole. The total field inside the conducting sphere is 0, because a conductor is an equipotential, thus 
the induced field inside must be —Ko, this field is needed to cancel the applied uniform field, so we find 


E; po —Eo 
If we treat the sphere as a conducting sphere we find that the dipole moment goes as 
p= 4n£9a^ Eo 


where a is the radius, thus E,.;¢ can be written as 


p 
Ese = aA 
self 41£9a? 
substituting this in E;;;,,, we find 
p 
E, er — E+ ——, 
‘h a Arena? 
we also know that 
E np P 
self — 3 £0 = 3 £0 
thus we find B 
P=np=na| E+ — 
np=n ( + =) 
there is a linear relationship between P and E, and the susceptibility which is given by 
= naf% 
M= TEAB 
if we solve for the atomic polarizability we find 
£0 Xe _ 3go K— 1 


~ nl+Xe/3 n x42 


and usin the ideal gas law we find 


a 3kT k—1 39 1 
— = —_——_ _ #/1.90 x 10°", 
4neo 4nPK+2 


using T = 293 K and 1 x 10°Pa for the pressure, along with the dialectric constant given. 


Problem # 5 

A hollow dielectric sphere, with dielectric constant £/£o = K, inner radius a and outer radius b, is 
placed in a uniform applied electric field Eo. The presence of the sphere changes the field. Find the field 
in the three defined regions, i.e. r « a,a « r « b and r » b. What is the field at the center of the spherical 
shell? What is the dipole moment of the dielectric medium? [for the last part, youOll find: 


. 4n£o(a? — b')b  Eo(1-- 2k) (x - 1) 
Po FB(K—1)2 —b3(2 + «)(1 +2) 


we find for a/47€ 


12 


V(3) 


Using the general solution for the potential in spherical coordinates we find thet the potentials are given 
by 


V = Y Ai Picos8 r«a 


1=0 
Ci 
V = C “rir ) Picos8 a<r<b 
1=0 


D 
y -— È -7p TPi cos 0 — Eorcos 0 
i-o" 


and the boundary conditions are given by 


( rad Baie oVi B 
lu de = Oe r=a 
. QV» E 9Va _ 
(iv) o = For r=b 


thus we have four unknown constants and four boundary conditions. Lets figure out what these con- 
stants are using the boundary conditions. From 1 we find 


C; 
IH L (so + ri) Picose 


which simplifies to 


C 
Aja! = Bya! + E I#1 |A=B+G l=1 


if 
141 A; =B,=C;=0 


from this and from the next boundary condition we know that / = 1 for all cases. From the second 


boundary condition we find (using / — 1) 
Bb -C 2 D— Eg? 
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from boundary condition number 3 we find 


2C € 
£j)À — £|B — — =— 
m | a? | £0 
which simplifies to 
2 
Anna — 2 
a 
and from our last boundary condition we find 
2D 2«KC 
Lt Eo = KB — —A- 


so now we have four equation and four unknowns, we can do this by brute force or just use use 
Mathematica. The four equations are summarized as 


C 
A = B+% 
a 
Bb’ +C = D-Egb? 
2 
E 
a 
2D 2«C 
E OR ee 
Mathematica yields 
dm 9b Eok 
2a? (x — 1)? — D$ (2 - K)(1--2x) 
ees HEU 2e) 
2a? (k— 1)? — £5 (2 -X)(1-4-2x) 
Ea 3a b’ Eo(K — 


1) 
2a3(« — 1)? — b3(2 + «)(1+2k) 
(a? —b?)b? Eg(1-:2Kk)(x — 1) 
2a? (x — 1)? — 55 (2. -- X) (1-- 2x) 


plugging these constants into the general expression for all potentials yield 


9b? Egk — 9b? Egk 


Ty SS SSS 0 
2a3(K — 1)? —b3(2++«)(1 4.2K) 2a3(K—1)2—B3(2+K)(1+2K)) 7 


V = 


and for the other two we find 
v 3b? Eo(1--2K) z 3a°b° Eo(« — 1) z 
= EEE ee eee —— ———— 
^ 2À(k-1?—P(24Xx)142x) 2a3(K — 1)? —3(2 + K) (142k) (2 +2 + 22)372 


and the final potential is given by 


(a? — b3)b3 Eo(1--2x) (x — 1) z 


= — 
2a (K—1)2-B3(2+K)(1 42K) (2@+yr+ep2 


74 


The electric field in all three regions is given by 
E = -VV 


SO 
9b? Egk 


—»]"^x-iy-5PQ4xXürzgj. '*^ 


Ei 


the electric field inside the dialectric is given by 


TEM 3b? Eo(1-- 2k) 3d 3a? b^ Eg(K — 1) X 
2  2a4(x-1?2-PQ-Xx)u-2x) 2a (k-1? - P2 4-Xx)(142x) 
| 3xz R 3yz ^ 32 i 1 N 
oe oe SO ooo sd 8 Oe 
(x - y? 72)5/2 (x? + y? +22)5/2 (x? - y? z2)5/2 (x2 - y? 4 222 


and finally for the electric field on the outside we find, using 


_ (a? —b?)bPEo(1+2«)(k—1) 
~ 2a3(«K— 1)? — B3(2+«)(1+2k) 


3AZ 3yz 32 1 


spit (xy ay - 


Bi Ss A E Lu. ar LO c ee E e cas 
3 GLEpyM 2A Qeiyisjm GRiyuu 


Eo 


The field in the center of the dialectric is given by 


MS 9b? Egk ; 
emer | 2a3(« — 1)? — b3(2 +k) (1 + 2k) 
We know that the potential in terms of the dipole moment is given by 


|| pf — p;cos0 
— 4n£gr)  4m£or? 


and looking at the potential on the outside we can see that the potential from the dipole moment is just 
the first term in the potential V5, the Eo is from the applied field, we are only interested in the potential 


from the dipole, i.e 
 |4neo(a? — b))b^Eo(1--2X)(K— 1) 
Po) F@B(K—1)2 —b3(2 +«)(1 +2) 
Problem # 6 


Consider two insulating media with dielectric constants €; /€9 = kK; and €2/€ 9 — K» , placed together 
with a planar interface between them. In the region of ?1 there is an infinite line charge with linear charge 
density Kı, parallel to the interface at perpendicular distance 














a. Use the method of images to find the electric field in both media. 
We can see from diagram (2) that the potential in the top region is given by 


/ 











$ 2TE1 z re pa 
and the potential at the bottom is given by 
A" 
Vo = In|z — 
? 2TE2 z al 


we also know that the boundary conditions are given by 
V, = Vp z=0 


from continuity. Our othert boundary condition is given as 


From boundary condition number 1 we can see that 


aec 
£2 


and from condition 2 we find 


Ac mU 


lugging in A" from the last equation into the previous one yields 
piugging q 


EE CL 


£1 
£2 


M (Scu) 
K2 + K1 


plugging this into either equation yields 


W= nal- ( 2K» Ja 
€2 +£1 K2 + K1 


76 


thus 








plugging these solutions into our potentials we find 


K2— K 1 
me-a +0 (= ) In z -- a| 


V, = 
: K5 -- Kj / 21£4 











2T1£4 


and the electric field at the top is given by 


nr (2) 1 1 
Ea a| #4 | 
2n£i|z— a| Ko +k, / 2n£i |z4- a| 


and the potential at the bottom is given by 
1 2K 
xl : J^ni-a 
2TE2 (K2 +K] 


and the electric field is given as 
1 1 1 
Ey -|-— À 
TXKoc-Ki/ |z—a| 


Show that the force per unit length on the line charge is 











dF u A? (Ki — K2) 
dL B 4n£ a(K, — K5) 
We know that 
F-—qE dF =dqE 
but 
dq — Àdl 
and since we know that there is no self induced force on A we must use the A’ term as the one con- 
tributing to the force, so we find 
dF ZZ 2 (Ki = K2) 
4TE] (Kı + K2) 


Note that the sign of the force depends on the difference Kı — K2 . Explain physically what's going on 
here. 


We can see from the expression for the force that if k; is greater than k2 the force will be pointing in 
the Z direction and if we have the inverse it will point in the -z direction. 


Problem £ 7 
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Two long cylindrical sheets of metal (radii r; and rz with r2 > r; ) are arranged coaxially. The plates 
connected to a battery that maintains a potential difference V between the sheets. The region between the 
conductors is filled with a material of conductivity o and permittivity €. 


a) Determine the capacitance per unit length of this system. 


We know that the capacitance and the potential are given by 
pu 
C—-— V--[Edl 
V 
we can find the electric field by using Gauss's law, i.e 
[£42 e da — 2xrL£ 


we find the electric field to be given as 


Q f 
2nrLe 





knowing the electric field allows us to find the potential 


"n Q Q r2 
yes dr = ——In| — 
» 2nrLeg g 2TLe : (2) 


knowing this allows us to find the capacitance per unit length 








b) Use Ohm’s law to calculate the electric current per unit length between the conducting shells. 


We know that Ohm’s law can be expressed as 
I= [J-da=o [E-da 


E-da== Q=dAL 
[Ea 


using Gauss’s law we find 


thus we find that the current per unit length is given by 


I oA 216 


= — V 
L € ln(r2/r1) 


c) Suppose the battery that maintains the potential difference V is suddenly disconnected from the 
circuit. Show that charge will leak off the two plates of this capacitor as an exponential function of time. 
What is the exponential time constant? You should obtain a simple function of € and o. 


We know that once the battery is disconnected the current will be given as 


j. 49 .. 90 
UE - 4x 
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and this is a simple integral of the form 


Ord t 
dQ == Ca 
Q Q 0€ 
which yields 
Qf = Qie *' 


and the exponential time constant is given by 
o 
k= — 
€ 


d) Show that the total energy dissipated by Joule heating as the capacitor discharges completely equals 
the electrostatic energy that was originally stored in the capacitor. 


We know that the work can be expressed as 


w= | Pat 


we know that the power given for Joule heating can be expressed as 





P—PR 
thus we find 
W= f P Rdt 
but we also know that do 
Oo o Oo o 
]2——2-Qgj—e t! I = Qie 
dt Q £ n = Q £ : 
we also know that 
v m(ẹ) 
R z= — Z> 
I 21o0L 


thus we find that the work is given by 


\ 2 B iom N2 
ge ee eres f gorge (OE qum s 
£ rij 2nxoL Jo 20 \ € rij 21oL 


this simplifies into 
Q? r2 1| Q r2 1 
W = —In{ — ] = = | —In[ — = -QV 
mE A anie Na A a" 


we know that the work from a capacitor is given by 














1 1 
W — 2CV? = -QV 
2 72 


thus we have just shown that these two expressions are equivalent. 
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Chapter 5 


Magnetostatics 


5.1 The Lorentz Force Law 


The basic problem in classical electrodynamics we have a collection of charges q1,42,q3... which are 
known as the “source” charges, and we want to know the force they exert on a charge Q. Up till now we 
were dealing with charges that were at rest electrostatics and the time has come to consider the forces 
between charges in motion. Imagine that we have two wires seperated by a few centimeters and when I 
turn on a current so that it passes up one wire and down the other and the wires jump apart. How do you 
explain this? Whatever force accounts for the attraction of parallel currents and repulsion of antiparallel 
ones is not electrostatic in nature. The force responsible for the observed phenomenom is the magnetic 
force, we can understand this using Figure 1 


wire 1 wire 2 


From Figure | we can see that a cross product will allow us to solve for the force. The magnetic force 
in a charge Q , moving with a velocity v in a magnetic field B, is 


This is known as the Lorentz force law. In the presence of both electric and magnetic fields, the net 
force on Q would be 
F — Q[E 4 (v x B)] (5.2) 


If we take a look at the Lorentz force law we will find some bizarre particle trajectories. 
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5.1.1 Cyclotron Motion 





A charged particle moving about a magnetic field in circular motion, with the magnetic field providing 
the centripital acceleration is represented in the figure above. A uniform magnetic field points into the 
page The charge Q is moving counterclockwise, the magnetic force points inward and has a magnitude 
QvB. We can write 


y? 


Feent = Macen — m — 
R 
in order to sustain circular motion 
2 


y 
Feen = F nag = QvB = Mp 


2 
ey 


Equation 5.3 is known as the cyclotron formula because it describes the motion of a particle in a 
cyclotron. We have also assumed that the particle moves in a plane perpendicular to B. If it starts out with 
some additional speed v || parallel to B, this component of the motion is unaffected by the magnetic field, 
and the particle moves in a helix. 

Another feature of the magnetic force law is 


Magnetic forces do no work 


For if Q moves an amount dl = vadt , the work done is 


which can also be written 


dWmag = Fmag -dl = Q(v x B) -vdt =0 


this follows because (v x B) is perpendicular to v, so (v x B) -v = 0. 


5.2 Currents 
The current in a wire is the charge per unit time passing a given point. Current is measured in coulombs- 
per-second, or amperes (A) 


1A=1C/s 
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A line charge À traveling down a wire at speed v constitutes a current 
I= v 


Current is actually a vector 
I= àv 


The magnetic force on a segment of current-carrying wire is given by 


Fma = | (vxB)dq= [ (vx Bat [amat 


if I and dl both point in the same direction, we can write this as 


(5.4) 





5.2.1 Surface Currents 


When charge flows over a surface, we describe it by the surface current density, K, defined by consider- 
ing a “ribbon” of infinitesimal width dl L, running parallel to the flow. If the current in this ribbon is dI, 
the surface current density is 
_ di 
dil 
this is simply the current per unit width — perpendicular — to — flow. This can also be written with a 
surface charge density 
K = ov 


and the magnetic force on the surface current is 


Ei l (v x B)oda = f (K x B)da (5.5) 





5.2.2 Volume Currents 


When charge flows thru a three-dimensional region, we describe it by the volume current density, J, 
defined by considering a “tube” of infinitesimal cross section da L, running parallel to the flow. If the 
current in this ribbon is dI, the volume current density is 


| dl 
— dal 


this is simply the current per unit area — perpendicular — to — flow. This can also be written with a 
surface charge density 


J (5.6) 


J= pv 


and the magnetic force on the surface current is 


Ens f (vx B)pdt = / (J x B)dt (5.7) 





Equation 5 can be written as 


1= | Jda = f J.da 
S S 
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which describes the current crossing a surface S. The dot product serves neatly to pick out the appror- 
iate component of da. In particular, the total charge per unit time leaving a volume V is 


fari [ina 


because charge is conserved, whatever flows out thru a surface must come at the expense of that 
remaining inside, this can be mathematically stated as 


9p 
uj 5.8 
J 3r (5.8) 
this is known as the continuity equation. For future reference, lets summarize what we know 
X )qivi ^v ( )Idl ~ ( )Kda — ( )Jdt 
izi line sur face volume 


this corresponds to q ~ Adl ~ oda ~ pdt for the various charge distributions. 


5.3. The Biot-Savart Law 


5.3.1 The Magnetic Field of a Steady Current 


dr 


The magnetic field of a steady current ius given by the Biot-Savart Law 


uo fIxr uo [EE 
B = —— — <T — 
(r) An j 7? An y^? 

The integration is along the current path, in the direction of the flow; dl’ is an element of length along 
th wire, and r’ , as always, is the vector from the source to the point r. The constant wo is called the 
permeability of free space. 

uo — 4n x 107 7N/A? 

These units are such that B itself comes out in newtons per ampere-meter (as required by the Lorentz 
force law) or teslas (T) 

IT—1N/A:m) 


The Biot-Savart law plays a role analogous to Coulomb's law in electrostatics. The inverse square law 
is common in both laws. 
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5.3.1.1 Wire segment 


The magnetic field a distance s from a long straight wire carrying a steady current / is given by 


Hol |. 
B = — (sinz — sin 
Ans (Sin 2 — sin 8) 


this is the general solution for any line segment, thus the magnetic field of an infinitely long wire would 
be given as. 
_ Hol 


IF (5.9) 


because 0) = —0; = 71/2. 


5.3.1.2 Circular loop of wire 


The magnetic fielda distance z above the center of a circular loop of wire with radius R, which carries a 
steady current / is given as 
I. GRE 
2 (R2 + 72)3/2 


For surface and volume currents the Biot-Savart law becomes 


A 


yt 7 
B(r) - ^9 T Ee i a aad a sk EE He 


4n y? 
5.4 The Divergence and Curl of B 


5.4.1 Straight-Line Currents 


The magnetic field of an infinitely straight wire is shown as 


According to Equation 7, the integral of B around a circular path of radius s, centered at the wire is 


I I 
fB a= f Ed = E TUE if 
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Notice that the solution is independent of s, thats because B decreases at the same rate that the circum- 
ference increases. Now suppose we have a bundle of straight wires. Each wire that passes thru our loop 
contributes ul, and those outside contribute nothing. The line integral will then be 


[pea = uolenc (5.10) 


where Jenc stands for the total current enclosed by the integration path. If the flow of charge is 
represented by a volume current density J, the enclosed current is 


lene = | J-da 


with the integral taken over the surface bounded by the loop. Applying Stoke’s theorem to Equation 
5.8 gives 


[ (0 xB)-da=w [ I-da 


and hence 
V x B = uoJ (5.11) 


We also know that the divergence of B is given by 
V-B=0 


the divergence of the magnetic field is 0. 


5.4.2 Applications of Ampere’s Law 


is called Ampere’s law (in differential form). It can be converted to integral form by the usual device 
of applying one of the fundemental theorems-in this case Stoke’s theorem 


[O xB)da= $B-di=m f J-aa 


Now, { J-da is the total current passing through the surface, which we call Jenc , thus 


The equation for the curl of B 


j5.a — uolenc 


This is the integral version of Ampere's law; it is generalized to arbitrary steady currents. Just as the 
Biot-Savart law plays a role in magnetostatics that Coulomb's law assumed in electrostatics, so Ampere's 
plays the role of Gauss's: 


Electrostatics: Coulomb  — Gauss 
Magnetostatics: Biot-Savart — — Ampere 


In particular, for currents with appropriate symmetry, Ampere’s law in integral form offeres a lovely 
and extraordinarily efficient means for calculating the magnetic field. 
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5.4.5 Comparison of Magnetostatics and Electrostatics 


The divergence and curl of the electrostatic field are 


V.E — P Gauss's law 
£0 
VxE — 0 no name 


These are Maxwell’s equations for electrostatics. The divergence and curl of the magnetostatic field 
are 


V-B — 0 noname 
VxB = upolenc Ampere’s law 


These are Maxwell’s equations for magnetostatics. Maxwell’s equations and the force law 
F= Q(E+vxB) 


constitute the most elegant formulation of electrostatics and magnetostatics. 


5.5 Magnetic Vector Potential 
5.5.1 The Vector Potential 
Just as V x E = 0 permitted us to introduce a scalar potential (V) in electrostatics, 
E-—-VV 
so V. B = 0 invites the introduction of a vector potential A in magnetostatics 
B=VxA 


The potential formulation automatically takes care of V - B — 0 (since the divergence of a curl is always 
Zero); there remains Ampere's law 


VxB-Vx(VxA)- V(V-A) — VA — uoJ 


if we let 
V.A—0 


with this condition to A, Ampere’s law becomes 
VA= —uoJ 


This is nothing but Poisson's equation. This is actually three Poisson's equations, on for each cartesian 
component. Assuming J goes to zero at infinity, we can read off the solution 


J r” 
awt [a 
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5.5.2 Summary; Magnetostatic Boundary Conditions 


Just as the electric field suffers a discontinuity at a surface charge, so the magnetic field is discontinous at 
a surface current. Only this time it is the tangential that changes. We can write this as 


Babove — Bbelow — uo(K x fi) 


where ni is the unit vector perpendicular to the surface, pointing “outward”. Like the scalar potential 
in electrostatics, the vector potential is continuous across any boundary 


Aabove = Abelow 


for V- À — 0 guarantees that the n ormal component is continuous, and V x A = B, in the form 


f^a- [54-9 


means that the tangential components are continuous. But the derivative of A inherits the discontinuity 


of B. 


0A 0A 
above _ %*below _ —uoK 


on on 


5.6 Multipole Expansion of the Vector Potential 


If you want an approximation formula for the vector potential of a localized current distribution, valid at 
distant points, a multipole expansion is in order. The vector potential of a current loop can be written as 


ul a 1 


A(r) - —— Y, —4 $ ()" P, (cos0)dr 
An =r 


or more explicetly 


yt fa ar zpr cosa +; je "(s cos? 0-5) a's. 


as in the multipole expansion of V, we call the first term (which goes like 1/r ) the monopole term, 
the second (which goes like 1/r?) dipole, the third quadrupole, and so on. It happens that the magnetic 
monopole term is always zeror, for the integral is just the total vector displacement around a closed loop 


fat = 


this reflects that there are no magnetic monopoles in nature. In the absence of any monopole contribu- 
tion, the dominent term is the dipole (except in the rare case where it too, vanishes) 


Lol / Mol A Wan 
Adip = Eaz fro os8'dl' = an 2f era 


this integral can be written in a more illuminating way if we invoke Eq 1.108 from Griffith's, with 


frat —— x n 


c=f: 
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then A 
uo mx fî 


where m is the magnetic dipole moment: 
m=/ f da — Ia 


where a is the area vector. The magnetic field of a pure dipole is easiest to calculate 1f we put m at the 
origin and let it point in the z—direction 





According to equation 5.10, the potential at point (7, 0,0) is 


uo msinO a 
Agip(r) 2 11— —9 


and hence 


Uom doen cue UA 
Bgip(r) =VxA= gr (2098 Of + sin 60) 


88 


5.7 Problems with Solutions 
Problem # 1 Polygonal electromagnets 


* 
N 


^ <= B(z)=B(s)cos $ 
^ S 


i 


yor 








-R 
Figure 1 Figure 2 


a) A square loop of wire, with size 2a x 2a lies in the xy-plane centered at the origin. Calculate the 
field everywhere on the z -axis produced by a current / through this loop. 


We know that a line segment produces a magnetic field a distance s away of 
uol, . Hol . 
B(s) = N—2sin0 = N— sin8 
(s) 4s 21s 
where N represents the number of segments, in our case this is just 4, so we find 


2uol 
B(s)— cro" sind 
TS 


we can see from the diagram that 


: a a a 
SO D NAE s= ya +z? 
thus we find 
pes ee 
Ts ./2a2 +22 
this gives us our magnetic field as a function of r, We know that the horizontal components will cancel 
by symmetry and we will be left with a vertical component in the 2 direction, i.e 


2uol ? 
B, = Bcosó — B5 — 2 
S 


T (a? +z?) V2a? +z? 


given that 


a 
cos 0 = — 
s 
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b) What is the line integral of the field along the z -axis from z = œ to z = —o? Could you have 
expected this result based on Ampere’s law? Explain. 


We must now find out what the integral yields 


fe d 2uola? uH dz 
Z= _ 
es T aan (a2 +22) (2a2 + 22)1/2 





we can make a substitution of the form 
z— v2atan0. dz — v2asec? 040 


plugging this into the integral gives 





2ugla? s v2asec? 040 — 2ugla? ^i sec 040 
X J—n/2 a2(1+2tan26)(2a2(1+tan26)'/2 m — J-ay2 a?(1--2tan20) 


we can now multiply both top and bottom by cos?0 to give 





2ugla? rf sec0d0 — cos?0  2ugIa? i. 2 cos 0d0 
X  J-zy2 aà?(14-2tan?0) cos?0 X  j-aj2 a? (2sin? 0 - cos20) 


knowin 
] cos?0 — 1 — sin^0 


allows us to simplify to 





2ugla? E ?  cos0d0 
m J—n/2 a2(1+sin*@) 


doing another substitution 
k —sinO dk=cos0d0 


which gives us 





2ugla? f! — dk Du. tes] 
T [zu n o e eD 
thus 
2ugla? dz 


sn $$ =] pr 
T /. (a2 -- z2) (2a? -- 22)? 
This is the line integral from —R to R where R — œ, what about the line integral of the semi-circle? 
We can solve this integral by understanding how the magnetic field behaves far way from the z axis, we 
know from the Biot-Savart law that the magnetic field goes like 


1 
Bo — 


R2 
If we let R — co we know that the integral along the semi-circle will go to zero, this exactly what we 
would expect from Ampere’s law, which just says that if you enclose a current with an infinite loop you 
should get the total current enclosed. 


c)Now find the magnetic field at the center of a n-sided regular polygon, letting the perpendicular 
distance from the center to any side be a. Show that the result approaches that for a circular loop of radius 
a in the limit n — co. 
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We know that the solution for the magnetic field of a line segment on a plane is given by 


pi Hol sin 
Ons 


B(s) 


now we need to find a relationship between the angle 0 and the number of sides of a polygon. 


Asides 4(20) — 2m 26-2 
5sides 5(20) — 2m 28-l 
6sides 6(20) — 2m == 
nsides n(20) — 2m 28-- 


using this defenition for 0 in the expression for the magnetic field yields 
uol sin(n/n) 
B(s) 2 ————— 
(s) 275 
we also know that 
: : T 
lim sin(t/n) — — 
n—oo n 


thus we find I 
Ho 

B(s) = — 

(5) 2ns 


this is only for 1 side, for n sides we find 


Hol Hol 
Rasen Se 


which is the magnetic field of a loop of current with radius s. 
Problem # 2 


Amperean loop 
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A large parallel-plate capacitor with uniform surface charge o on the upper plate and —o on the lower 


is moving with a constant speed v, as shown in Figure 3. 


a) Find the magnetic field bewteen the plates and also above and below them. 
We know from Ampere's law 


fB -da = Holenc 


where 
lene = OVL  da-2L12 


plugging this into Ampere’s law gives 








B2L — ugovL 
thus the magnetic field is given by 
B= HoOv 
B 
thus 
o 
B = -2 vo above top plate 
HoOV , 





= + y below top plate 


2 


and the bottom plate is analogous to the top plate except all the signs change due to the —o surface 


charge 








B = + “y above bottom plate 

B = a below bottom plate 
combining our results we find 

B = 0 above and below the plates 

B = woovy in between the plates 


b) Find the magnetic force per unit area on the upper plate, including its direction. 
We know that 
Fnae = 9(v x B) = [ow x B)da — [& x B)da 


where 
K; = ovx 


and we also know that the force is caused by the bottom plate, thus 











H0OV , 
B= 
2 y 
so we find that the cross product yields 
& y» £z 270 
o 
KxB=/ov 0 0 DET 2 
o 
o0 E o 
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thus the force per unit area is given by 


2:72. 
[E [(KxB)da = Hee th 








A 2 


c) At what speed v would the magnetic force balance the electrical force? 








We know that 
Fe = Finag 
but we also know 5 
Fe = du Finag = qu 
thus 
v= A > y= l = 
£00 v/uo£o 


thus we have shown that these two forces would be equal at a velocity of the speed of light. 
Problem # 3 








Amperian loop 





L 


Amperian loop 


Suppose a very long solenoid carries a current of 10 A and has 1000 turns per cm. Find the radial 
force per unit length f on one turn of the winding. Careful: the magnetic field for an infinite solenoid 
is discontinuous - how did we handle such discontinuities to determine the pressure on the surface of 
a charged conductor? Such forces pose severe design constraints on the construction of large, typically 
superconducting, electromagnets. 


We know that the force per unit length is given as 


_ Finag 
f= L 





and we can find Fingg using 
i5 fo x B)dq = fo x B)Adl = (ov x B)da — (K x B)da =1 f(a x B) 
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and since we know that K and B are orthogonal we find that the force per unit length is 


Fina o 
—— —IB 
L 
we can find the magnetic field by using Ampere's law, we also know that the magnetic field runs 
parallel to the axis and by the righthand rule we know that it points upward inside the solenoid and 


downward outside. For the loop inside the solenoid (loop 2) we find 


js : dl — BL — uonlesc — uonlL 


for the loop outside (loop 1) we find 
js : dl — BL — uonlojc — 0 


due to the fact that the loop does not enclose any current. Thus we find that the magnetic field for this 


solenoid is given by 
B-—0O outside 


B = pon! inside 
we also know that the magnetic field of an infinite solenoid is discontinous at a surface current, which 
is the same case as the discontinuity of the electric field at a surface charge. This means that the magnetic 
field is just the average of the magnetic field below and and the magnetic field above, 
B= Babove + Boelow -£ ponl, 
2 2 


thus we find that the force per unit length is 








F, mag — uon : A 
Ema 
plugging in the constants yield 
F, mag 


= |6.28— 





L 


Problem # 4 

Cycloid motion 

Griffiths’ Ex. 5.2 obtains a cycloid solution for motion of a charged particle in crossed electric and 
magnetic fields using the Lorentz force law in the laboratory frame. Read this example, and then solve the 


same problem using an alternate approach: 





94 


a) From relativity, but in the limit of velocities small compared to the speed of light, if frame S moves 
with velocity Bc with respect to (laboratory) frame S’ , the electric field is transformed as 


E — E- B x cB 


Taking B along Å and E along Z, as in the text, and assuming E < cB in magnitude, choose the simplest 
B so that E’ vanishes. 


We know that for E’ to vanish, B x cB — —E, thus 


which gives 
—cBp,2 -- cBp.$ — —E2 
so we find 


E, 
p= ox 


b) The two frames coincide at t = 0. Using the Galilei transformation for velocities (valid for B < 1), 
if the particle is at rest at t = 0 in frame S , what is its initial velocity in frame S’ ? 


we know that the Galilei transformations for velocities can be expressed as 


u =u-v 
or 
u’ ux — Bxc 
Uy — dy-— Bye 
u, = uz—B,c 
but we know that 
Bx = Bz =0 


thus we find that the initial velocity in the S’ frame is 


EE Ay 
u y= uy— Byc u=0 


c) The magnetic field transforms in this case as 


SO 


cB - cB- B x E 


Because E « cB and p < 1, approximate B’ to be equal to B. Now, in frame S’ solve for the motion 


of the particle. Using the Galilei transformation for positions, determine the motion of the particle in the 
lab frame. 
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If we assume that B’ = B than 
BxE=0 E=0 


sO we can see that 
Frag = q(v' x B’) 


and lets assume thatv' has velocity components 
v= (x,y,z) B'—(B,0,0) 
thus we find 


v xB = = Bżŷ — Bý2 


OS s o 
O w. 
OAN N 


the magnetic force is now given as 
Fj4, — q(Bz$ — Byz) — má — m(y$ 24) 
which reduces to 
qBi—my  —qBý = mz 


we can define the cylotron frequency as 


and we find 
j= z--ay 
if we differentiate the first to get 


which becomes 


y—oy-0 
this has a general solution of the form 
y(t) — Ai1cosot -- A»sinot J- A3 
ý = —A,@sinot + A2@cos Of 
ï = —A,@* cos wt — A207 sinwt 
VS A10? sin Of — AQ? cos (Of 


and we can see that this solution works for our differential equation. Now to find z/(f) we can use the 
first derivative of y' which yields 


ž = —QGy-A,;09?sinor — A;09?cosot 
—A ,Wcos Wt — A2@sin wt 
Z(t) = —Ay;sin@rt+A2cosot +A4 


V— 0e 
| 


now we have two expressions for the motion of the particle, 


y(t) = Ay;cos@r+A2sinwr+A3 
Z(t) = —A;sin@t+A2cos ot +A4 
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now we need to find what these constants are using our boundary coditions 


y'(0) 
v, (0) 


z(0)—0 boundary condtion #1 
—cp$ v.(0)—0 boundary condtion $2 


from these two boundary conditions we find 


0 — Ao2-—cA4 A2-—-A4 from 1 


W W 


thus we find that the solutions that describe this motion are given by 


1 cB. 1 cB cB 
t) = —— sin Ot t) = —— cos Ot + — 
y() - -5 zo)-- i 
Problem # 5 
6) Along with the necessary addition made by Maxwell, Ampere’s law reads 
1 OE 
c* ot 
Taking the divergence of this equation, show that electric charge is conserved globally. We thus see 
that the conservation of charge is tied fundamentally with the validity of Maxwell’s equations. 


V x B = uoJ + 


taking the divergence of both sides of this equation yields 


1d 
V: (V xB) =mV:J+ -z5 (V-E) 


but we know that 


this gives us 


HoV- J+ 25 =0 
but 
T 
H0€0 
and so we find án 
V-J= -3 


because the charge is conserved, whatever flows out through a surface must come at the expense of 
whatever is inside. Therefore this is the mathematical statement of local charge conservation. This is 
generally refered to the continuity equation. 

Problem # 6 On the gyromagnetic ratio 


A thin uniform donut, carrying a charge Q and mass M, rotates about its axis as shown in Figure 1. 
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a). Find the ratio of its magnetic dipole moment to its angular momentum. This is called the gyro- 
magnetic ratio ( or magnetomechanical ratio). 


We know that the magnetic moment is defined as 
m — yL 


where y is the gyromagnetic ratio and L is the angular momentum of the system. We also know that 
the magnetic moment can be defined as 
m — Ja — Inr? 


where J is the current and a is the area. The angular momentum of a hoop is given by 
L — Io 


where I is now the moment of inertia of a hoop, which is given by I = mr’, thus we find the angular 
momentum to be given as 
L = mro 


we also know that the current can be expressed as 


2 
1-2 gre 7-22 
T o 2x 


putting all this together we find that the gyromagnetic ratio is given by 


m Qonr? ne 
T-L-2nxm?o |2m 
b). What is the gyromagnetic ratio for a uniform spinning sphere? [This requires no new calculation; 
simply decompose the sphere into infinitesimal rings, and apply the result form part a).] 





From part a) we can see that the gyromagnetic ratio does not depend on the radius and thus decom- 
posing the sphere into infinitesimal rings will still yield the same ratio. The gyrommagnetic ratio is only 
dependent on the total charge Q of the system and the mass M. This can be expressed mathematically as 

Xm Y$5o Q 
EL Emro 2m 


€). According to quantum mechanics, the angular momentum of a spinning electron is sh, where h 
is Planck's constant. What, then is the electron's magnetic dipole moment, in A-m?? [This semiclassical 


98 


value is actually of by a factor of almost exactly 2. Dirac’s relativistic electron theory got the 2 right, 
and Feynman, Schwinger, and Tomonaga later calculated tiny further corrections. The determination of 
the electron’s magnetic dipole moment remains the finest achievment of quantum electrodynamics, and 
exhibits perhaps the most stunningly precise agreement between theory and experiment in all of physics. 
Incidentally, the quantity (eh/2m), where e is the charge of the electron and m is its mass, is called the 
Bohr magneton. | 


since we know that the magnetic moment is defined as 
m — yL 
and the angular momentum of a spinning electron is given by 
1 
L=-h 
2 
we find the magnetic moment of an electron to be given as 


m = Ou —|4.61 x 107” A m? 
4m 
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Chapter 6 


Magnetic Fields in Matter 


6.1 Magnetization 


All magnetic phenomenom are due to electric charges in motion, if you could examine a piece of magnetic 
material on an atomic scale you would find tiny currents: electrons orbiting around the nucleus and elec- 
trons spinning about their axes. For macroscopic purposes, these current loops are so small that we may 
treat them as magnetic dipoles. Ordinarily, they cancel each other out because of the random orientation 
of the atoms. But when a magnetic field is applied, a net alignment of these magnetic dipoles occur, and 
the medium becomes magnetically polarized, or magnetized. 

Unlike electric polarization, which is almost always in the direction of the field, some materials acquire 
a magnetization parallel to B (paramagnets) and some opposite to B (diamagnets). A few substances 
called ferromagnets, retain their magnetization even after the external field has been removed. 


6.2 Torques and Forces on Magnetic Dipoles 


A magnetic dipole experiences a torque in a magnetic field just as a electric dipole does in an electric field. 
Lets calculate the torque on a rectangular current loop in a uniform magnetic field B. 








We have centered the loop at the origin, and a tilt of angle 0 from the z axis toward the y axis and the 
magnetic field B points in the z dircetion. The forces on the two sloping sides cancel(they tend to stretch 
the loop but not rotate it). The forces on the "horizontal" sides are likewise equal and opposite(so the net 
force on the loop is zero), but they do generate a torque: 


N = aF sin®X 
The magnitude of the force on each of these segments is 


F = IbB 
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and therefore 
N = JlabBsinO& — mBsinO& 


m 


where m — Iab — IA is the magnetic dipole moment of the loop. This gives the exact torque on any 
localized current distribution, in the presence of a uniformfield; in a noninuform field it is the exact 
torque (about the center) for a perfect dipole of infinitesimal size. Notice that this is analogous to the 
torque given in electrostatics. In particular, the torque is again in such a direction as to line the dipole 
parallel to the field. It is this torque that accounts for paramagnetism. Paramagnetism normally occurs 
in atoms or molecules with an odd number of electrons, where the “extra” unpaired electron is subject to 
the magnetic torque. 
In a uniform field, the net force on a current loop is zero: 


F=1f (axe) =1( fat) xB=0 


the constant B comes outside the integral, and the net displacement $ dl around a closed loop vanishes. 
For an infinitesimal loop, with dipole moment m, in a field B, the net force is 


or 


F = V(m x B) (6.2) 


which is analogous to the its electric counterpart F = V(p x E). 


6.3 Effect of a Magnetic Field on Atomic Orbits 


Electrons not onl;y spin; they also revolve around the nucleus- for simpliciuty, let’s assume the orbit is a 
circle of radius R. Although technically this orbital motion does not constitute a steady current, in practice 
the period T = 27R/v is so short that unless you blink awfully fast, its going to look like a steady current 


e ev 
I = — Z — 
T  2nxR 
Accordingly, the orbital dipole moment (ITR?) is 
m = -3eR 


(The minus sign accounts for the negative charge of the electron). Like any other magnetic dipole, this 
one is subject to a torque (m x B) when the atom is placed in a magnetic field. The orbital contribution to 
paramagnetism is small. There is, however, a more significant effect on the orbital motion: The electron 
speeds up or slows down, depending on the orientation of B. For wheraa the centripetal acceleration v?/R 
is ordinarily sustained by electric forces alone 


1 e v? 


np RD UR 
in the presence of a magnetic field there is an additional force, —e(v x B). For the sake of argument, let's 
say that B is perpendicular to the plane of orbit, as shown in figure 6.3, then 


2 2 


e E v 
+ evB = neg 


4T1£0 R2 
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Under these conditions, the new speed v is greater than v 


Me 





evB = (PF v) S (v 4- v)(v — v) 
or assuming the change Av — v — v is small 
eRB 
Apc 
á 2Me 


When B is turned on, the electron speeds up. 
A change in orbital speed means a change in the dipole moment 


eR 


4m, 


B (6.3) 





1 
Am = — 5 el Av) Ra =— 


Notice that the change in m is opposite the direction of B. (An electron circling the other way would 
have a dipole moment pointing upward, but such an orbit would be slowed down by the field. 


6.4 Magnetization and the Field of a Magnetized Object 


In the presence of a magnetic field, matter becomes magnetized. Upon a microscopic examination it will 
be found to contain many tiny dipoles, with a net alignment along some direction. We describe the state 
of magnetic polarization by the vector quantity 


M = magnetic dipole moment per unit volume 


where M is called the magnetization; it plays a role analogous to the polarization P in electrostatics. In the 
following section, we will not worry about how the magnetization got there- it could be paramagnetism, 
diamagnetism, or even ferromagnetism- we shall take M as given, and calculate the field this magnetization 
itself produces. 


6.4.1 Bound Currents 


Suppose we have a piece of magnetized material; the magnetic dipole moment per unit volume, M, is 
given. What field does this produce? Well, the vector potential of a single dipole m is given by Equation 


5.12 ] 
uo m x fî 


r) = ——— 

(r) 4n r? 

where r’ is the seperation vector. In the magnetized object, each volume element dt’ carries a dipole 
moment Md", so that the total vector potential 


A(r) Li [ete 


— 4x 7"? 


That does it, in principle. But as in the electrical case, the integral can be cast in a more illuminating 
form by exploiting the identity 
E. 


pt. 
Vcn 
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with this, 


A(r) - D f Me’) x (vi) dv 


Integrating by parts, uisng product rule 7, gives 


A(r) - £2 |f iv x M'a- fv x ka ar! 


r! 


We can write the latter as a surface integral 


A(r) = m J ly x M(r)Ja — ^ | i [M(r^) x da/] 


r! 4r) r 
The first term looks just like the potential of a volume current 
(64) 
while the second term looks like the potential of a surface current, 
K,—M xf (6.5) 


where fi is the normal vector. With these definitions, 


po f Ji) s Ed / 
A- P f od etm | da 





What this means is that the potential (and hence also the field) of a magnetized object is the same as 
would be produced by a volume current Jp = V x M throughout the material, plus a surface current K, = 
M x fion the boundary. 


6.5 The Auxiliary Field H 


6.5.1 Ampere's law in Magnetized Materials 


In the previous section we found that the effect of magnetization is to establish bound currents J, = V x M 
within the material and K, = M x fi on the surface. We are now ready to put everything together: the 
field attributable to bound currents, plus the field due to everything else- which I shall conductocall free 
current. The free current might flow through wires imbedded in the magnetized substance or, if the latter 
is a conductor, through the material itself. In any event, the total current can be written as 


J=Js+Jşf (6.6) 


In view of Equation 6.4 and 6.6, Ampere’s law can be written 


1 
ee de SNM, 


or, collecting together the two curls 


Vx (—5-) =J; 
Ho 
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The quantity in paranthesis is designated by the letter H: 


1 
H = —B -M (6.7) 
Ho 
In terms of H, then, Ampere’s law reads 
VxH-J; (6.8) 
or, in integral form, 
fima rae (6.9) 


where If is the total free current passing through the Amperian loop. 

H plays a role analogous to D in electrostatics: Just as D allowed us to write Gauss's law in terms of 
the free charge alone, H permits us to express Ampere's law in terms of the free current alone- and free 
current is what we control. Bound current, like bound charge, comes along for the ride- the material gets 
magnetized, and this results in bound currents 


6.5.2. Boundary Conditions 


The divergence of H is given by 
V.H— -V.M (6.10) 


The magnetostatics boundary conditions can ne written in terms of H and the free current. From Equation 
6.10 it follows that 


ij L Nu 3 L 
H bove = Hbelow SSM above Mbelow) 
while Equation 6.8 says 


| | 
H bove E Hbelow 


In the presence of materials these conditions are sometimes more useful then the corresponding boundary 
conditions on B 


=Ky xf 


ie jp _ 
B above . B below . 0 
And I 
Babove Y Boelow = Ho(K x fi) 


6.6 Linear and Nonlinear Media 


In paramagnetic and diamagnetic materials, the magnetization is sustained by the field; when B is removed, 
M disappears. In fact, for most substances the magnetization is proprtional to the field, provided the field 
is not too strong. 

M — y,,H (6.11) 


The constant of proportionality Xm is called the magnetic susceptibility; it is a dimensionless quantity that 
variues from one substance to another- positive for paramagnets and negative for diamagnets. Materials 
that obey the above relationship are called linear media. We also know 


B — uo(H - M) — uo(1 - X,,)H 
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for linear media. Thus B is also proportional to H 
B — uH 


where 
uU — uo(1-- Xu) 


u is called the permiability of the material. In vacuum, where there is no matter to magnetize, the suscep- 
tibility Xm vanishes, and the permiability is uo. That's why uo is called the permiability of free space. 


6.7 Problems and Solutions 


Problem # 1 Forces on magnetic dipoles 


A uniform current density J = Joz fills a slab straddling the yz plane, from x — —a to x — —-a. A 
magnetic dipole m = mpoX is situated at the origin. 


Amperean loop 





a). Find the force on the dipole , using Equation 1. 


The force on a magnetic diple is given as 
F = V(m- B) (6.12) 


We need to use Gauss’s law and an Amperean loop to find what the magnetic field is, we can see that 
the magnetic field will point in the $ dircetion, We know it cannot point in the Z direction or the X direction. 


[ 9:4 |da 


we know that dl = y$ and da — xyz, thus 


B = uoJox$ 
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and so we find that the force is 
F = V(m-B) =[0] 
because m and B are orthogonal. 
b). Do the same for the dipole pointing in the y—direction: m = mo. 


We find that the force with this new dipole orientation is given by 


i ~ : - 
F = V(m: B) = V (mof : uoJoxĝ) = > (mouoJox) — 


€). In the elesctrostatic case the expression F = V (p - E) and F = (p- V)E are equivalent (prove it), 
but this is not the case for the magnetic analogs (explain why). As an example calculate (m - V)B for the 
configurations in (a) and (b). 


We know that in electrostatics this statement is true 
V(p-E) =(p-V)E 


we know that we can expand the left hand term using Equation 4 in the from of Griffith’s, suing product 
rules, i.e 
V(p-E) =px(VxE)+Ex(Vxp)+(p-V)E+(E-V)p 


but we know 
px (V x E) =Ex(V xp) = (E-V)p=0 


the first term goes to 0 by definition, if E is a uniform constant electric field, the second term goes to 
0 because p is also a contstant vector, in a linear dialectric medium p is proportional to E. The last term 
also goes to 0 because p is a costant vector, and when you take the derivative of a constant vector you get 
0. Thus 
V(p-E)=(p-V)E QED 


The magnetic analog is given as 
V(m-B) = m~x (V x B)+ Bx (V xm)+(m-V)B+(B-V)m 


but we also know that 
Bx (V xm) = (B-V)m=0 


and 
m x (V x B) 40 


thus we find 
V(m-B) =m x (V x B)+(m-V)B 


for the configurations in a) we find 


9 E - 
(m-V)B= mo (uoJox)$ = 


and for the configuration in b) we find 
9 n 
(m. V)B = mos, Hodox)¥ - [0] 
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The next set of problems concern the magnetic trapping of ultracold atoms, specifically atoms of "Rb 
(one of the rubidium isotopes). In one of their internal states (the F = 1 hyperfine state), such atoms have 
a magnetic dipole moment of m = mg/2 where mg = eh/2m, is the Bohr magneton. Magnetic trapping 
is achieved under the condition that the atomic magnetic moment remains anti-aligned with the magnetic 
field as the atom moves around in space, thus producing a trapping potential of the form U (r) — m|B |. 


Problem # 2 Spherical quadrupole magnetic trap 


First, let’s consider a simple kind of magnetic trap known as the spherical quadrupole magnetic trap. 


Inverse Helmholtz, Coil Quadrupole Field 


a J 


ME «y)=00,0.9) 
à à ya) 


I 








a) Consider two electromagnet coils sharing the same axis and separted by d, each of radius a and 
N turns, each coil carrying current J running in the opposite direction as that of the other coil. Placing 
the origin at the midpoint between the two coils, calculate the magnetic field B(x, y,z) in the region near 
the origin, keeping just the lowest order terms. Note that to determine the field for radial displacements 
(i.e. in the x or y directions), you can consider the symmetry of the problem and apply the condition of 
zero-divergence for the magnetic field. 


An expansion about the origin of the magnetic field B(x, y,z) yields 
B(x,y,z) = Bot f'(B:) =ax+by+cz 


By(x, y, z) Bo f'(By) — ax by  c'z 
B,(x,y,z) Bo+ f (Bx) — ax - by 4 c"z 


we know that Bo is O at the origin. We also know that 








a=b=b=d =c'=b"=c"=0 
thus we can write the expanded field as 
B(x,y,z) ^ a"x - by - cz 


and since we know that V - B — 0 we find 
ð y 9 9 
— —b' — =0 
34 +t T (y) + = e(2) 
but since we know that this configuration has rotational symmetry we can say 


IH / 
a —b 


3,09 = ay) 


we find 
c= —2a" 


now we just need to find out what c is. From our expansion we can see that 
f' (Bz) ^ e 


and we know that the magnetic field from a loop of wire at some point z away from the loop is given 
by 
| Hol a? 
72 ay 


but for our configuration, and the fact that we have two loops with currents in opposite directions we 


find 
| Hol a = ee eee 
Ee e cm e| 


B(z) 


and we know that 


f _ Hol 3a^(z — d/2) 3a? (z - d/2) 
PEO) E G+ (4/22 A Gina] 


but since we want to know what the magnetic field is at the origin we find 


; 3ugla? d 
fae) | aa 


thus we find that 


= 3ugla? | d | 
2 | (a2 + (d/2)?)>/ 
and : 
a! opo oS Sula aa 
2 4 | (a+ (d/2)?)°? 


thus the magnetic field at the origin is given by 


d 
| (xX + y¥ — 222) 


z + (d/2)?)5/? 





b) You should find a spherical quadrupole field with field gradients in the radial and axial directions 
determined by properties of the coils which create the field. A spherical quadrupole field can be written as 


B(x, y,z) = B'(x& - y$ — 222) 


where B’ is the radial field gradient. Now put in some real numbers. Consider coils of radius a = 2.5 
cm separated by d — 2a/ v/3 (the optimal separation), each of 20 turns running J = 92 amperes. What are 
the axial and radial gradients (G/cm)? 


We know that the radial field gradient is given by 


3uoI Na? d 
Bo = T Z | hern = 1561G/cm 


(a? + (d/2)?)°/? 
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and the axial field gradient is 


3uoINa? d 
/ == 
B ial m 2 ha USER =| —3120G/cm 


Problem # 3 Majorana losses 


In the spherical quadrupole magnetic field described above, the potential U(r) = m|B| is clearly mini- 
mized at a point, and, thus, particles with their magnetic moment anti-aligned with the field will be trapped 
by such a field. However, a problem with the spherical quadrupole trap is the fact that it is always lossy. 
Magnetic trapping requires that the atomic magnetic moment always follow the local orientation of the 
magnetic field. If the field is changing its direction faster than an atom’s orientation can keep up, the 
atomic orientation may change, causing the atom to be ejected from a magnetic trap. This type of loss is 
given the name Majorana loss. 


a) Consider a "Rb atom moving with velocity vand impact parameter b in the x — y plane of a spherical 
quadrupole trap. Give an estimate for the critical impact parameter bg below which the atom has a high 
probability of being lost, and above which the atom has a low probability of being lost. To do this, consider 
the maximal rate at which the orientation of the field is changing - this quantity has units of s~!. Compare 
this with the Larmor precession frequency of the atom in the local magnetic field (this frequency is zero at 
the origin and larger elsewhere). 


We know that the orbital frequency and Larmar frequency for this problem are given by 


" 
Qorp— — Qra, — yB 
bo 


setting these two expressions equal to each other yields 
bo 


but we know that 
Mm ^m 
thus we find 


2mv 
QB' 


b) Now we will estimate the Majorana loss rate for atoms at a certain temperature trapped in a spherical 
quadrupole trap. To do this, consider that the treatment above determines an ellipsoidal volume within 
which atoms at a given velocity may be lost from the trap. A loss rate can be calculated as the product nov 
where nv is the density of the atomic gas, © is the area of the ellipsoid (which differs slightly if approached 
from different directions; just take some characteristic area), and v is the typical velocity of atoms. To 
estimate the volume occupied by the gas (you'll need this to determine the density n), consider how large 
is the volume in which the trap potential U (r) is smaller than kpgT where kg is the Boltzmann constant. 
Similarly, to obtain a typical velocity, consider the velocity at which the kinetic energy is kgT . Having 
made all these estimates, find the loss rate per atom given the temperature T and the radial field gradient 
B. 





bo — 


Since we know that the lost rate is given by 


Lrate = NOV 


109 


and we know that the energy is given by 
U (r) 2 m|B| 2 m|B'(x& - y$)| 2 mB'r 


and if the thermal energy is equal to the potential energy we get 





kgT — mB'r 
rearranging this for r gives 
| kgT 
mB’ 
and the volume is given by 
V 2 3 a 
=ne —— 
3 a 
but since 
bo <r 
we find the volume to be given as 
4 
V = irr 
3 r 


and using our expression for r we find that the volume is given by 


4 f(kgT V? 
V = -T| — 
3 eon ) 
and we know that the loss rate is given by 
N 
Lrate = NOV = yo 
we know that the velocity can be found using 


1 2kpT 
“mv =kgT v=4/—— 
2 m 


and also 
o= nb? 


3b5\/ P 3Nup?(B') 


thus 


Lrate = 





4 (E) Q(kaT? 


c) Again, let's put in some numbers. Consider a gas of *Rb trapped in a spherical quadrupole trap 
with a radial gradient of 120 G/cm. What is the Majorana loss rate per atom at temperatures of 1, 10, and 
100 uK. 


Using 





Lyate = 3us (B^? 
N il 54 
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we find 
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Problem # 4 The TOP trap 


Losses from spherical quadrupole traps made evaporative cooling of atoms in magnetic traps ineffec- 
tive until some manner was employed to avoid these losses. As you found above, the problem becomes 
ever worse as the temperature of the gas gets lower, thus precluding the ultra-low temperatures needed 
to reach Bose-Einstein condensation. One solution, employed by the team of Cornell and Wieman in 
Colorado for the first achievement of Bose-Einstein condensation in a gas, was the TOP (time-averaged 
orbiting potential) trap. Consider that a spherical quadrupole trap is formed, and then an additional bias 
field is applied of the form 


Brop(x,y,z) = Bocos@t& + Bosin@ty 


Atoms held in a time-varying magnetic trap of this type will experience a time-averaged potential at 
every point so long as (1) the field rotates at a rate faster than the motional time-scale of atoms so that an 
atom experiences a local average of the magnitude of the field, and (2) the field rotates at a rate slower than 


the Larmor precession frequency so that atomic magnetic moment remains anti-aligned with the magnetic 
field. 


a) Assume that these conditions are met. What is the time-averaged magnitude of the magnetic field 
near the origin? You should find that this magnitude has a minimum at the origin, and varies quadratically 
away from the origin with a fixed ratio of trap curvatures in the different directions. 


We know that the total magnetic field is given by 


Bro = Bi (x& + y¥ — 22%) + Bocoswt& + Bosinwty 
= (B'x+Bocosmrt)&+ (B’y + Bocos wt)¥ — 2B’ z% 


and the magnitude is given by 
Brot = \/ (B'x)? + (B’y)? + (2B'z)? + 2B'Bo(xcos ot + ysinar) 


thus magnitude of the time-averaged magnetic field is given by 


o 22/o 
(Brot) = "x j (B'x)? -- (B'y)? 4- (2B'z)? -- 2B'Bo(xcos wt + ysin wt )dt 
This can also be written as 
o [27/9 
(Brot) = zf Bov 1+ edt 
2T Jo 
where € is given as 
BW 2B' 
€= | — (x +y +427) + —(xcos ft + ysin of) 
Bo Bo 


111 


we can now Taylor expand about the origin to find 
2 


€ € 
MNES EARS T y e 


and we know that 


2 

Ba 2B! 

g£ —|[—) (2 42) 9 — (xcos or + ysinar) 
Bo Bo 


and we know that the only thing that survives from this in the time averaged field is 


B! * 

Dek, 

£^ £z E (xcosar + ysinar) 
Bo 


where the first higher order terms in x, y, and z vanish in €?. We can write the time averaged expanded 
field as (after integration) 





b) Now put in some numbers. Suppose on top of the 120 G/cm radial gradient spherical quadrupole 
field, the time-varying TOP field has a magnitude of 10 G and rotates at a rate of 0 — 21 x 7.5 kHz. What 
are the oscillation frequencies of a "Rb atom in this trap? (Note: if a one dimensional potential has the 
form U(x) — (1/2)mQ?x? , then a particle of mass m will oscillate harmonically in this potential with 
oscillation angular frequency 2). Do the field settings chosen satisfy the two conditions for adiabaticity 
described above? 


We know that the potential energy in three dimensions can be written as 
L O23 
U(r)= xin r^ — upgB 


where B is the time averaged magnetic field from part a). Thus we find that the oscillation frequency 
can be expressed as 





B! B l 
Q, = (=) Beet 
0 





c) Finally, the critical temperature T, for Bose-Einstein condensation for atoms in a harmonic trap is 
determined by the geometric mean of trapping frequencies à and the number of atoms N by the relation 


3 
kpT. 
N=6&(3) | —— 
«e (E) 
where kg is the Boltzmann constant and € is the Riemann Zeta function (€(3) 21.202). This equation 
can be interpreted as stating that once the number of quantum states available to a gas at a temperature 
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T. roughly equals the number of particles, multiple occupancy of quantum states occurs and quantum 
statistical effects such as Bose-Einstein condensation ensue. Consider the above TOP trap with 2 x 104 
atoms. What is the Bose-Einstein condensation transition temperature? 


The critical temperature can be expressed as 


$ NI hà 
"^ kge(3)!/3 


where Q is defined as the geometric mean 


1 


à — (0,0,0,)!/5 — 213- 


thus 


T ~ [LEk] 


Problem # 5 Refrigerator magnets and atom mirrors 


Most of you have probably seen large, flexible refrigerator magnets. These are flat magnetic materials, 
up to a couple of inches on a side, that stick to an iron refrigerator door by inducing magnetization in the 
door with the magnetic fields emanating from the surface of the magnet. 


a) Suppose you have such a magnet with a circular shape, and that magnet is uniformly magnetized, 
with a magnetization out of the plane. What is the magnetic field on the surface of the magnet at its center, 
as a function of the radius of the magnet? On the basis of this observation, explain why this is not how 
refrigerator magnets are magnetized. 


Since we are assuming that this is a ring of charge then we already know what the magnetic field is, 
we have seen this problem before, the magnetic field of a ring of charge is given by 


2 
Bi) = _* _, 
2 ( R24 E /2 

but since we are only concerned about the magnetic field at z = O we find that the magnetic field is 
given by 
_ bol, 
~ 2R 

This would not make a very good magnetic due to the fact that as the radius increases the magnetic 
field decreases. If we say that the thickness of the magnetic is h and also a magnetization of M we can 
write the current as 


B(R) 


I = Mh = K,h 
and find that the magnetic field due to this configuration is given by 


uoMh, — uoKph, 
B(r) = 27 2— T. 2 








and this field has the same charactersistics. Thus if we have a magnet that is twice as big, the magnetic 
field will be twice as small. 


The magnetization in such a refrigerator magnet is actually a periodic linear array of oppositely ori- 
ented magnetization. Let us consider the field produced by such a magnetization. Specifically, consider a 
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magnet that fills the space —d < z < 0, being infinite in extent in the € and $ directions. Let the magneti- 
zation M be oriented in the plane of the magnet, varying as 


M — Msinkx& 


where M gives the maximum value of the magnetization and k the spatial wavevector of its modulation. 























b) Show that this magnetization can be represented by bound surface currents at z — O (top layer) and 
z — —d (bottom layer). 


We know that 
Kp =Mxn 


so for the top surface 
K,(z=0) =Mxz=Msinkx(-y) ñ=? 
and at the bottom 
K,(z= —d) =Mxz=Msinkx(¥) ñ= -—ĉ 


c) Consider just the field due to the top layer of bound currents. In the region z = 0, there are no 
currents, and therefore V - B — 0 and V x B — 0. Thus, we can define a scalar magnetic potential V; for 
that obeys V?y,, — 0 and B = —VV,. Considering these relations, the boundary condition at z — oo, and 
the y-independence of this problem, show that Vm can be written as 


Vin(t) = x Vocos(kx 4- 9)e ^ 


where the — refer to solutions for V,, either above or below the layer of currents. The values of Vo and 
Q are still to be determined. Hint: recall solutions to the Laplace equation that are separated by Cartesian 
coordinates. 


We can solve this problem by using seperation of variables in cartesian coordinates 
Vm =X (x)Z(z) 


we need to solve Laplaces equation 


QV, OV, 
v2 De m m — 
V, 32 za 0 
which becomes y y 
X Z 
— =k SEZ > 
ox? (x) Oz (2) k20 


which is just a second order differential equation with a decaying solution and an oscillating solution, 
ie 
X(x)—Acoskx--Bsinkx | Z(z) — Ce E De * 


and our boundary condition z — œ Vm — 0 gives 
C=0 Acos(kx) +Bsin(kx) = A‘cos(kx +6) 


where we find 
X(x) =A’ cos(kx +0) Z(z) =De™ 
and we can finally write the scalar potential as 
Vo cos(kx +)e*|_ above the surface 
~Vocos(kx+)e* below the surface 


Vm 
Vm 


where Vo absorbs A’ and D. 


d) Now obtain an expression for the magnetic fields from the top layer of currents, still containing 
unknown quantities Vo and 9. Use Ampere’s law to determine those quantities and thus to fully determine 
the fields. 


We know that the magnetic field can be found with 
B(z > 0) = —W yn 


which can also be written as 


B(x,z > 0) = — (Ses ea) = kVo(sin(kx + o)e“& + cos(kx + o)e 2) 


and for below the surface we know that the magnetic field can be found with 
B (z d 0) = —Win 


which can also be written as 


B(x,z <0) =— (Ses 4 Pa) = kVo(—sin(kx + o)e”& + cos(kx + o)e2) 


to find what Vo and $ is we can use the integral form of Ampere’s law, which is given by 


n : dY — uolenc 


but we know that dl = dx$ on the top of the loop and Tenc = Kpdx thus we find, we know that the 
verticle line segments cancel, thus we find 


l l 
2| Vo(sin(--6)& — keos(Er--0)2) dii =p | M sinkxdx 
0 0 


which becomes 
2kV J (sinkx + o)dx = pM / sinkxdx 
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thus we can see that 
LoM 
W — — — 0 
pce D 


thus the magnetic field on the top layer is 


M 

Bix,z > 0)= Se (sin(kx)&+ cos (kx)8) 
M 

B(uz « 0)= Se (= sin(kx)& + cos (kx)8) 


e) So far you have determined the field from the top layer. Now add the fields from the bound currents 
on both layers of the magnet. What is the magnitude of the field on the surface of the magnet (z = 0)? 


We know that the contribution to the magnetic field from the bottom layer is going to be equal and 
opposite as for the top layer (current is running in oppoiste direction) than the magnetic fields for the 
bottom layer are given by 


M 

B(x,z > 0)= — Be (sin(kx)8 + cos(kx)2) 
M 

B(x,z < 0)= — Fe (~sin(kx)& + cos(kx)2) 


To find the total magnetic we must add these two 


M M 
Br =B,+Br= T sin(kx)& — cos(kx)2) — ST (sin (kx) — cos(kx)2) = 0 


Aside from binding to refrigerators, such periodic magnetizations can also be used as mirrors for beams 
of atoms with magnetic moments. See, for instance, Roach et al., Physical Review Letters 75, 629 (1995). 
2) 


Problem # 6 The magnetic scalar potential 


Maxwell’s equations for the magnetic field (Ampere’s law, in particular) clearly demonstrate that the 
magnetic field cannot be defined via a scalar potential. Further, there are no magnetic charges, and, hence, 
no magnetic charge density. Nevertheless, such quantities are useful even though they are solely fictional. 
Here we consider the use of the magnetic charge density pm and the magnetic scalar potential V, as devices 
to determine the fields from magnetized materials. 


a) We start with the field of the ideal magnetic dipole. In class, we found out that the electric field of an 
ideal electric dipole and the magnetic field of an ideal magnetic dipole are quite similar in form, the only 
difference being the singular field at the position of the dipole. Considering the relation between these two 
fields, argue why one can write the field from a magnetic dipole as 


B(r) = -VV,(r) -- uoM(r) (6.13) 


and determine the magnetic scalar potential V,, in this case. 


We know that since electric properties are analogous to magnetic properties, then the defenition of the 
magnetic dipole seems reasonable. We know that we can compare the first term of the magnetic dipole to 
the first term of the electric dipole. The second term is neccassery so that (like in electrostatics) we can 
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account for the discontinuity at r = 0 just like in the case of electrostatics. This is only in the absence of 
free currents. 


Since we know that the magnetic field and electric field have similar characteristics, we can see that 
the electric field dipole can be expressed as 


p NOS: p 
Ei = $ mA 
dip Tes cos 6f + sin 60) Pe 


we can compare the magnetic dipole to this dipole to find 





uoM(r) analogous to P sy) 
and so di 
A (2cos0f --sin00) analogous to a (2cos OF + sin 86) 
thus we can see from this relationship that 


uom A R ^ 
—Win = 45,3 605 0f 4- sin 00) 


and since we know that the scalar potential for the electric field dipole is given by 
1 f- cos 6 
Vaäip(r,0) = — -E =? 


— 4mn£) r? 4TEor? 
and the analogous formula for the magnetic dipole scalar potential is given by 


uof-m  |uomcosO 


b) Now consider the field produced by a magnetization M. Using Eq. 1, and following our treatment 
of the electric field from polarized objects, show that the magnetic scalar potential in this case would be 
defined as 








Vor) — o f M ya 


— 4x Js |r- r'| 





4% pm(r') 3 


4x Jv |r— r'| 
Here V denotes the volume of the magnetized object, S its surface (with normal vector fi), and py = 
—V. M and oy — M - f. 





Using the expression for our scalar potential 





uo f.m / 
= — r=r-r 
" 4n D? 
and the total potential is given by 
£-M(r 
4nJy r 
observing that 


we have j 
v= a M. V! (+) dt 
4n Jy r 


and integrating by parts and using product rule number 5, gives 


v-2 fv (Sav [Sv pan 


or using the divergence theorem 
H 1 |; ,. M0 1 1 1 
V = — ọ -M.da +— J -(—V -M)dī 
4n $ r x 4n Jv x ) 
and we can define 
Qu —-V.M oyv=M-i 


thus we find 


Vr) 


uo (Ow(r) uo f pu(r) 3 
=o da 4- — d 
b ete oF n v [Ir — r!| 3 





The reason why this analogy works so far is that there are, as of yet, no free currents in the systems 
we're considering. Now we add those currents back in, and we find, of course, that the magnetic field can 
no longer be given just as the gradient of a potential. 


c) Give expressions for the B and H fields in the presence of both a magnetization M and a free current 
density J, still using the magnetic scalar potentials described in parts (a) and (b). 


In the presence of free currents we must write the magnetic field as 
B(r) =| —VVm + uoM (r) + uoĦ 
where we know that 
V x B = uV x M + uV x H = uo(Jp + Jf) = uoJ 


and we also know 
Vx H = " 
this can also be written as 1 
J= —(V xB) — Jp Jf 
Ho 
Problem # 7 


In solving this problem, keep in mind the boundary conditions for B and H. On a historical note, such 
cavities were considered in the early days of electromagnetic theory as operational ways of measuring the 
B and H fields inside magnetic materials. 








Suppose the field inside a large piece of magnetic material is Bo, so that Ho — (1/uo)Bo — M. 


a) Now a small spherical cavity is hollowed out of the material. Find the field at the center of the 
cavity, in terms of Bo and M. Also find H at the center of the cavity, in terms of Ho and M. 


We know that the magnetic field for this configuration is given by 
B = Bo- B' 


where B' is the magnetic field induced by the magnetization of the material, which was given in 
example 6.1 


B’ = -uoM 
3/0 
thus the magnetic field is given as 
2 
B = Bo — 740M 
and the H field is given by 
H = Ho — H’ 
where i : 
H’ = —-M=--M 
Ho 
thus we find 


1 
H=Ho+ 3M 


b) Do the same for a long needle-shaped cavity running parallel to M. 


We know that we cann treat this problem like a solenoid problem, where we have found that the 
magnetic field is given by 
B’ = uon! = oK, = mM 


and so the magnetic field due to this configuration is given by 


and the auxillary field is given by 


H = Ho — H’ 
where u 
H——-M-0 
Ho 
thus 


c) Do the same for a thin wafer-shaped cavity perpendicular to M. 


We know that 





B= By —B’ 
but B’ can be found by assuming this is a loop of wire, which we have solve in problem 1 
2R 
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but 
K, =Mxfi= Msin6 


where if we take the limit that this thin wafer gets thinner 
h—0 B'—0 


thus we find 


and for the H field we find 
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Chapter 7 


Electrodynamics 


7.1 Electromotive Force 


7.1.1 Ohm’s Law 


To make a current flow, you have to push on the charges. How fast they move, in response to a given 
push, depends on the nature of the material. For most substances, the current density J is proprtional to 
the force per unit charge, f 

J — of (7.1) 


The proportionality factor o (not to be confused with the surface charge) is an emperical constant that 
varies from one material to another; it's called the conductivity of the medium. For most purposes metals 
can be regarded as perfect conductors with © = co. 

In principle, the force that drives the charges to produce the current could be anything- chemical, 
gravitational, or trained ants with tiny harnesses. For our purpose, it's usually an electromagnetic force 
that does the job. In this case Equation 7.1 becomes 


J 2 o(E-4 v x B) (7.2) 
Ordinarily, the velocity of the charges is sufficiently small that the second term can be ignored 
J =oE (7.3) 


Equation 7.3 is called Ohm's law, though the physics behind it is really contained in Equation 7.1, of 
which 7.3 is just a special case. Ohm’s law can also be written as 


V=IR (7.4) 


The constant of proportionality R is called the resistance; it’s a function of the geometry of the arrange- 
ment and the conductivity of the medium between two electrodes. Resistance in measured in ohms © : an 
ohm is avolt per ampere. 

For steady currents and uniform conductivity 


1 


V-E=—v-J=0 (7.5) 


and therefore the charge density is zero; any unbalanced charge resides on the surface. It follows, in par- 
ticular, that Laplace’s equation holds within a homogeneous ohmic material carrying a steady current, so 
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all the tools and tricks of Chapter 3 are available for computing the potential. As a results of the collisions 
between the electrons moving through a medium, the work done by the electric force is converted into 
heat in the resistor. Since the work done per unit charge is V and the charge flowing per unti time is /, the 
power delivered is 

P=VI =FR (7.6) 


This is called the Joule heating law. With Z in amperes and R in ohms, P comes out in watts (joules per 
second). 


7.1.2 Electromotive Force 


There are two forces that drive current around in a circuit: the source f,, which is ordinarily confined to 
one portion of the loop (a battery, say), and the electrostatic force, which serves to smooth out the flow 
and communicate the influence of the source to distant parts of the circuit: 


f=f,+E (7.7) 


The physical agency responsible for f, can be any one of may different things. Whatever the mechanism, 
its net effect is determined by the line integral of f around the circuit 


z-jra-jr.a (7.8) 


(Because ¢ E- dl = 0 for electrostatic fields, it doesnt matter whether you use f or f,.) £ is called the 
electromotive force, or emf, of the circuit. An electromotive force can also arise if a loop of wire is 
moving through a magnetic field, in this case the emf is defined as 

d® 

E£ = —— 7.9 

T; (7.9) 
This is known as the flux rule for motional ems, where ® is the flux of B through the loop and is defined 
as 

D= |B-da (7.10) 


Apart from its delightful simplicity, it has the virtue of applying to nonrectangular loops moving in 
arbitrary dircetions through nonuniform magnetic fields; in fact, the loop need not even maintain a fixed 
shape. The flux rule is a nifty short-cut for calculating motional emf’s. It does not contain any new physics. 
Ocassionally you will run across problems that cannot be handled by the flux rule; for these one must go 
back to the Lorentz force law itself. 


7.2 Electromagnetic Induction 


In 1831 Micheal Faraday reported on a series of experiments, including three that can be characterized as 
follows 


Experiment 1. He pulled a loop of wire to the right through a magnetic field. A current flowed in the 
loop 


Experiment 2. He moved the magnet to the left, holding the loop still. Again current flowed in the loop 
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Experiment 3. With both the loop and the magnet at rest, he changed the strength of the field (he used 
an electromagnet, and varied the current in the coil). Once again, current flowed in the loop 


The first experiment, of course, is an example of motional emf, conveniently expressed by the flux rule 
Equation 7.9. I don’t think it will surprise you to learn that exactly the same emf arises in experiment 
number 2, all that really matters is the relative motion of the magnet and the loop. Indeed, in light of 
special relativity it has to be so. Faraday had an ingenious inspiration 


A changing magnetic field induces an electric field 


It is this “induced” electric field that accounts for the emf in Experiment 2. Indeed, if (as Faraday found 
emperically) the emf is again equal to the rate of change of the flux, 


d® 
£= p E-dl= —-— 7.11 
f dt va 
then E is related to the change in B by the equation 


f E-dl EH d (7.12) 

.dl= — | —-da ] 
ot 

This is Faraday's law, in integral form,. We can convert it to differential form by applying Stoke’s 


theorem 


oB 
VAS ae (7.13) 


Note that Faraday’s law reduces to the old rule ¢ E- dl = 0 in the static case (constant B) as, of course, it 
should. 

In Experiment 3 the magnetic field changes for entirely different reasons, but according to Faraday’s 
law an electrif field will again be induced, giving rise to an emf —d®/dt. Indeed, one can subsume all 
three cases into a kind of universal flux rule: 

Whenever (and for whatever reason) the magnetic flux through a loop changes, an emf 


d® 
dt 


Will appear in the loop. 


7.3 Inductance 


Suppose you have two loops of wire, at rest. If you run a steady current /, around loop 1, it produces a 
magnetic field B;. Some of the field lines pass through loop 2; let P2 be the flux of B; through 2. You 
might have a tough time actually calculating By, but at a glance at the Biot-Savart law 


HO dl, X r 
B, — —1;  —— 
Pam. p? 


reveals one significant fact about the field; it is proportional to the current I,. Therefore, so too is the flux 
through loop 2: 


&) = [By -da; (7.14) 
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Thus 
®, = Mr], (7.15) 


where Ma; is the constant of proportionality; it is known as the mutual inductance of the two loops. 
There is a cute formula for the mutual inductance, which can be derived by expressing the flux in terms 
of the vector potential and invoking Stoke’s theorem 


9;— [Bid — [ (5 A): dai $ A1 -dh 


the vector potential can be written as 





uol fdh 
A; = — $ — 
4n r 
and hence ; " 
Holy 1 
®, = —— — }-dl 
2 4m f(g r' ) : 
Evidently 
Ho dl dl 
My = — : 
up (16 


This is the Neumann formula; it involves a double liune integral. It is not very useful for practical 
calculations, but it does reveal two important things about mutual inductance: 


1. Ma; is a purely geometrical quantity, having to do with sizes, shapes, and relative positions of the 
two loops. 


2. The integral in Equation 7.15 is unchanged if we switch the roles of loops 1 and 2; it follows that 


M2; =M\2=M 


Suppose now that you vary the current in loop 1. The flux through loop 2 will vary accordingly, and 
Faraday’s law says this changing flux will induce an emf in loop 2 


| dà» _ dli 


E> = ——— = pete 
dt dt 


(7.17) 
What a remarkable thing: Everytime you change the current in loop 1, and induced current flows in loop 
2- even though there are no wires connecting them. 

Come to think of it, a changing current not only induces an emf in any nearby loops, it also induces 
an emf in the source loop itself. Once again, the field (and therefore also the flux) is proportional to the 


current 
=LI (7.18) 


The constant of proportionality L is called the self-inductance (or simply the inductance) of the loop. As 
with M, it depends on the geometry (size and shape) of the loop. If the current changes, the emf induced 


in the loop is 


dI 
£= -L— 7.19 
di (7.19) 


Inductance is measured in henries (H); a henry is a volt-second per ampere. 
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7.3.1 Energy in Magnetic Fields 


It takes a certain amount of energy to start a current flowing in a circuit. The work done on a unit charge, 
against the back emf, in one trip around the circuit is —£ (the minus sign records the fact that this is the 
work done by you against the emf, not the work done by the emf. The amount of charge per unit time 
passing down the wire is J. So the total work done per unit time is 


dW dl 
—_ = —€] = LI— 
dt dt 


If we start with zero current and build it up to a final value I, the work done is 


we l LP (7.20) 


It does not depend on how long we take to crank up the current, only on the geometry of the loop (in the 
form of L) and the final current /. We can also write the work done in terms of the magnetic field as 


1 
W= z) B^d« (7.21) 
2uo Jall space 


In view of this result, we say the energy is "stored in the magnetic field,” in the amount (B? /2uo) per unit 
volume. 


7.4 Maxwell's Equations 


We are now ready to write down all of Maxwell’s equations (with corrections), they are summarized as 


(i) V-E=? (Gauss’s law) 
(ii) V-B=0 (no name) 
71:22 
(iii) VxE--2E (Faraday's law) d 
(iv) V x B — uoJ 4 uoo E (Ampere's law with Maxwell's correction) 
Together with the force law 
F — q(E--v x B) (7.23) 


they summarize the entire theoretical content of calssical electrodynamics (save for some special properties 
of matter). Even the continuity equation 


op 
Vf 7.24 
J 3: (7.24) 
which is the mathematical expression of conservation of charge, can be derived from Maxwell's equations 
by applying the divergence to number (iv). I have written Maxwell’s equation in the traditional way. I 


think it is logically preferable to write them as 


(i) V-E=2 (Gauss’s law) 
(ii) V-B=0 (no name) 

1.2 
(iii) VxE+ a =0 (Faraday’s law) ee) 


(ivy) VxB— Leo — uoJ  (Ampere's law with Maxwell’s correction) 


with the fields (E and B) on the left hand side and the sources (p and J) on the right. This notation 
emphasizes that all electromagnetic fields are ultimately attributable to charges and currents. Maxwell’s 
equatins tell you how charges produce fields, the force law tells you how fields affect charges. 
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7.5 Maxwell's Equations in Matter 


Maxwell’s equations in the form of 7.25 are complete and correct as they stand. However, when you are 
working with materials that are subject to electric and magnetic polarization there is a more convenient 
way to write them. In terms of free charge and currents they are summarized as follows 


() V.D — p; (iii) vxg--2 (7.26) 


(ii) V-B = 0 (iv) vxu-p 2 


Some people regard these as the "true" Maxwell’s equations, but they are in no way more "general" than 
7.25; they simply reflect a convenient division of charge and current into free and non-free parts. And 
they have the disadvantage of hybrid notation, since they contain both E and D, both B and H. They must 
be supplemented, therefore, by appropriate constitutive relations, giving D and Hb in terms of E and B. 
These depend on the nature of the material; for linear media 


P=€0X¥-E and M=y,,H 


SO B 
D=cE and H=— 
u 


where £ = £ọ(1 + Xe) and u = uo(1 +Xm). Incidentally, you'll remember that D is called the electric “dis- 


placement”; thats why the second term in the Ampere/Maxwell equation (iv) is called the displacement 
current. 


7.6 Boundary Conditions 


In general, the fields E,B,D, and H will be discontinuous at a boundary between two different media, or at 
a surface that carries charge density © or current density K. The explicit form of these discontinuities can 
be deduced from Maxwell’s equations 7.27, in their integral form 


(i) f pda = Ofenc (iii) fe-a=—-5 [Baa (7.27) 


(ii) $B-da = 0 (iv) f ma nace t. [0-da 
S P tis 


where (i) and (ii) are defined over any closed surface S, and (iii) and (iv) are for any surface S bounded by 


the closed loop 7. For (i) we find 
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for (11) we find 


for (iii) we find 


and for (iv) we find 


Equation 7.29-7.32 are the general boundary conditions for electrodynamics. In the case of linear 
media, they can be expressed in terms of E and B alone 


(i) EL —e,Et = op (ïi) E!-E! =0 (7.32) 
1.5 4 

Gi) BŁ-BŁ}=0 = 0 (iv) —H ——H!=k;xa 
Hai M 


In particular, if there is no free charge or free current at the interface, then 


(i) sE} -eÆ = 0 (iii) El-E!=0 (7.33) 
1 l 1 

(ii) BŁ-B}=0 = 0 (iv) —H -—Hl-0 
HU 1 M 


7.7 Problems and Solutions 


Problem # 1 
Notice the following parallel: 
VxD = 0 VxE=0 g&E=D-P 
VxB = 0 VxH=0 &H=B-uUM 
thus, the transcription D — B, E — H, P — uM, €9 > uo turns an electrostatic problem into an analo- 


gous magnetostatic one. Use this observation together with your knowledge of the electrostatic results, to 
rederive 


a) The magnetic field inside a uniformly magnetized sphere. (Eq. 6.16) 
The magnetic field of a uniformly magnetized sphere is given by Equatio 6.16 
2 
B = -uM 
3/0 


The electric field of a uniformly polarized sphere is given by Equation 4.14 


1 
E = ——P 
3£9 
if we let 
E—H P-uoM £o go mH =B- uM 
we find 


1 1 
H--2M mH =- -mM 
3 Ho 3/0 
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which simplifies to 
1 
B — uM = —34M 


thus 
2 
B = -u9oM 


b) The magnetic field inside a sphere of linear magnetic material in an otherwise uniform magnetic 
field. (Prob. 6.18) 


From example 4.7 we know that the electric field inside a linear dialectric sphere is given by Equation 


4.49 
3 


IFE 





Eo 
If we let 
E-H ¢=14+%e—-€m=14+%m Eo— Ho 


thus we have 3 





H= H 
2 Xm | 
if we let 
P = £ox,E — uox5H — uM 
therefore 
M — 4, H 
we also know 5 
oH =B—poM —H=——— 
uol 1 F Xm) 


and now to find what Ho is we must use 
UoHp = Bo — UM =Bo M=O0 outside 
thus 
Ho 


and so we find that the magnetic field inside a sphere of linear magnetic material in an otherwise 
uniform magnetic field is given by 


Ho 


3(1+%m) 
B= ^m p 
aur 


c) The average magnetic field over a sphere, due to steady currents within the sphere. (Eq. 5.89) 


We know that the average magnetic fields over a sphere due to steady currents is given by Equation 


5.89 5 
uo zm 
Baye — — — 
ave An R3 
we can derive this using our knowledge of what the average electric field due to a similar electrostatic 
problem. The average electric field is given by Equation 3.105 as 


1 p 


Eigse A 
avg 4n1£o R? 
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we also know that 


where we know that the volume is given by 


thus we can rewrite the average electric field as 


1 
Eaye — ——P 
ave 360 


which is identical to part a) of this problem which yielded 


2 
Bave = 3M 


but we know that this can also be written as 


M=— >B = 2 
m 
V ave 3y "o 


where this can be written as 


Problem # 2 


A metal bar of mass m slides frictionless on two parallel conducting rails a distance / apart (Figure 1). 
A resistor R is connected across the rails and a uniform magnetic field B, pointing into the page, fills the 
entire region. 





a) If the bar moves to the right at speed v, what is the current in the resistor? In what direction does it 
flow? 


Since we know that 


M 
| 
zx 
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from Ohm's law, where V is now the induced emf, which is given by 


| do 


V=£=-— 
dt 


and ® is the magnetic flux which is given by 
p= / B -da = Bxl 


where the normal are vector is into the page. The emf is given by 


dx 
£ — — — Bi — —Bi 
di E 


thus the current is given by 
Blv 


I= acs counter-clockwise 
b) What is the magnetic force on the bar? In what direction? 


We know that the Lorentz force is given by 
F5 — q(v x B) = qvB 


but since we know that 
qv — Il 


we find 
w?Bg? 





Fg = IIB = — (x) 
there is a magnetic force in the —X direction. 
C) If the bar starts out with speed vo at time ¢t = 0, and is left to slide, what is its speed at a later time t? 


Since we know that 





dv 
Fz = ma = m— 
3 dt 
we can solve for the velocity, i.e 
vl? B? dv 
=, = m— 
R dt 


which can be written in a more suggestive way 


t |2p? vd 
apap 


0 Rm vo y 
which is a simple integral with a solution of 


B PB 


—kt 
v = voe k= —— 
Rm 


d) The initial kinetic energy of the bar was, of course, Imv,. Check that the energy delivered to the 


resistor is exactly imve. 
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We know that the energy is given by 


E= | Pat 


where P is the power that is obtained from the loss of kinetic energy of the bar, the power is defined as 


ai = 7 ydg 


P=PR= 





where k has been previously defined, so we find that the total energy delivered to this transistor is 


EOS aaa IPB? Rio 5 3 
Qs = R opp" 2" 


Problem # 3 


Two tiny wire loops, with areas a; and a», are situated a displacement r apart (Figure 2). 





ad : : ers 


a) Find their mutual inductance. [Hint : Treat them as magnetic dipoles, and use Eq. 5.87.] Is your 
formula consistent with Eq. 7.23? 


Using Equation 5.87 
uo 1 


Bgip = 4, 52 (m-£)f —m] 


using Biot-Savart law and the definition of the flux thru loop 2 from the magnetic field due to loop 1 
o,-— Ji B; -daz 
we can define the mutual inductance of the two loops as 
Py = Moh 


where M»; is the mutual inductance. Using this we can now find 
$ 1 
Mz; = — = — J Bı -da 
21 n UR f 1 :4a2 


given the definition of the magnetic field due to a dipole we find 


LEN 


— —— — |3(m;ı -ĉ)f—mı| -a 
= Fag Bom ) 1] -a2 
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knowing 

mı = al, m? = ah 
gives us 
zl 
EP TE 
we can also show that the is the same if we look at it from the perspective of the other loop, 


M» [3(ai -f)(az-f) — al - a] 


oi 1 
My =—=-— | Bo-d 
12 b al 2:08] 


given the definition of the magnetic field due to a dipole we find 


1 uo 1 Lx 
Mij5— L2, OCDE mU] ‘al 
knowing 
mı = al; m = aol 
gives us 
1 
2 = Fe 5 (ao) (ai -f) — 85: ai] 


which is excatly the same, thus we can conclude 


Mzn =M =M 


b) Suposse a current 7; is flowing in loop 1, and we propose to turn on a current h in loop 2. How 
much work must be done, against the mutually induced emf, to keep the current /; flowing in loop 1? In 
light of this result, comment on Eq.6.35. 


Since we know that the total work done per unit time is given by 


—=-£] 
dt 
where £ is the induced emf and Z is the current going thru the loop. We need to find 
dW = —-€\hdt 


but we know that the induced emf on 1 is given by 


d® 
Est ©, = [ By-da; = Mh 
dt 
and so we find Jii 
2 
E£ = M— 
j dt 


putting this into our equation for the work yields 
dW =Mhdh >W =Mhih 
We found the mutual inductance from part a), substituting this in gives us 


1 7 i 
= mp : f) (azh : £) 5 a; ` azh] 
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where we have defined 


which gives us 





this looks very similar to the interaction energy of two magnetic dipoles. 
W=-U 
where U is the interaction energy of two magnetic dipoles 


uo 1 A 4 
U = —— |m; -m — 3(m; - f) (m? -Ê 
45 3m mo —3(m; -£) (mg -F)] 
Problem # 4 
The current in a long solenoid is increasingly linear in time, so that the flux is proportional to t :® = æt. 
Two voltmeters are connected to diametrically opposite points (A and B), together with resistors (R and 
R2), as shown in Figure 3. What is the reading at each voltmeter? Assume that these are ideal voltmeters 
that draw negligable current (they have huge internal resistance), and that a voltmeter registers JPE : dl 
between the terminals and through the meter. [Answer : Vi = QR1/(Rı + R2); Vo — —aR»/(R| 4- R2). 
Notice that V; 4 V2 , even though they are connected to the same points!] 


— [(ind) 






Ps 


I 






Solenoid 


Since we know that these resistors are in series we can write the current as 














— Vind = Vind 
Rr Ri +R 
but we know that the induced emf is given by 
d® 
Vina = E = S -~Q 
thus the current running thru this loop is given by 
= Q 
|» Ri Ra 
We can find what the potential is across resistor 1 using 
e aR, 
V. =- f E.dl— —IR; = 
1 ; 1 R +R 
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and to find the potential across resistor 2 we use 











b R 
v= | Bedi ire e 
a Ri + Ro 
thus 
oR, OR? 
V = 2=— 
Ri +R2 Ri +R2 
Problem # 5 
Suppose 
l q x 
E(r,t) = “ee E B(r,t) =0 


(the theta function is defined in Prob. 1.45b). Show that these fields satisfy all of Maxwell’s equations, 
and determine p and J. Describe the physical situation that give rise to these fields. 


If we look at 


Using product rule #5 
V-(fA) =f(V-A) +A: (VF) 


where f = 0(vt —r) and A= $ , thus we find 


ES loon) (v5) rE Veo) =? 


Equation 1.99 in Griffith’s gives 


and 


l => 


-(V(O(vt—r)) = = (380-0) = 


Since we know that 


Thus we find that 





The second Maxwell equation is 


The third equation is 


Since B = 0, we find 
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For the last Maxwell equation 
9E 
V x B = uoJ FPES 


We know that 





VxB=0 
Thus i 
J = 790. 
We can find 
gE d q [ de d(w-r) 5 
mu mec D Pme ————. | = ——vd(0t — 
or 4n£or? or He 4nr? Ls —r) dt Anl Mp 
Thus 
m M. = 
J= TR vé(vt — r) 
Problem # 6 


Assuming that ““Coulomb’s law” for magnetic charges (qm) reads 


F= Ho Im, Im " 
4n r 
where 7’ is Griffith’s script r, which is the seperation vector. Work out the force law for a monopole gin 
moving with velocity v through electric and magnetic fields E and B. 


Looking at the paper written in 1961 by Robert Katz from Kansas state university we find that Maxwell’s 
equations can be written as 











oB 

V-D = p VxE=-— (Fa) 
ot 
oD 

V-B = Pm VxH= (P+) 
ot 

he has defined the following properties 
4nr? 4nr? Arr? Anr2 


set up in this way magnetic poles satisfy the continuity equation; we demend the conservation of poles. 
We can immideatly write the magnetic analog of the Lorentz force equation as 


1 
Fin = q5;[H - v x D] 2 q,(B - v x B) 2 q4(B— 5v x E) 
c 


where the 1/c? was added in order to make the units work out, and also due to the motivation that E — cB. 
We have shown why these are the units needed in the next part of this explanation. 


We claim to find this using simple unit analysis and Equations 7.43 from Griffith's. We know that the 
Lorentz force for electric charges is given by 


F; = q,[E - v x B] 
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Thus by analogy we should find 
Fm = qm[B € ci(v x E)] 


we need to find a constant c, that would make the units work out. What are the units of qm? We know that 
the units of the force must be in Newtons and we know that the units of the magnetic field are 


B Fin = qmB 
— MEE — 
A m — Qm 


where the units of qm must be (Am). Now looking at the second term we find 
Fin = qmc1 (V x E) 


If 
Nm 
qm œ> Am vxE—> AS 
Thus the units of cı must be 
2 
S 


airs 
m2 


In order to have the force be in Newtons, a good choice for this units would be the speed of light c, thus 
Fin = qm[B + : v x E] 
now for the sign we must look at Equations 7.43 in Griffith's, we can see that 


OB OE 
— VxB= è eo— 
3: x Lode + Ho OX 


we can see that when V x B — V x E we gain a minus sign and thus we can conclude that 


V x E = —-uoJa — 


1 
EF, =m\B-—v XE 
In| zV x E] 


Problem # 7 


In a perfect conductor, the conductivity is infinite, so E=0 (for J = oF), and any net charge resides 
on the surface (just as it does for an imperfect conductor, in electrostatics). 


a) Show that the magnetic field is constant (OB/dt = 0), inside a perfect conductor. 
We know that 


oB 
VxE=-— 
* ot 
Since we know that E = 0 we find JB 
— =0 
ot 


This will only be true if the magnetic field is constant. 
b) Show that the magnetic flux through a perfectly conducting loop is contant. 


We know that jé 
£-- [&-a- - 57-0 
dt 
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Because E = 0, and this will only be true if ® is a constant. 


A superconductor is a perfect conductor with the additional property that the (constant) B inside is 
in fact zero. (This “flux exclusion” is known as the Meissner effect.) 


c) Show that the current in a superconductor is confined to the surface. 
Using Maxwell’s equation 
9E 
VxB = Hod - uo£o- =0 


and since we know that E = 0 than 
J=0 


and any charges that are moving must reside on the surface. 


d) Superconductivity is lost above a certain critical temperature (T), which varies from one material 
to another. Suppose you had a sphere (radius a) above its critical temperature, and you held it in a uniform 
magnetic field Boz while cooling it below 7;. Find the induced surface current density K, as a function of 
the polar angle 0. 


Using Example 5.11 in Griffith's we find that the magnetic field inside of a rotating magnetized sphere 
is given by 


2 
B= 3!nore 


This is the electric field before the temperature falls below the critical temperature T > Te. Once we 
decrease the temperature below the critical temperature we know that the magnetic field inside the sphere 
must be zero due to the “Meissner effect" and thus we know that we must induce a surface current that 
will oppose the external field. Thus 


Bexternal = Bot Binduced = —Boz 


Thus we know that 
2 x ` 3 Bo 
B = -ugoroZ— —Bo2 coro = —-— 
3 2 uo 


and since we know that the surface current is defined as 


K — o(0 x r) — orosin 00 


Kcu s 09 
2 Lo 
Problem # 8 


A familiar demonstration of superconductivity is the levitation of a magnet over a piece of supercon- 
ducting material. This phenomenom can be analyzed using the method of images. Treat the magnet as 
a perfect dipole m, a height z above the origin (and constrained to point in the z direction), and pretend 
that the superconductor occupies the entire half-space below the xy plane. Because of the Meissner effect, 
B — 0 for z < 0, and since B is divergenceless, the normal (z) component is continuous, so B; — 0 just 
above the surface. This boundary condition is met by the image configuration in which an identical dipole 
is placed at —z, as a stand-in for the superconductor; the two arrangements therefore produce the same 
magnetic field in the region z > 0. 


and using the result above we find 
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a) Which way should the magnetic dipole point (+z or —z )? 
We know that the image dipole must point in the —z direction so that at z = 0 the perpendicular 
component of the magnetic field must be 0. Mathematically 


m=mMZzZ mMm? = —moZ 


where m, is the original dipole and mz is the image dipole. 


b) Find the force on the magnet due to the induced currents in the superconductor (which is to say, the 
force due to the image dipole). Set it equal to M; (where M is the mass of the magnet) to determine the 
height h at which the magnet will “float.” [Hint : refer to Prob 6.3] 


We know that the force on a dipole due to a magnetic field is defined as 
F = V(m-B) 
since this is an image problem we know that this can also be written as 
F = V (m; -B2) 
which says that there is a force on the original dipole to to the magnetic field from the image dipole, which 
is defined as Lo 


Bo = Ane [3 (m5 : £)f m;] 





we know that 
mj — mi m» — —mj2 f=? 


we can now write the magnetic field from the image dipole as 


Ho 








By = =~ [-3(m2-2)2 +m 2] = — mid] 
Thus we find the force to be given as 
Ho ( mi 
F—--—Vi|-— 
2T r? 
Thus we find ; 
m 3uom1 
21^ 
we know that 
rEZ 


which is just the speration distance between the two dipoles, thus the force is given by 


3224 s 





If we let z = h we find that the height is found to be 
1/4 
h= 3uom; / 
321M, 
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c) The induced current on the surface of the superconductor (the xy plane) can be determined from the 
boundary condition on the tangential component of B (B — uo(K x 2). Using the field you get from the 


image configuration, Show that 


3mrh A 
Kult. s c 
mGA, 


where r is the distance from the origin. 





rhat 





We can find what the induced current if we know what the magnetic field is on the xy plane, thus we 
have to find what the total magnetic field is at z = 0 this can be achieved by figuring out what the magnetic 
field is due to the two dipoles, we know that the magnetic field from a dipole is given by 


Ho 2a 
Bgip — az? Dn £)£ — m] 


thus for dipole 1 which I will consider my real dipole we find 





Ho AKA 
= 3(m,-r,)r, —m 
1 anri (m; Â 1] 
and like wise for the image dipole 
H0 AURA 
B5; — ——3(m5-5)65-m 
2 m (mo - £j) — m] 


thus we know that the total magnetic field is given by 


o [1 xeu 1 EE 
Br = B, + Bo = al = [3(m) Ñ Ñ —m]] + -= [3 (m2 -Ñ )É —m;| 
4n |r; r> 
since we know that 
T] — r2 m = -mM = —m]ıĉ 


thus we can write the above expression as 








o3 ^ AX ^ ^ ^ 
Br = m; : Yj ) E — (ni, - £5)Y 
T Anci 1: Ej )É, — (má, - 5) E] 
or we can also write this as 
r= PLG -eÀ 
4nri 


we can see from the diagram that 
m,- Fr; = sinO@f—cos0zZ my, r = sin Of + cos 02 


or from the diagram on the right we also find 


A 


Z.Y; — —cos0 Z2.Y5»—cos0 
thus the total magnetic field can now be written as 


3uom| a Ne 
Br; =— An? [cos 0(f; 4- )] 





but we know that 
rı +r = 2a = 2sinĝ6f rj—wv r? +h? 


and since we can see from the diagram that 


cos§ = — = 5 sinO — iu r 
ri vyr^th* ri Vre+he 
using this we now find that the total magnetic field is 
3uoh 
2n(r? +h?)5/2 
To find out what the surface current is we use 
B — uo(K x 2) 


in order to solve for K we must take the cross products of both sides give respect to 2 we find 
2xB—2x(Kx2) 

using the following relationship 

we find 


we know that 


because they are orthogonal. We know that 


f ô 2 
2xB= 0 0 1| %Hohrm ; 
|... 3pgohrm 0 0 2n(r? -++ h2)5/2 
2r(r? +h?) 
thus we find that the surface current is given by 
3hrm E 

K — ——————À— 

2n(r?2 + h2)5/2 


Problem # 9 


140 











Figure 1: A parallel plate capacitor is connected on one side to a voltage source that increases the voltage 
across the capacitor linearly in time. 


Figure 1 shows a parallel plate capacitor that is connected at one end (x — 0) to a source whose voltage 
increases slowly and linearly with time as dV /dt — k. Assume that edge effects are negligible in this 
problem, i.e. that a,b > s. 


a) Find the surface current K in the top capacitor plate (positive voltage) as a function of position x 
along the capacitor. 


Since we know that 








Q(t) =CV 
then dQ(r) us 
t 
= C— = Ck 
dt dt 
thus 
f=Ck 
the capacitance from a parallel plate capacitor is given as 
EG 
c= 
s 
thus the current is given as 
xb€ok 
‘om 
s 
In order to satisfy the boundary conditions we must write this as 
_ (a—x)beok 
s 


and we know that the surface current is defined as 





b) Find the magnetic field B inside the capacitor by two different ways: using the integral Am- 
pere/Faraday law for two different loops such that one has a flux only of real current and the other only of 
displacement current. From your expression for B (or by other means) determine the vector potential A 
inside the capacitor. 


We know that Maxwell ’s equation is given by 
9E 
VxB — uoJ -- uo£0- 
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and the integral Ampere/Faraday law can be written as 


9E 
VxB= js.a- ] i: d2 oo [ 5: da 
A A Ot 
the Amperian loops used are given in the following diagram 


loop with real current loop with displacement current (inside the plates) 


ff OG Y» 
KG) j / 
A ---- 








For finding the magnetic field using the real current we can use the diagram on the left, we can see that 


since the loop is parallel with 
9E 


ae 


Thus we have 
fB di — [wos -da = uoles — uoK(x)l 
A 


thus we find that the magnetic field is given by 


£ok 
pcc (a—x)$ B=0 outside of the plate no z component) 
s 


now to find this using the displacement current we know that 
J =0 no surface current enclosed 


thus we find 





OE €ok 
$ B-ai= poco f Z -da = "(a1 — x2) (y1 — 72) 
A Ot s 


where we can also write the left hand side as 
o£ok E 
B(x1) ~B(x2) = =P (x1 —2)9 
where this can also be written as 


£ok 
p. "080 





(a—x)¥ x =x x-a 


which is the same thing we found from the real current, which it should be. To find what the vector 
potential is we can use 





& $ 2 
2 od 9 9 Y,  uo£ok A 
A, Ay A, 


Thus we can see that a suitable solution of the vector potential could be given as 


Eur 1, 


A, = —-)2 A,-0 
z -lax ; X M 
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Thus the magnetic field inside the plate and the vector potential can be written as 





c) Find B and A outside the capacitor. 


For part be we concluded that 


B=0 


for the vector potential outside we can guess a solution and see if it satisfies the boundary conditions we 
can use 


£ok x £ok 1 loa, 
Aabove = -A (az - zx) — PO" (ax z y* T 3:8 
does this satisfy our boundary conditions? 
VxA=0 
we can see that this is satisfied, 
V-A=0 


this is satisfied if we assume that we are using the Coulomb gauge which is defined as this. What about 


Aabove(Z = 0) — Apelow(z = 9) 
we can see that this is also satisfied. And finally 


dA above (Z = 0) E 9Abeloy(z — 0) 


on on 


and we can see that this will also be satisfied by our vector potential. Thus we can conclude that 


— —uoK 


€ok £ok 1 
B =0 A E99 (az — ax)& is (ax— =x? 


above above E 3 5* 





Problem £ 10 Magnetic pressure on a solenoid 
Consider a long solenoid that produces a magnetic field of magnitude B in its interior. 


a) Considering the magnetic energy per unit length in a solenoid of radius R and R + dR, show that 
there is an outward pressure on the surface of the solenoid (i.e. on the many turns of wire that resemble a 
continuous sheet of current) of magnitude P = B? / (2uo). 


We know that the energy stored due to a magnetic field is 


— 2uo an | Pe ms - 
where the volume of this solenoid is given by 
dV = xh|(R -- dR)? — R?] 
thus the energy per unit length is given by 


dW nB? 


qas R+) -R?] 
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We also know that the energy is related to the force by 


R--dR 
dW — [ 9a» [ aac 2mp f rdr = ThP|(R + dR)?” — R?] 
R 


and we can see that this can also be written as 


dW B? 
—— = Pn|(R+dR)” — R7] = —n|(R+dR) — R?] 
h 2uo 
thus we find that the pressure is given by 
B2 
P = — 
2uo 


b) Suppose a public utility proposes to moderate the ebb and flow of intermittent power sources by 
storing energy in a large inductor, reaching a maximum field of, say, 14 Tesla. What is the storage capacity 
of this inductor, in kilowatt-hours / (meter)? ? What would be the magnetic pressure on this solenoid in 
atmospheres? 


We know that 5 " 
B B 

W = —V x = — 

2uo V 24w 


and since we have been given the magnetic field and we know what uo is we find 


kWh 
W =7.8x ioi. =|21.7—— =770 atm 
V m? m? 


c) Alternately, suppose the energy is stored instead in a large capacitor with, say, a dielectric constant 
¢/€9 = 3 and a maximum allowable field of 1.5 x 105 V/m. This storage would not require, for instance, 
cryogenic cooling of a superconductor to avoid resistive heating of a current-carrying solenoid. What is 
the storage capacity of this capacitor, also in kilowatt- hours / (meter)? ? For reference, the energy content 
of gasoline is about 100 kilowatt-hours / (meter)? . 


We know that the energy stored in a capacitor with a dialectric media is given as 


£0 2 £0€£; 2 
W--—[eE'dv--——E-V 
a I 
thus 
W  eog,E? 
VvV 2 
We know that » 
Le £,—3 
£ 
thus w- 3 
— — -ggE? 
yu 


given that we know what E is we find 


W 24 kWh 


144 


Chapter 8 


Conservation Laws 


8.1 Problems and Solutions 


Problem # 1 


a) Consider two equal point charges q, seperated by a distance 2a. Construct the plane equidistant 
from the two charges. By integrating Maxwell's stress tensor over this plane, determine the force of one 
charge on the other. 








We can see that this configuration will cause a repulsive force on the top charge due to the bottom 
charge. If we pick a point located a distance r from the origin a distance r’ from the charge and at z = 0 
we see that the electric field here only has a component (the vertical components cancel by symmetry). 


E, — E, 20 Ey(r) —2Esin0 
Max well's stress tensor is given as 


1 1 1 
T;j = £0 CE m (ss - 58) (8.1) 
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The indices i and j refer to the coordinates x, y, and z, so the stress tensor has a total of nine components. 
We can write the stress tensor explicitily in matrix notation as 


z T Ty Ty; 
T- Ty Ty Ty: 
Tx Ty Tz 


One can form the dot product with a vectro a: 


(a. T); — Y dii 
i=x,y,Z 
If the force per unit volume is given as 
I oS 
f— V-.T-—£ouo— 
oUo or 
then the total force is given by 
= d 
F= f da - euo. | Sd (82) 
dt Jy 


where we used the divergence theorem to convert the first term to a surface integral. In static situations the 
second term drops out and we are left with 


F= f Tan (8.3) 
We can see from our figure that if we were to use the xy plane as our boundary surface then 
da — rdrdQ(2) 


We can write Maxwell's stress tensor as 
1 1 1 
Ta = zt(E; E, Ei) Ty — 5& (E, E;-E;) Tex = 5€ (E; -E;-E,)  — (84) 


Since we know the net force is in the Z dircetion it suffices to calculate 
(Í. da); = Tz,da, + T,yday + T,,da, 


since we know that 
day — day, — 0 


this becomes 
(T-da); — T;.da; (8.5) 
We can see that 


€ 
Ta =- E E, =E,=0 da, = rdrdọ 


Thus the force is given Equation 3 and Equation 4 as 
£0 
F = -2 f E?rdraġ 


we find that the electric field at r and 7' is given by 


E,—2Esin0 sin0—-. dignis 
y= = =z 
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The electric field for a uniform charge is given by 
q 1l. 
E(r)- ——— 
E greo r2" 
since we only have an electric field in the $ direction this becomes 


2q" Y 
— ———.sin0$ 
> Greg Y 
If we integrate over the entire xy plane and also a hemisphere extending out to infinity, where E = 0 at 


infinity and thus the integral contributes nothing allows us to write the force as 
2 
£o (2% f e ae cess 
p-2-— d —  —5sin0 | rd 
2 7*1. Vinea? ner 
2 pæ, 
£M c di 
4T£9 JO r'6 
From the diagram we can see that r’ = va? +r? thus 


This simplifies to 


ps ff ar = |e 
4n£o Jo (a?-4- r2)? 4ngp| 4(r?+a7)? |, 
Evaluating this integral gives 

2 





q ^ 
=| — 2 
16n1£9a? 
Which is exactly what we would expect using 
QU. m luis PM UE 
= 2 2a = 
4n£or? 161£9a? 


b) Do the same for charges that are opposite in sign. 






xy-plane 





positive charge and negative charge 
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We can see that for this configuration we find the net force is attractive and if we consider the top 
charge to be —g than it will experience a net force in the —2 direction. We can see that the electric field at 
r and r' is now given as 

E,—E,—0 E;—2Ecos0 


since we know that the net force is still in the Z direction then we can choose 
da = rdrdà(2) 


and thus we it works to solve for 
(T- da), =T,,da, 


where e 
T.— DU da; — rdrdà 
and so we find the force to be 
£o 2m P 2q 1 $ 
F = — d — — —. cos 0 d 
A Q a (aos rdr 


where 
2 


cos 0 = cos? 0 = a 


this now gives 





ga J ae 
= — dr 
4T£o Jo (@ +r?) 
doing a u substitutin gives 


d 
u=8@ +r rdr = = 


gives 








2,2 poo 2 
r= 74 / ud =|-—_* 9 
871£9 Ja? 16n1£09a 


which is exactly what is expected. 
Problem £ 2 


Consider an infinite parallel-plate capacitor, with the lower plate (at z — —d/2) carrying the charge 
density —o, and the upper plate (at z = +d/2) carrying the charge density +o. 
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a) Determine all nine elements of the stress tensor, in the region between the plates. Display your 
answer as a 3x3 matrix: 


Tx Try Tz 
T Tzy Tz, 


since the Maxwell stress equation is given by 
1 2 1 1, 
Tij - Eo (ss, — 59uE ) + i (so. = 38 ) 
and assuming that the magnetic field inside the capacitor is B = 0 we can write the nine components of 
the stress tensor as 


: EE, — jôu E? EXE, EE, 
T = £o E,E, EyEy — 58yyE* E,E, 
EE, EE; EE; — 4ôzE? 
this can be written as 
1 2 2 2 
E 3 (E E E?) - EE, E,E, 
T = & E,E, 5 (E, —E; — Ez) E,E, 
EEx E.E, } (E2 — E? — E?) 





b) Use Equation 8.22 


F= fT -da— eoi | Sdr 
S dt Jy 


to determine the force per unit area on the top plate. Compare Eq 2.51 


Given the static situation 


S=0 da=-z 
we find that the force is given by 
F o 1 0 0 ] 
f— Aa = oR —] O |f 
fh ree sov, Nc wp 


since we know the force is only in the -Z direction we find 


2 

= o 
f= T = ——z 
28g, 
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which is identical to Equation 2.51. 


c) What is the momentum per unit area, per unit time, crossing the xy plane (or any other plane parallel 
to that one, between the plates). 


According to Newton’s second law. the force on an object is equal to the rate of change of its momen- 
tum 


dp h 
F= mec 
dt 
where we can now write Equation 2 as 
dPmech gi d ik 
—— = ọ T- da — £ouo— | Sdt 8.6 
di i ouo), (8.6) 


The second integral represents momentum stored in the electromagnetic fields themselves: 


Pmech — A0£0 n Sdt 


while the first integral is the momentum per unit time flowing through the surface. Equation 6 is the 
general statement of conservation of momentum in electrodynamics. If we let Pyech be the density of 
mechanical momentum, and fem the density of momentum in the fields 


Tem -— uo£oS 


then equation 6 in differential form says 


d 5 
3; mech T)-VT 


and so we call -T the momentum flux density. We have found this to be 


o?, 


qe 2m 
29^ 


Pfux density 


d). At the plates this momentum is aborbed, and the plate recoils (unless there is some non-electrical 
force holding them in position). Find the recoil force per unit area on the top plate, and compare your 
answer to (b). [Note : This is not an additional force, but rather an alternative way of calculating the same 
force- in (b) we got it from the force law, and in (d) we get it by conservation of momentum.] 


The recoil force on the top plate should be equal and opposite to the force exerted on the plate 


Frecoil B _ 4Pmech _ 0. 


a dadi — Paux density ~ ~ 2€9 


which is the same as (b), like Griffith’s implied in the problem. Same problem just two different ways of 
looking at it. 


Problem # 3 
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A charged parallel-plate capacitor (with uniform electric field E = Ez) is placed in a uniform magnetic 
field B = BX as shown by the figure 


a) Find the electromagnetic momentum in the space between the plates. 


The electromagnetic momentum is defined as 


Pem = Ho£0 / Sdt = £o | (E x B)ac 
V 


thus we find 
pem — £oEBsinOV$ 


and since 0 — x/2 we find that the momentum is given as 
Pem = €gEBdAy 


V=dA 


where the volume is given by 


b) Now a resistive wire is connected between the plates, along the z axis, so that the capacitor slowely 
discharges. The current through the wire will experience a magnetic force; what is the total impulse 
delivered to the system, during the discharge? 


1- f rai 


The magnetic force on a segment of current-carrying wire is given by 


We know that the impulse is given as 
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thus the impulse can now be written as 


I= f [ 17axma- [ raxma 


since we know that | = dz and B = BX gives us 
I= f 1dB(2x &)dt — |-dBgit|g — —dB(q(s») — q(0)$ — q(0)aB$ 
0 


where q(t) — It and q(»») — 0, we assume that there is no charge left at t = oo. We also know that Gauss's 
law gives 
[9:4 dene. AP qu 
£0 
Thus we find that the impulse is given by 


I = e9EBdA¥ 
dpem 
I= | Fdt= dt = 
/ i I( dt ) Pem 


c) Instead of turning off the electric field (as in (b)), suppose we slowely reduce the magnetic field. 
This will induce a Faraday electric field, which in turn exerts a force on the plates. Show that the total 
impulse is (again) equal to the momentum originally stored in the fields. 


which is the same thing as 








Jj ux 


———*- 


E(z)y. hat 






| 77 — —  Amperean loop 


Since we know that Faraday's law in integral form is given as 
B 
j E-dl=— = -da 
ot 
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and in differential form 2B 
VxE=-— 
* ot 

we know that a changing magnetic field in the $ direction induces a Faraday electric field in the $ direction. 


Using the Amperean loop shown in the figure we find 


b d b d 
jEa-[ Edi. f E-a- f E-a- f E-dl 
a 0 a 0 


we can see that the second and last term cancel and since the electric field is in the $ direction we find 


(E(0) 4- E(d)) — 2 (B=) — B(0))dl = ABMS 





and the force is simply 





"ys 


F — qE — q(E(0) 4 E(4)) -«( 


using this in the impulse equation yields 
oB(0 
I= n = qd | a ) di$ — qdBg 


q = €0AE 


I= &EBdAĵ 


which is the same thing as in the momentum stored in the fields. 


Problem # 4 


but we have shown that 


thus the impulse is 


Imagine an iron sphere of radius R that carries a charge Q and a uniform magnetization M = Mz. The 
sphere is initially at rest. 


a) Compute the angular momentum stored in the electromagnetic fields. 


Since we know that the magnetic field outside of a uniformly magnetized sphere is given by the mag- 


netic dipole equation 
uom 


Bgip — dnm [2cos Of + sin 66] 
r 


and we know that the electric field inside this conducter is zero so we are only concerned with the electric 
field outside of this sphere which can be modeled as the electric field of a point charge at a large distance 


away 
qd . 


— f 
4n£or? 


the angular momentum density is given as 
lem =r X fem — £o[r x (E x B)] 
where we have defined the momentum density 


fem — uo£oS — £o(E x B) 
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thus we find 
uoqm 


1612£or? 
the total angular momentum stored in the fields is given by 


ExB= sin 0 


Lem = femV 


where 





we know that 0 in cartesian coordinates is 


^ 


© = cos O cos 9$ -- cosOsin 9$ — sin 0Z 


thus the angular momentum density is 


Hogm 


lem = -Ter sinð (cos cos ox + cos 8 sin oy — sin 8z) 


and the total angular momentum is 


> Hoqm ^ 1 2n T -3 ^ 
Yes i, lemar = ET | zdr [ dà [ sin? 6402 





which becomes 


ias 4(2n)uoqm 
(3)16:?R 
and if 
m= AM 
3 
gives 


Lem — -uoMqR?2 


i 


b) Suppose the sphere is gradually (and uniformly) demagnetized (perhaps by heating it up past the 
Curie point). Use Faraday's law to determine the induced electric field, find the torque this field exerts, 
and calculate the total angualr momentum imparted to the sphere in the course of the demagnetization. 


s-rsin 0 





If we let 
dl=—2nsh da=ns72 


fE a= J-Z aa 


thus the Faraday electric field is given as 


thus we can use Faraday’s law 


2 
Efar = -Zr = -im R sin0 
but we know that 
gnia a 
MEE. ddp 2 3 1 di 
the torque is given as 
N=rxF 


where the force is from the induced Faraday electric field 
F = qEfar. dE = dqEfar = OdaEfar 


and the total torque is given by 
N 
f aN = | rxdF= f RAF(® xô) s —RdF (6 ) = | tar sinez 
0 
da — R?sin0d0dó 


1 M _ 1264nR* M 
=55R e no [ ag |” sin? 0402 — icu m 


if we now define 


thus 





thus the total torque is given d: 
2 dM 
N = ZugqR? —42 
904^ m 
and thus the total angular momentum is given as 


2 " 
Lax Nat g dttoR Mi 


€) Suppose instead of demagnetizing the sphere we discharge it, by connecting a grounding wire to 
the north pole. assume the current flows over the surface in such a way that the charge density remains 
uniform. Use the Lorentz force law to determine the torque on the sphere, and calculate the total angular 
momentum imparted to the sphere in the course of the discharge. (The magnetic field is discontinuous at 


the surface... does this matter?) 


We know that the force will now be given by 
F = q(v x B) — (K x B)da 
and we know that the torque is now given by 


N=rxF=(rx(KxB))da 
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using the BAC-CAB rule for the cross product of a cross product we find 
N = (K(r-B) —B(r-K))da 
we know that r- K = 0 due to orthogonality, thus we find 
N=K(r-B)da 
we know that the magnetic field will just be the average of the internal field and the external field 


1 | uom ^. nm Hom cos OR? . uomsin OR? » 
B = - | =— = eee 
5 | an (cos 0f — sin 00) au Ane 0 


and 


pom 
r-B = ——cos0 
2nr? 


now we need to find what the surface current is, we know that 
2x T 
q= foda = of ao | R?sin0d0 — o2xR?(1-- cos0) 
0 0 


and 





dq d 2 
— — —DmnR*o(l 0 
Ji 3 xR^o(1--cos9)] 
we also know that Q 
rius 4nR? 


thus we find the current to be given by 


I= dQ (1+cos® 
dt 2 


but we know that the surface current is given by 
I 


K=—6 
rae 
but 


1, =2nRsin@ J=—— 2 


dQ (1+cos® 
dt 


thus the surface current is 


_ dQ ] 4- cos 0 6 
dt \ 4nRsin® 


thus our total torque is given by 





_ dQ ] -- cos0 A /uom A 
ME c) Ca 29) dab 
but we know that 
da = R?sinOd0dó — 0 — (—sin02) 
thus this integral be 
d 20 T 
N= TE : ao | (sin® cos @ + sin®@ cos” 6)d0z 
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since we know that 
m M 
AnR3 3 
thus we find 


M 2x T 
N = HOM dO do [ (sin@ cos @ + sin® cos” 0)d6z 
67 dt Jo 0 


which gives 


L. pude. t. ; 2 
N = -uR u= | (sin cos + sinOcos* 6)d6 
3 dt Jo 
which becomes à do 
N = ZugMR? — 
gem 


and thus the total angular momentum is 


2 2 
Lren NiE g'oMR Q2 
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Chapter 9 


Electromagnetic Waves 


9.1 Problems and Solutions 


Problem # 1 
Suppose 


inĝ . ^ . 0 
E(/,0,6,1) = A~— |[cos(kr — wt) — (1/kr) sin(kr —ot)]6, with P7 
r 
(This is, incidentally, the simplest possible spherical wave. For notational convenience, let (kr — œt) = u 


in your calculations.) 


a) Show that E obeys all four of Maxwell's equations, in vacuum, and find the associated magnetic 
fields. 


In regions of space where there is no charge or current, Maxwell's equations read 


(i) V-E = 0 (iii) vxg--2 
(ii) V-B = 0 (iv) V x B= peo 


for (i) we find 
V-E= ye ese: (1/kr)sinu]ó 
r 


in spherical coordinates this gives 


A Q 
V-E= gap n — (1/kr) sinu] = [0] 
to find the magnetic field we must use 
oB 
VxE=-— 
5 ot 


where letting 


Vo iy andl er ee sinu] 
r 
1 di ~ lð 
VxE= Fane [ag sine r— ap lvoe 
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which is given by 








A2sin0cos0 1 Asin® d 1 
VxE = [eos — sinu]? — m 4; eosu — tL sinu] 
A2cos0 1 Asinð 1 
VxE = eee T cosu — — sinu]f — ea —ksinu + — sinu — kcosu 6 
r? kr kr 


we know that this is equal to 


1 
—ksinu + ga ne — kcos 





cosu — — sin ulf — 


A2 cos 0 1 . AsinO 
r? | kr | 


2 9B 
i ua Ot 


thus the magnetic field is given by 


A2cos0 1 AsinO 1 1 be 
B= | > [cosu — z sinu]êdi — f zd |- ksinu + — sinu — Leos! Odt 
r kr r 








kr2 
solving this integral gives 


2A cos 8 


. 1 ~ AsinO 
5 sin u + — cosu| f+ 
or kr 


1 1 
|- —keosu- L5 cosu-t - Zr 





now for the divergence of the magnetic field 
1 ð / ,2Acos0 1 
V.B = Ay G TOS b inut cosu] ) 


] d (* sin? 0 





1 1 
TE a [—keosu+ oy cosu += sina ) 


which is simply 0 
V.B—0 


now to for the last of Maxwell's equation we use 





oE 
V x B = uo£0 — 
Uoto àt 
if we call 
2A cos 0 1 AsinO 1 1. 
v = ——=— sinu + —cosu| ve = —kcosu + — cosu + -sinu 
or kr or kr r 


Thus the cross product gives 





1 d Qv, A 
V xB = - |> (rva) - — 
SB r ED x 
2 j L1 5889 k? sinu — I sinu — z COS Uu — ! in osi 
or a o fo kr? r? r 
ang dv, 2A 1 
Vr a. . uu em. id 
d ae inus E cosu! (—sin0)K 
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thus summing up these two contributions yields 








Asin®@ 1 2 9E 
yxp [ksinu+ Ecosu $ = meo 
r ot 
now we know that 
Æ 102 in ~ lAsin0|o ol Asin®@ 1 
moy ay Be kosu (1/kr)sinu| $= 3e [Pasina+ 2 cos] = = [sinu + eosu 


thus we can see that both sides are equal. 


b) Calculate the Poynting vector. Average S over a full cycle to get the intensity vector I. (Does it point 
in the expected direction? does it fall like r^, as it should?) 


We know that the Poynting vector is given by 


s- l(ExB) 
Ho 


given the magnetic field found in part (a) and using the electric field given gives (minus a lot of algebra) 


f 6 Ô 
Sh 0 0 A828 [cos(kr — wr) — (1/kr) sin(kr — ot) 
nm 2Acost [cosu -- i: sinu] sae |-kcosu + z cosu+ Esinul 6 0 
which will give 
A? sin” 0 1 1 
S = -r leosu— (1/kr) sinu] -kosu + pgcosu+ sinu! f 
2A? sinOcos@ T A 
Ju. POSES (1/kr)sinu] osu sina] 0 


Simplifying the f component gives 


A^sin9? [/2 1 \sindu 1 x 


we know that the intensity is given by 
I= (S) 
but we also know that 
1 
(sinucosu) = (sin2u) = (cos2u) — 0  (cos?u) — (sin? u) = 5 


we find that the intensity will now be given by 


A?ksin9? 
I LI S == —— f 
( ) 2ugor? 5 


c) Integrate I- da over s spherical surface to determine the total power radiated [Answer : 4TA? /3uoc] 
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we know that 





A?k 2h 4A?2mk _ 4nA? 
Pod e d f 39 — = 
f a d PSAE eadk 


Find the fields, and charge and current distributions, corrsponding to 


Problem # 2 


1 qt 
V(r,£) 20, A(r,t) = -——-—f 
(r, ) , (r, ) 4T£0 Ar 
We know that 3À 
q 4 
E-—--VV——z-|——s 
ot Ameor 
and the magnetic field is given by 
B=VxA=[0] 


to find the charge distribution we can take the divergence of the electric field 


r 
p=€0V-E= $v. (3) = q& (r) 


and finally the current density 


1 9E 
J= zV xB- ez; = [0] 


since we know that the magnetic field is O and also that the electric field is constant in time. 
Problem # 3 


Use the gauge function A = —(1/4reo)(qt/r) to transform the potentials in Prob. 10.3 and comment 
on the result. 


Using Equation 10.7 we find 


A'=A4+VA=0 Vi= y .9* 
ot 
thus 


B — V x A' — [0] 


ðA’ q 
E=-VV'— =| — fî 
dt  |4n£or? 


J=[0] 


These are the same conclusions, this follows from the fact that we can add any VÀ to A, provided we 
simultaneously subtract JA /dðt from V. 
Problem # 4 


we can also see that 


and the current distribution is again 


Find all the elements of the Maxwell stress tensor for a monochromatic plane wave traveling in the z 
direction and linearly polarized in the x direction 


~ po i 2 1 z : 
E(z,t) = Ege g B(z,t) = z Eod 9 
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or taking the real part 
1 
E(z,1) — Eocos(kz — ot 4- 6)& B(z,t) — —Eocos(kz — ox -- 6)y 
c 
Does your answer make sense? (Remember that T represents the momentum flux densitu.) How is the 
momentum flux density related to the energy density in this case? 
1 2 1 1 2 
Tij = £0 EjE; — -ÓjjE + — BiBj — -óijB (9.1) 
2 Ho 2 


The indices i and j refer to the coordinates x, y, and z, so the stress tensor has a total of nine components. 
We can write the stress tensor explicitily in matrix notation as 


= Ty Try Ty; 
T- Ty, Ty, T. 
TX Try Tz 
We can write Maxwell’s stress tensor as 
1 1 
Tem ;&o(E? — Ej — E?) tas 55 -BO (9.2) 
1 1 
Des 5€» (Ey — EL - ED) t 57 (By - BE - B5) (9.3) 
1 1 
T. = 580 (E; —EL-ED- E - B? — B°) (9.4) 


and 








and since we know 


1 2 1,5 
Ta = 30k: = Dae? (9.5) 
l mz, lm 
M NEST 
T = — 580k = 35 (9.7) 
thus Maxwell’s stress tensor in matrix notaion is given by 
lo p2- 1 p2 
T puc my 1 : 1g i 
T= 0 —5€oE; zu P» l i sd 


where the off diagonal elements are all zero. Plugging in our expressions for the electric field and magnetic 


field gives 
0 


2 0 0 
T—-/|00 0 
0 0 —esoEe cos (kz — ot +ò) 
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this makes sense because the wavefront is moving in the z direction. In this case the momentum density is 
equal to the negative of the momentum density. Griffiths defines the energy density on Equation 9.55 as 


u — €9Eg cos? (kz — at 4- 8) 


which is 
u = —Ī = —T; = £0Eĝ cos? (kz — œt + 8) 
Problem # 5 


a) Suppose you imbedded some free charge in a piece of glass. About how long would it take for the 
charge to flow to the surface? 


The equation that governs how charge dissipates through a conductor is given by Equation 9.120 as 


py — e 9/9'p,(0) 


thus any initial free charge densiy p (0) dissipates in a characteristics time 


a 
II 
alm 


where 1 
= P glass ~ (10!° — 10H) 





Oglass 
where P giass 18 the resitivity and Ogjass is the conductivity and the permitivity is 


Eglass © 4.9 — 7.580 


we find that the characteristic time is given by 


0.433 sec ~ 109 Egiass © 4.980 


a 
Q 





Oglass 


6198 sec 


a 
Q 





104 estass 7.080 


glass 
b) Silver is an excellent conductor but it's expensive. suppose you were designing a microwave exper- 
iment to operate at a frequency of 10!° Hz. How thick would you make the silver coatings? 


We know that the distance it takes to reduce the amplitude of a wave by a factor of 1/e (about a third) 
is called the skin depth which is defined as 


(9.8) 
where K is defined as 1/2 
EET ns 
«sey [re | (9.9) 


thus if the thickness of our silver sheet is given by 8 we will have a coating with a thicknes that will stop 
the microwaves by at least 1/3. Given that 








Usilver ~ HO = 1.59 x 1078 Esilver = £0 


silver 
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since we also know that @ = 27f where f is the frequency we find 
ke T5710 m^ 


thus we should make this 
d 726.36 x 1077m 


we should probably make it a little thicker than this to stop all EM waves, i.e 


width z 2d z 1.27 x 109m 


c) Find the wavelength and propgation speed in copper for radio waves at 1 Mhz. Compare the corre- 
sponding values in air (or vacuum). 


We know that the wavelength and speed of propgation is given by 











2x o 
A= — ae 
k ^ k 
where k is defined as T 
£ o2 
C= ete (—) +1 
2 £Q 
for copper 
1 
= Pcopper © 1.68 x 1078 Q — 2nf =21 x 10657! 
Ocopper 
thus 


Ocopper © 5.9 x 107 
we also know that €copper © (2 — 20)€o 
0 £0 

we can know that 

o2 0 \2 

(—) siS 1+(=) >] 

EO € 

and thus we can write 


eu |o oou 
k z: Q4 | — A — F 4| — 
2Ve 2 
where u © Ug we find 
60 
ke ALS S 1.5 x 10m7! 





thus the wavelength is 


2 
es T AAT x 1074m 


and the velocity is given as 


v= — x 419 m/s7! 


=le 


Problem # 6 


A plane wave travels in the z direction in a conductor with real conductivity. 
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a) Find the instantaneous and time-averaged power loss per unit volume due to resistive heating for 
any z. 


The electric field and magnetic fields in this case are given by Equation 9.138 in Griffiths as 


E(z,t) = Ege "*'cos(kz — Qt J- Og)& (9.10) 
B(z:) — Boe "cos(kz — Qt -- óg 4- 0) 


The power loss due to resistive heating is given as 


dP — d(VI) 
we know that 

dV — Edx 
and 

dI =J - da = Jdydz 
thus 
dP = JEdxdydz 

where 


J=6E 


thus the power loss per unit volume is given by 

P 

wc GE? = oEje 7 cos” (kz — wt + 8.) 
and the time average power loss is given as 


P —2Kz p2 
zc E 
y 2 


because we know that 
(cos*(kz — ot + 8z)) = = 


b) Find the total power loss per unit area between z = Qand z — oo. 


To find the total power loss per unit area we must integrate the power per volume 


= p 2672 


c) Find the time-averaged Poynting vector at any z. 


The Poynting vector is defined as 


1 
S — - (Ex B) 
u 
using Equation 10 and Equation 11 we find 
BoE, 
S= 70 eK: cos(Kz — ot + dg) cos(kz — ot + dg 4- 9)]Z 
u 
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if we let x = kz — @t + dg and 6 = 77/4 we find 


BoE 
S = 0 Oe 2Kz 


BoE 
0 Pee 2 
u 


cos(x) cos(x 4- 0)]Z cos^ xcos 9 — cosxsinxsin $|Z 


The intensity is given as 


F COS OZ 2 
since we know that 
tand=— cosd= p. 
A / k2 -++ K2 
and since we know that 
K=vVk +K? 
we find 
K coso = k 
thus the intensity is given as 
k Ej os, 
I= (S) = ——e “2 
Qu 2 


d) Compare the value of your result for parts (b) and (c). Are the results in accord with one another? 
Explain. 


at z — 0 the magnitude of these two functions should be the same. For the Poynting Vector we find 


_ KBR 


S(z = 0) 2o 


and for the energy per unit area per unti volume we find 











A 4k ° 
we need to show 
k o 
uo 2K 
this can be done using 
1/2 5 1/2 
& 
kx — o2 8 Le (—) +1 1+(=) ad ege 
2 £0 £0 2 £0 
thus we find 
kK o 
LO Dig 


which is what we were looking for all along thus these two quantities are equal at z = 0 
Problem # 7 A parallel-plate capacitor at high frequencies 


When one designs an electronics circuit with conventional devices like resistors, capacitors and in- 
ductors, one is wise to keep in mind that “everything changes at high frequencies”. To illustrate this, let 
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us consider a capacitor composed of two closely spaced and parallel circular plates, with radius a and 
separation h (h < a). 


a) Consider first the “quasi-static” situation, by which we mean the low-frequency response of the 
capacitor. If the capacitor is driven with a single-frequency sinusoidal voltage V = Voe’™ , what are the 
electric and magnetic fields, E(r,t) and B(r,t), produced inside the capacitor? You can write these as 
complex quantities, but keep in mind that the physical values of fields, etc., are given by the real parts of 
the complex functions. 


We know that 


N> 





= Mecos wr + isin Or )Z 
thus the real part is given by 

E(r,t) = Meos or 
we also know from Maxwell’s equation 


19E 
c? ot 


l 
= [E-da= pB-al 
C 


da—mnr dl-2mró 


—VxB 


and using Stoke's theorem 


we know 


thus we have 


nr? dE 
T B2xr 
where 
dE Vo. 
o —,9sin cor 
thus we find the magntic field to be given as 
r ` ^ 
B(r,t) = Oo a sin wrod 


b) From the formula relating the magnetic energy and self-inductance of an inductor and from your 
answer above, what is the self-inductance for this capacitor? In determining this energy, you will have to 
consider a specific volume over which to integrate the energy density - take this to be the inner volume of 
the capacitor. 


We know that the energy is given by 


1 
W = -LP 
2 
Equation 7.34 gives the energy of the magnetic field as 


1 
w=] Ba 
2uo all space 


thus 
1 


1 
= | Bdt = -LÊ 
2uo all space 2 
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the left hand side becomes 


1 2 
x. i (ows sin on) do diccns 


which gives 


1 b : 2 hn OVo 2 sc a 3 
no) (ovoz zsinor) 2Thrdr = T (se | sin o f r-dr 





thus j 3 
Ta OVo a) 
=—— | — cor 
Aoh (23) v 
to find the inductance z 5 
— | — ot = —LI 9.11 
Auoh ( 2c ) "A 2 Men 
if we let > 
dV Ta^9 dV Vo. 
I=C— C= —— — = —0— sin Ot 
T h dt Be 
thus 2 
2 
£ 
P= (5 2T sinar) 


substituting this into Equation 15 gives a inductance of 


uoh 


8T 


c) From the standard capacitance C and the inductance L that you just calculated, determine a resonant 
frequency 0 for the capacitor as you would for a standard L — C circuit. The frequency 0 is a good 
benchmark for what qualifies as “low” or quasi-static vs. “high” frequency operation of the capacitor. 


We know that for n L-C circuit the angular frequency is defined as 


1 
o = —— 
VLC 
we know that 
ud 
LC = — 
8c? 
thus the resonant frequency is given by 
a 


d) In (a), you will have determined the magnetic field generated by the time-varying electric field from 
Faraday’s law. But, of course, the time-varying magnetic field also generates an electric field according to 
the Ampere-Maxwell law. Taking the magnetic field discussed in parts (a) and (b), calculate this correction 
to the electric field of a capacitor. 


The Ampere-Maxwell law gives 


using Stoke’s theorem we find 
d 
E-dl=—< [Beda 
f dt 


da = hdr 


where 


and using the magnetic field from part a) we find 








E(0)h — E(r S sf oor a z Sin wt (hdr) 
where 5 
oV 
3: sf QVo7 2 z Sin @t (hdr) = acm r? cos Qt 
defining 
E(0) =0 
we find 
2 
oV 
(r)2— Tha r^ cos Mrz 


this is the correction to the electric field due to the changing magnetic field. Thus the total electric field is 


Vo e^ Vo » 
Etot = ($- Anz COs MLZ 





e) You should find that, at this level of approximation, the electric field inside the capacitor is parabolic 
as a function of radius r. At what frequency @ does the electric field at the outer edge of the capacitor a 
go to zero? How does it compare with the frequency € determined above? 


If we set the electric field equal to zero we find 


Vo oV, 5 


= 
h^ Am 





we find that the frequency that makes this go to zero is 


which is smaller by V2. 


Beyond the approximation discussed above, a complete solution for the electric and magnetic fields 
inside a capacitor can be found that is consistent with MaxwellOs equations. As often in problems with 
cylindrical symmetry, the solutions for the spatial variations of E and B are Bessel functions. These 
functions also serve to characterize the resonant modes inside a cylindrical metal can. 
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